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SECTION  I 


INTRODUCTION  AND  SUMMARY 


1.1  BACKGROUND. 

In  the  early  1980s  the  Defense  Nuclear  Agency  began  to  investigate  the  long 
wavelength  infrared  backgrounds  expected  to  result  from  high-altitude  nuclear 
detonations.  The  gamut  of  activities  included  theoretical  efforts,  laboratory 
measurements  of  spectra  and  fluorescence  efficiency,  predictive  computer  codes, 
and  research  toward  an  upper  atmospheric  uranium-release  simulation  in  October 
1982,  H.  C.  Fitz,  Jr.  [Chief,  Atmospheric  Effects  Division  (RAAE)]  established  a  six- 
member  working  group  [Uranium  (Oxides)  LWIR  Review  Committee,  chaired  by  i.  L 
Kofsky  of  PhotoMetrics]  to  formally  review  all  the  technical  aspects  of  that  LWIR 
program.  One  major  function  of  that  committee  was  to  consider  the  effectiveness 
of  field  measurements  of  the  LWIR  from  uranium  oxides  produced  by  controlled 
releases  of  uranium  vapor  from  rockets  in  the  thermosphere  Reports  of  the 
committee  reviews  were  provided  in  detailed  memoranda  by  the  chairperson  to  the 
RAAE  program  manager  An  excellent  overview  of  the  committee's  activities  and 
the  then-current  status  of  the  uranium  oxides  radiation  program  was  prepared  by 
Kofsky  [Ref.  40;  Section  2  &  Appendix],  More  recently,  a  report  [Ref.  60]  by 
Visidyne,  Inc.  personnel  describes  efforts  directed  principally  to  the  design  of  an 
experiment  for  the  release  of  uranium  vapor  in  the  upper  atmosphere  in  sufficient 
quantity  to  permit  an  evaluation  of  the  impact  of  high-altitude  nuclear  weapon 
debris  emissions  on  advanced  space  systems. 

1.2  SUMMARY. 

In  addition  to  the  formal  meetings  of  the  Uranium  (Oxide)  LWIR  Review 
Committee,  the  writer's  more  extended  support  of  the  committee  was  the 
preparation  of  numerous  informal  working  memoranda  on  relatively  narrow-in- 
scope  technical  issues  which  had  been  raised,  Those  selected  topics  involved 
aspects  of  neutral-neutral  and  ion-neutral  atomic  and  molecular  collisions  for  u- 
release  conditions,  i.e.,  velocity-dependent  momentum-transfer  cross-sections, 
stopping  power  and  range  versus  energy,  scattering  angles,  diffusion  coefficients, 
and  reaction  rate  coefficients.  This  information  would  ultimately  be  used  to  make 
improved  estimates  of  the  geometrical  extent  of  the  LWIR  radiating  region  for  U- 
release  conditions.  Those  memoranda  have  been  extensively  reorganized  and 
synthesized  into  the  body  of  this  report. 

in  Section  2  we  consider  the  topic  which  is  basic  to  the  other  sections,  i  e  , 
momentum-transfer  cross-sections.  Initially  our  general  approach  was  (a)  to 
consider  some  general  formulas  for  atom-atom  and  ion-atom  collisions,  (b)  to 
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apply  and  extend  those  formulas  to  [(0,0+MO,N2,02)]-coiiisions,  since  there  is 
more  known  about  such  collisions  than  those  for  uranium  species,  and  (c)  to 
consider  procedures  for  uranium-species  collisions 

More  specifically,  for  atom-atom  collisions  at  low  velocities,  we  obtain 
velocity-dependent  momentum-transfer  cross-sections  from  the  work  of 
Hirshfelder,  Bird,  and  Spatz  [Ref.  3!]  by  fitting  the  Lennard- Jones  ( 1 2,6)-potent ia! 
to  the  ground-state  potential  of  the  molecule  corresponding  to  the  colliding  atoms. 
This  cross  section  is  then  extended  to  lower  relative  velocities  by  using  the 
asymptotic  results  of  Monchick  and  Mason  [Ref.  59]  and  of  Klhara,  Taylor,  and 
Hirshfelder  [Ref.  38]  for  an  attractive  (r6)-potential.  Thus,  there  are  procedures 
available  for  computing  momentum-transfer  cross-sections  for  (O+O)-  and  (U+0)- 
collisions  based  on  the  (r6)-potential 

The  situation  is  more  complicated  for  atom-molecule  collisions.  Our 
approach,  for  lack  of  a  recognized  simple  alternative,  is  to  consider  collision- 
integrals,  which  are  momentum-transfer  cross-sections  averaged  over  the 
relative  velocity  distributions  in  a  Maxwellian  gas  at  temperature  T.  We  appeal 
to  collision-integral  results  in  the  literature  for  [0+(0,N2,02)]-collisions.  from 
which  we  infer  an  approximate  scaling  factor  to  relate  [0+(N2,02)]-collisions  to 
(0+0)-collisions.  That  same  scaling  factor  is  then  applied  to  (U+0)-collisions  to 
get  cross  sections  for  [U+(N2,02)]-collision$. 

For  ion-atom  collisions,  at  sufficiently  low  velocities  where  (at  least  for 
some  collision  pairs)  the  polarization  force  dominates,  the  momentum-transfer 
cross-section  is  based  on  the  polarization  potential  (r~4).  Thus,  a  procedure, 
which  seems  to  be  on  a  relatively  firm  basis,  is  available  for  computing 
momentum-transfer  cross-sections  for  (0++0)-collisions  . 

However,  for  (U+*0)-collisions,  we  deduce  evidence  that  the 
polarization  potential  appears  to  be  inappropriate,  owing  to  its  smallness 
relative  to  the  (r~6)-potential  at  the  "orbiting  radius."  Perhaps  the  polarization 
potential  should  be  replaced,  in  principle,  by  the  ( !2,6,4)-potential  or,  in  practice, 
by  the  (12,6)-potentlal  (fitted  to  a  ground-state  U0+  potential).  As  for  atom-atom 
collisions,  such  a  cross  section  can  be  extended  to  lower  relative  velocities  by 
using  the  mentioned  asymptotic  results. 

To  proceed  with  [U*+(0,N2, (^-collisions,  in  the  middle  of  our  U-release 
studies,  we  radically  shifted  our  approach  by  detouring  to  study  Ba+  collisions 
with  atmospheric  species,  hoping  for  useful  insights.  Indeed,  that  study, 
documented  in  Appendices  a,  b,  C,  and  D,  revealed  the  following  very  import anr 
facts  (better  appreciated  after  reading  Appendix  C)  which  were  basic  to  the 
writer's  subsequent  U-release  studies 
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a.  The  experimental  mobility  of  Ba*  In  N2  (at  room  conditions)  agrees  well 
with  the  theoretical  value  based  on  the  polarization  (r_4)-potential. 
b  The  (r_6)-contribution  dominates  the  (r~4)-contribution  in  the  ( 1 2,6,4)- 
potential  for  (Ba++0)-co1lisions  (analogous  to  the  situation  for  ((U,U+K)]- 
collisions) 

c.  An  effective  value  of  the  parameter  (/Jr,)w  -  which  appears  in  the 
asymptotic,  low-energy  momentum-transfer  cross-section  based  on  the 
(r&)-potentiai  -  is  postulated  to  be  determined  by  requiring  the  equality 
of  mobilities  (at  STP  conditions)  based  on  an  effective  (r6)-potential  and 
on  the  (r_4)-potential.  (At  non-STP  conditions,  the  two  mobilities  will 
differ.) 

d.  values  of  the  ambipolar  diffusion  coefficient  [/},  E(#  ]  computed  from 

an  effective  (r-6)-potential  are  -  for  ionospheric  ion  temperatures  ranging 
from  11 35  K  to  750  K  -  about  79%  to  85%  of  the  ambipolar  diffusion 
coefficient  [4  p*.]  based  on  the  polarization  potential  An  overall 

comparison  of  Dn  E(f  (r_*)  pot  with  (sometimes  ambiguous  or  uncertain) 

experimental  values  [z\  exp  |  derived  from  12  Ba-releases  is  about  the  same 
as  for  po,  ^ ,  so  that  a  strong  preference  cannot  be  established  on  this 
basis. 

Not  having  the  experimental  mobility  of  U+  in  N2  (or  any  gas),  we 
qualitatively  extrapolate  from  (Ba++N2)-collisions  to  [UM0,N2, (^-collisions 
and  postulate  that  (a)  the  mobilities  of  U*  in  atmospheric  species  at  STP 
conditions  can  be  predicted  by  using  the  (r-4)-potential  and  (b)  the  effective 
values  of  the  collision  parameter  (rf,e0)vz  can  be  determined  by  requiring  equality 
of  the  mobilities  (at  STP  conditions)  based  on  an  effective  (r~6)-potential  and  on 
the  (r_4)-potential.  These  postulates  allow  us  to  use  all  the  general  formulas 
developed  in  the  Appendices  for  Ba+  (and  Ba)  collisions  and  apply  them  to  U"  (and 
U)  collisions.  Those  applications  are  documented  in  the  remaining  portions  of 
Section  2. 

In  Section  3  the  overall  goal  is  to  obtain  simple  range-energy  relations  for 
low-energy  uranium  species  traversing  the  ambient  atmosphere.  Such  relations 
are  obtained  by  integrating  an  appropriate  stopping-power  formula  involving  mean 
values  of  the  momentum-transfer  cross-sections  developed  in  Section  2 

To  obtain  a  stopping-power  formula  for  a  single  particle  traversing  a 
Maxwellian  gas,  we  adapt  ideas  from  Banks  (Ref.  6]  who  was  interested  in  the 
related  problem  of  elastic  collisions  and  energy  transfer  between  gases  which 
have  separate  Maxwellian  velocity  distributions.  Banks  [Ref  6]  shows  that  the 
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exact  equation,  derived  by  Deslogue  [Ref.  17] ,  for  the  energy  transfer  rate  Detween 
gases  of  arbitrary  temperature  and  particle  mass,  can  be  separated  into  three 
fundamental  factors,  each  of  which  depends  on  a  different  aspect  of  the  collision 
process  and  gas  composition,  i.e.,  (Da  ratio  of  particle  masses,  (2)  the  difference 
in  the  gas  thermal  energies,  and  (3)  a  nonequtlibrium  collision  frequency  for 
energy  transfer  which  entails  a  mean  nonequilibrium  momentum-transfer  cross- 
section. 

Banks  [Ref  6]  acknowledges  that  his  decomposition  of  Deslogue's  formula 
[Ref.  17]  is  not  rigorous  since  he  did  not  use  the  proper  averaging  techniques  of 
kinetic  theory  needed  to  arrive  at  an  exact  expression.  Whatever  shortcoming  that 
fact  imposes  probably  is  not  very  important  relative  to  (a)  our  additional 
approximations  and  (b)  the  uncertainty  in  the  basic  cross  sections  used  in  our 
applications.  To  provide  insights  concerning  the  numerous  approximations 
required  to  obtain  simple  analytical  expressions,  we  first  sketch  Banks' 
development  in  general  terms,  followed  by  more  details  with  equations  and  our 
adapted  formulas,  as  given  in  Section  3 

In  Section  4  we  present  (a)  formulas  (developed  in  Appendices  A  and  D)  for 
the  ordinary  and  ambipolar  diffusion  coefficients,  based  on  both  the  polarization 
potential  and  the  effective  (r*)-potential  discussed  in  Section  2  and  (b) 
evaluation  of  the  diffusion  coefficients  at  four  altitudes  [200(25)275  km]  in  a 
quiet  MSIS-83  atmosphere. 

in  Section  5  we  compare  some  of  our  momentum-transfer  cross-sections 
developed  in  Section  2  with  cross  sections  for  certain  uranium-oxide  reactions, 
specially  atom-transfer  reactions.  (This  writer  expects  momentum-transfer 
cross-sections  should  exceed  atom-transfer  cross-sections.)  As  a  by-product  of 
these  considerations,  rate  coefficients  are  derived  for  selected  uranium  oxide 
reactions  for  which  velocity-dependent  cross-sections  are  either  known  or 
assumed,  based  on  extrapolations  and/or  averaged  values  of  measurements. 

As  previously  mentioned,  the  work  presented  in  Appendices  A,  B,  C,  and  D 
was  developed  during  the  auxiliary  study  of  Ba+  collisions.  Early  in  that  study  the 
writer  found  that  the  (ordinary)  ion  diffusion  coefficient  was  written  differently 
by  various  authors.  Thus,  a  first  step  was  to  reconcile  those  different 
expressions,  as  is  done  in  Appendix  A.  We  also  specialize  the  formulas  for  several 
parameters  [(a)  the  mean  momentum-transfer  cross-section,  (b)  the  collision 
time,  (c)  the  mobility,  and  (d)  the  diffusion  coefficient]  for  three  interaction 
potentials  [(a)  polarization  (r-4),  (b)  (r6),  and  (c)  (hybrid)  effective  (r6)].  These 
formulas  are  written  explicitly  for  a  binary  gas,  except  that  for  the  diffusion 
coefficient  which  is  valid  for  a  multicomponent  gas. 


4 


in  Appendix  B  we  elaborate  on  ion-neutral  collisions  in  a  multicomponent 
gas  by  (a)  commenting  on  literature  formulas  for  collision  time,  (b)  presenting 
explicit  formulas  based  on  both  the  polarization  potential  and  effective  (r-6)- 
potential,  valid  in  a  multicomponent  gas,  for  mobility,  collision  time,  and 
diffusion  coefficient,  and  (c)  comparing  the  Ba+  collision  times  computed  from 
various  formulas  for  a  number  of  altitudes  in  a  selected  atmosphere. 

In  Appendix  C  we  discuss  Ba  and  Ba+  collisions  with  atmospheric  species. 
(It  is  these  considerations  which  influenced  our  treatment  of  uranium  collisions  > 
For  (Ba*+0)-collisions,  we  have  two  simple  methods  for  estimating  momentum- 
transfer  cross-sections.  The  first,  based  on  a  polarization-potential  interaction, 
would  probably  be  the  natural  choice.  However,  estimates  indicate  that  the  r-6 
term  is  much  more  important  than  the  r4  term.  Thus,  we  are  led  to  seek  a 
compromise  between  these  two  conflicting  approaches.  Several  observations 
obtain: 

(a)  The  experimental  mobility  of  Ba+  in  N2  [Powell  and  Brata,  unpublished  [Ref. 
58]]  can  be  predicted  to  within  about  four  percent  by  use  of  the  polarization 
potential,  as  shown  in  detail. 

(b)  A  similar  situation  is  postulated  to  hold  for  (Ba+*0)-collisions 

(c)  There  are  reservations  about  the  applicability  of  the  polarization  potential 
at  the  higher  velocities. 

(d)  There  is  concern  about  the  (r~6)-potential  giving  a  larger  momentum- 
transfer  cross-section  for  (Ba+0)-col  1  isions  than  that  given  by  the 
polarization  potential  for  (Ba++0)-collis1ons  at  (room)  thermal  energies 

In  view  of  these  considerations  we  suggest  a  compromise  by  introducing  what  we 
will  call  an  effective  (r~6)-potential,  determined  by  requiring  an  equality  of 
mobilities  (at  STP  conditions)  based  on  (a)  an  effective  (r_6)-potential  and  on  (b) 
the  polarization  potential.  This  requirement  is  equivalent  to  an  equality  of  mean 
momentum-transfer  cross-sections  (at  STP  conditions)  based  on  the  same  two 
potentials,  since  the  mobility  is  (inversely)  proportional  to  the  mean  momentum- 
transfer  cross-section.  The  effective  (1 — 6)-potential  is  also  applied  to 
[Ba++(N2,C>2)]-collisions. 

In  Appendix  D  we  introduce  the  ambipolar  diffusion  coefficient, 

Da* Dj{\+Tel Tt)  (where  D<  is  the  ordinary  diffusion  coefficient,  and  Te  and  T,  are 
the  electron  and  ion  temperatures),  and  present  limited  comparisons  of  (a) 
theoretical  results  for  the  coefficient  of  ambipolar  diffusion  of  Ba+  along  the 
direction  of  the  geomagnetic  field,  l\  (taken  to  be  equal  to  Oa),  with  (b) 
experimental  results  for  some  barium  release  experiments.  The  motivation  is  to 
see  whether  such  experiments  can  provide  a  basis  for  preferring  either  the 
polarization  potential  or  the  effective  (r_6)-potential  in  computing  the  ambipolar 
diffusion  coefficient.  Such  comparisons  have  many  uncertainties,  including  those 
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for  selection  of  a  model  (neutral -particle)  atmospnere  and  the  assignment  of  the 
ratio  of  the  electron-  and  ion-temperatures.  An  overall  comparison  of  0„  E(( 

with  A.  exp  (derived  from  12  Ba-reieases)  is  about  the  same  as  for  ^  so 

that  a  strong  preference  cannot  be  established  on  this  basis,  unless  one  should 
have  reasons  (unknown  to  the  writer)  to  select  or  exclude  certain  experimental 
data. 


In  Appendix  E  we  present  the  relative  speed  for  a  monoenergetic  beam 
traversing  a  Maxwellian  gas,  based  on  two  approximate  formulas  and  the  exact 
formula;  one  of  the  approximate  expressions  is  used  in  Section  3.  Appendix  F 
presents  the  differential  and  integral  Maxwellian  energy  distributions 
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SECTION  2 


MOMENTUM-TRANSFER  CROSS-SECTIONS 

A  very  useful  and  convenient  survey  of  momentum-transfer  cross-sections 
is  given  by  Stoeckly,  Stagat,  and  Kilb  [Ref.  65],  but  values  directly  applicable  to  U- 
release  species  and  conditions  are  not  given.  Hence,  our  approach  will  be  (a)  to 
present  some  general  formulas  for  atom-atom  and  ion-atom  collisions,  (b)  to  apply 
and  extend  those  formulas  to  [(0,0+)+(0,N2,02)]-collisions,  and  (c)  to  consider 
procedures  for  uranium-species  collisions. 

2. 1  ATOM-ATOM  AND  ION-ATOM  COLLISIONS;  GENERAL. 


2.1.1  Atom-Atom  Collisions. 


2. 1.1.1  Lennard-Jones  ( 1 2.6)-Potential.  For  the  Lennard-Jones  ( 12,6)- 
potential,  which  is  written  in  alternative  forms, 


<#>(/-)  =  4  e0 


lr  \b 


r 

-2\f 


( 1 a,b) 


D  =  rm/  '2*  =  0.8909  rm 


( l  c) 


Hirschfelder,  Bird,  and  Spatz  [Ref.  31]  tabulate  the  reduced  momentum-transfer 
cross-section  5° ’Or)  (in  units  of  xD2)  for  the  reduced-energy  range  0.1  srs  100 
(from  which  selected  values  will  be  presented  later).  Here,  c0  is  the  depth  of  the 
potential  at  its  minimum,  rm  is  the  radius  at  which  the  minimum  energy  (-$) 
occurs,  and  D  is  the  radius  at  which  the  potential  is  zero.  The  cross  section  nD2 
is  interpreted  as  that  for  rigid  spherical  particles  of  collision  diameter  D.  Hence, 
for  such  spheres,  5<l,Gr)  would  be  unity.  The  reduced  energy  *•  is  defined  as 


(2a) 

(  2b  I ) 
(2b2) 

where  ft  is  the  reduced  mass  for  the  colliding  particles, 


**  -  ^cfi  I  £<> 

^CM  =  2  AiT2 

Me V)  =  5.1 88x1  O'3 /t(A)  ^ 


A  = 


(3a) 
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/i(A)  = 


44 
4  +4  ’ 


/%  =  1.66xi  Cr^g, 


(3b,c) 


and  g  is  the  relative  velocity.  To  relate  g  to  r  we  write 

-|l/2 

\  /  r  r  \  .  \  r  r 

9 - 


2rf  at 

1/2 

=  1.388x1 06 

L  A  . 

/l(  A) 

(4a, D) 


The  momentum-transfer  cross-section  <r(^),  based  on  the  Lennard- Jones  (  1 2,6)- 
potential,  is 

(5) 


2. 1.1. 2 


Asymptotic  Extension  to  l  ow  Energies.  To  extend  these  cross 
sections  to  lower  velocities,  we  use  asymptotic  values  derived  Dy  Moncmck  and 
Mason  [Ref.  59]  who  state  that  the  reduced  cross-section  (in  units  of  nD1), 


5^=2 


-(-l)41"1 


2  1 +A 
takes  the  asymptotic  form 

1  l+(- 1)4 

2  1+^ 


f  0-cosA  /)ddd , 

Jo 


i^A) 


-i-i  r 


L2 


nd 
-  M  92 


12  m 


A  fin) 


at  low  energy  for  which  the  potential  becomes,  in  general, 
dKr)  =  t\d\\ 


D_ 

r 


(6) 


(?) 


(8) 


in  Eq.  (7),  AA\n)  is  a  pure  number  whose  value  depends  only  on  the  power  n  of  the 
potential  and  on  whether  the  potential  is  repulsive  (subscript ♦  )  or  attractive 
(subscript  -).  Reference  59  gives  the  values  of  AfXn)  in  its  Table  I,  taken  from 
Kihara,  Taylor,  and  Hirschfelder  [Ref  38]. 

For  the  lennard-Jones  ( 12,6)-potential,  at  low  energy  we  have 
0(r)=>-4co 
n  =  6, 


(9) 


Thus, 


d  =  4c„ 


and  we  need  A']( 6),  given  in  Table  I  of  Ref.  59  as 
A"{ 6)  =  0  4342 


( 10a, b) 


(II) 
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Thus,  from  Eqs.  (7),  (3),  ( i 0),  and  ( 1) ),  we  have 


5™  *3.1 1 7x1 04 


t)/3 


fiKA)  g 2 


( 1 2a) 


By  combining  Eqs  (5),  ( 1  c),  and  ( 1 2a)  (cf  Appendix  A  3  1  2),  we  have 


(&^en)u\r.  [rle^ev) 

/ig1  )  fiiA)g2  ) 


6  >-pot 


=c 


\l/3 


=  r(  r%€0Ke V) 


\l/3 


T  =  9.204x1  O'4, 
C"  =  36.08 


l  fiA)  7km/ sec  1 

C  -  7.773x1 04 


2.1.2  Ion-Atom  Collisions. 


( 1 2D. o 

( 1 2ci) 

( 1 2e.f ) 
dig; 


2. 1.2. 1  Polarization  Potential.  At  sufficiently  low  velocities  where  the 
polarization  force  between  an  ion  and  neutral  dominates,  the  velocity-dependent 
momentum-transfer  cross-section  for  sinah'-^harqed  ions  [Ref  7,  Pt  A, 
p  218]  is 


<J a  =  2.2\  Jr  [a  t>z  /  /i  g2]”7  =  259^0'*  [ar/  /i(A)\yl  g~' ,  ( 1 3a. 0) 

where  a\s  the  polarizability  of  the  neutral  particle  (with  a  in  units  of  iO-24 
cm3),  e  is  the  electronic  charge,  and  /i  and  g  are  as  defined  in  Section  2  1  1 

This  cross  section  is  only  a  little  larger  than  that  for  "orbiting  collisions" 
[Ref  55,  p.  72], 


q0*ntf0*2n{ae  /  fig1]'2  ^nXlae1  /  fwY2  ( 1 3c  1  • 

4,(1  O'4  cm)  =  8635  [ a 7  /z(4)f4  ( '  oc2  > 

For  an  impact  parameter  equal  to  b0,  the  ion  spirals  into  a  circular  orbit  of  radius 

r9  =  t>ji 2  (13d) 

From  Ref  55,  p.  72:  "Orbits  for  which  b  <  &c  pass  through  the  origin  if  no  repulsive 
core  is  present,  whereas  those  orbits  for  which  t?  i  t>0  come  no  closer  than 

2  " 


9 


2. 1.2.2  (l2.6.4)-Potential  The  (l2,6,4)-potential  [cf.  Ref  65]  is 


<P(r)  =  ^ 


(1 +y) 


<r  f 

*  m 

\r) 


r„ 


\6 


-M'-fl  -3d-/) 


r. 


+rb  +  r4], 


( !4a! 


t  i 


(14a2) 


where  the  dimensionless  parameter  y  is  a  measure  of  the  relative  contributions 
of  the  4  and  To,  terms  For  y  =  1  we  recover  the  Lennard-Jones  (12,6)-potentiai 
given  by  Eq.  ( 1 )  and  for  y  =  0  we  have  the  ( 1 2,4)-potential.  To  determine  y ,  as 
was  done  in  Ref.  65,  we  require  the  Ta,  term  to  match  the  polarization  potential, 


whereupon 


polarization ' 


<?at 

'2r 


A  > 


Y  =  1- 


tf2 


a 


r 4 

^  co  m 


=  1-4.806 


a ■ 


eSeV )  (c)4’ 


with  in  units  of  10~8  cm. 


(14b) 


(14c) 


we  will  be  interested  in  the  momentum-transfer  cross-sections 
corresponding  to  the  full  ( 12,6,4)-potential.  Such  cross  sections  (in  units  of 
nrj„  )  have  been  computed  for  y-  values  of  0,  0.25,  and  0  50  by  Mason  and  Scnamp 
[Ref  51,  Table  2;  Ref  52;  Ref  56,  Table  1-3],  Similar  cross  sections  (in  units  of 
x D1)  have  been  computed  for  the  ( 1 2,6)-potential  -  which  implies  y=  1  for  the 
( 1 2,6,4)-potential  -  by  Hirschfelder  et  al.  [Ref.  3 1  ].  The  calculations  [Ref  3 1  ]  for 
y=  1  were  for  a  different  set  of  /r-values  than  were  those  by  Ref.  51,  hence,  the 
first  step  in  using  these  results  was  to  interpolate  Ref.  31  results  to  obtain 
results  for  the  same  set  of  rs  as  in  Ref,  51.  Next,  for  later  use,  we  interpolated 
between  the  results  for  y  =  0.5  and  y=  1  to  obtain  results  for  y  =  0.75.  Finally, 
we  interpolated  again  to  obtain  results  at  y-  values  of  interest  for  (0++0)-  and 
(U++0)-collisions,  i.e.,  at  respective  y-  values  of  0.6869  and  0.960.  The  results 
from  Ref.  51  and  Ref.  31  (the  latter  multiplied  by  (2,/6)-2  )  and  the  interpolated 
results  are  given  in  Table  1. 
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Table  l .  Reduced  momentum-transfer  cross-sections  for  ( I  2,6,4)-potential 
(in  units  of  jr/^j). 


tan*'  k 

K *  v=0 

0.25 

0.50 

0.6896° 

0  75 

0.960e 

v=  1 

0.2a 

0.0 

0.1 

0.2 

0.2027  la  5.50  la 

2.7072 

4.9 18a 

4.3068 

3.95c 

3.84C 

3.57C 

4.322b 
3  553b 
3  52c 

0.4a 

04 

0.422799  3.7 1 9a 

3.432a 

3.l37a 

2.93C 

2.87C 

2.64° 

2  673° 

0  ^oc 

_  •  V*  S/ 

0  6a 

0.6 

0  6841 4a  2.892a 

2.730a 

2  506a 

2.39C 

2.35C 

2  22c 

2.309° 
2  19° 

0  7a 

0.8 

0  84229a  2.394a 

2.333a 

2.258a 

2.I5C 

2.12c 

2.01C 

2.002° 

1  99c 

0.8a 

09 

0.95 

1.0 

1 .0296a  1 ,974a 

1 .950a 

t  926a 

1.91C 

1 .9  lc 

1.89° 

1  988b 

1  958b 

1  906° 

1 .89° 

a  Ref.  56,  Table  1-3.  d  ForO++0,  see  Fig.  5. 

b  Ref.  31,  Table  VIII,  multiplied  by  (21/6)-2  =  0.79370.  e  For  IT +0,  see  Fig  9 
c  Interpolated. 

2.2  [(0,0+)+(0,f^,02)] -COLLISIONS. 

2.2.1  (0+0)-Collisions 

Values  of  J^Gr)  for  the  range  0.  U  tr  i  2  [HB-48]  are  reproduced  in  Table  2 
To  relate  these  data  with  results  in  Ref.  65,  we  apply  the  data  to  (0+0)-collisions. 
The  potential  energy  parameters  [Ref.  65,  p.  40]  are 

^  =  5.22eV,  /^- 122x10-*  cm.  (15a,b) 

so  that  with  use  of  Eq.  ( 1  c)  we  have 

D  =  1.09x1  O'®  cm.  (15c) 

From  Eq.  (4)  the  relative  velocity  g  is  related  to  tc  by 
1.1 2 1x10  V* 


(16) 


The  momentum-transfer  cross-section  ai.g),  Eq.  (5),  is  given  in  units  of  10'16 
cm2  jn  Table  2  and  plotted  in  Fig.  I.  These  values  for  (0+0)-collisions  agree  wen 
with  those  plotted  in  Fig.  4-3  of  Ref.  65. 


Table  2.  Velocity-dependent  momentum-transfer  cross-sections 
from  Ref.  31  for  Lennard-Jones  (!2,6)-potential  and 
application  to  (O+O)-  and  (U+0)-collisions. 


JC 

tf'KA')? 

_ 0 

+  0 

U 

+  0 

(Tb 

cm/sec 

0.1 

5.445 

3.54  E+05 

20.2 

3.16  E+05 

45  9 

0.2 

4.476 

5.01  E+05 

16.6 

4.47  E+05 

37.8 

0  4 

3.368 

7.09  E+05 

12.5 

6.32  E+05 

284 

0.6 

2.909 

8.68  E+05 

10.8 

7.74  E+05 

24.6 

0.8 

2.523 

1  00  E+06 

9.36 

8  93  E+05 

21  3 

0.9 

2.505 

1.06  E+06 

9.30 

9.47  E+05 

21.1 

0.95 

2.467 

1.09  E+06 

9.16 

9.73  E+05 

20.8 

1.0 

2.401 

1.12  E+06 

8.91 

999  E+05 

20  3 

1.2 

2. 1 007 

1.23  E+06 

7.80 

1.09  E+06 

17.7 

1.4 

1.8714 

1.33  E+06 

6.95 

1.18  E+06 

158 

1.6 

1.6974 

1.42  E+06 

6.30 

1.26  E+06 

143 

1.8 

1.5656 

1.50  E+06 

5.81 

1.34  E+06 

13.2 

2.0 

1.4801 

1.59  E+06 

5.49 

1.41  E+06 

12.5 

Fig.  1 

Fig  6 

a  Reduced  cross-sections  [J01^)],  in  units  of  n  D2,  from  Ref  31,  Table  vm 
Potential  parameters  are  given  by  Eq.  (15)  for  (0+0)  and  by  Eq  (22;  for 
(U+0). 

b  Cross  sections  are  in  units  of  10" 16  cm2. 

To  extend  these  cross  sections  to  lower  velocities,  we  use  Eqs.  (5)  and  ( 1 2) 
to  give  us 


<7(0+0)|(/._«)_pot=  1.00  x  1 0'1 1  g~2'1  cm2  (17a) 

=  4.64x!0-'5c^;ic  cm2  ,  (17b) 

For  g  =  3.54x1 0s cm/sec,  Eq.  (17)  gives  20x1  O'16  cm2,  in  excellent  agreement  with 
the  lowest-energy  value  in  Table  2.  We  will  adopt  Eq.  (17)  as  the  appropriate 
relative-velocity  momentum-transfer  cross-section  for  (0+0)-collisions  for  the 
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relative  velocity  range  less  than  about  5x!05  cm/sec,  shown  as  the  dashed  line 
(s3xltf  cm/ sec)  in  Fig.  I 


Figure  1 .  Velocity-dependent  momentum-transfer  cross-sections  for  (OO)- 
collisions,  based  on  lennard-Jones  ( l2,6)-potential  and  asymptotic 
(i — 6)-potential 


2.2.2  [0+(N2,02)]-Collisions. 

To  eventually  obtain  cross  sections  for  [U (N^ ,02 )] - co  1 1  i s i ons ,  we  will 
appeal  to  collision-integral  results  for  [0+(0,N2,02)]-collisions,  from  which  we 
infer  an  approximate  scaling  factor  to  relate  [0+(N2,02)]-collisions  to  (0+0)- 
col  1  isions  That  same  scaling  factor  will  then  be  applied  to  (U-^O)-col lisions 

We  would  prefer  to  see  momentum-transfer  cross-sections  for 
[0 +(0,N2,02)]-col  1  isions,  but  what  we  find  to  be  available  (without  an  exhaustive 
search)  are  collision  integrals ,  i.e.,  momentum-transfer  cross-sections  averaged 
over  the  relative  velocity  distribution  in  a  Maxwellian  gas  at  temperature  T  Such 
a  cross  section  (or  collision  integral  for  diffusion)  is  defined  by  Eq.  (A8c)  in 
Appendix  A.  Yun  and  Mason  [Ref.  69]  have  computed  and  tabulated  such  integrals 
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(called  (/2t,,l))  in  their  notation  and  expressed  in  units  of  ICr16  cm2  but  without  a 

necessary  factor  of  *to  make  them  a  conventional  cross  section)  for  the  various 
species  in  dissociating  air  for  the  temperature  range  1 000  i  AK)  i  15000  (it 
appears,  notwithstanding  an  otherwise  description,  that  Ref.  69  accounts  for  a 
large  number  of  excited  states  in  computing  averages  for  atom-atom  collisions, 
whereas  we  have  limited  our  considerations  to  ground-state  interactions  )  We 
will  regard  these  collision  integrals  as  pseudo  velocity-dependent  momentum- 
transfer  cross-sections  by  expressing  them  as  a  function  of  the  average  relative 
velocity  between  the  two  members  of  each  of  the  three  pairs  of  particles  of 
interest,  i.e.,  O+O,  O+Nj,  and  0+02.  In  such  a  gas  the  average  relative  velocity  [Eg. 
(A8b)]  is 


1.45x1 0-1 


where  the  reduced  mass  number  is  8,  10.18,  and  10.67  for  the  three  stated 
pairs.  Values  of  g,k  are  given  in  Cols.  2,  3,  and  4  of  Table  3.  Collision  integrals 
presented  as  a  function  of  temperature  in  Ref.  69  (Tables  III,  VII,  and  VI)  are  given 
here  (but  multiplied  by  n)  in  Cols.  5,  6,  and  7  of  Table  3  as  a  function  of  the 
temperature  and  in  Fig.  2  as  a  function  of  the  average  relative  velocity. 


Figure  2.  Pseudo  momentum-transfer  cross-sections  versus  mean  relative 
velocity  <g>m[Qk T / n fiY1,  inferred  from  temperature-dependent 
collision  integrals  [Ref.  69]  for  [0+(0,N2,02)]-collisions  at  gas 
temperatures  in  the  range  1 000 ^  7"(K)  ^  1 5000 
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Table  3  Pseudo  momentum-transfer  cross-sections  versus  mean  relative 

velocity,  inferred  from  collision  integrals  [Ref  69]  for  [0(0,N2,02>!- 
collisions  for  gas  temperatures  in  the  range  1000  s  T0<)  s  15000 


r 

<?,*,  km/  sec 

n 

io-'6 

cm2 

1 000  K 

0+0 

0+N, 

0+0, 

0+0 

0+N, 

0+0, 

1.0 

1.62 

1.44 

1.40 

19.4 

27.9 

23.6 

1.5 

1.99 

1  76 

1.72 

177 

25.1 

21  6 

2.0 

2.29 

2.03 

1.99 

16.6 

23.2 

20.3 

2.5 

2.56 

2.27 

2.22 

15.6 

21  7 

19  3 

3.0 

2.8! 

2.49 

2.43 

14.9 

20.6 

18  5 

3.5 

3.03 

2.69 

2.63 

14.3 

19.7 

17  8 

40 

3.24 

2.87 

2.81 

138 

189 

173 

45 

3.44 

3.05 

2.98 

13.3 

18.2 

16.8 

5.0 

3.62 

3.21 

3  14 

1 2.9 

176 

164 

5.5 

3.80 

3.37 

3.29 

12.6 

17.1 

16.0 

6.0 

3.97 

3.52 

3.44 

12.3 

16.6 

15.6 

6.5 

4.13 

3.66 

3.58 

12.0 

16.2 

153 

7.0 

4.29 

3.80 

3.71 

1 1.7 

15.8 

15  0 

75 

444 

394 

3  84 

1 1  5 

154 

146 

8.0 

4.59 

4.06 

3.97 

11.2 

15.1 

145 

8.5 

4.73 

4.19 

4.09 

1 1.0 

14.8 

143 

9.0 

486 

4.31 

4.21 

10.8 

14.5 

14.1 

9.5 

5.00 

4.43 

4.33 

10,6 

14.2 

13.9 

10.0 

5.13 

4.54 

444 

10.5 

13.9 

13  7 

1  1.0 

5.38 

4.77 

466 

10.2 

13.5 

13.3 

12.0 

5.62 

4.98 

486 

9.89 

13.0 

13.0 

130 

584 

5.18 

5.06 

964 

12.7 

127 

140 

6.07 

5.38 

5.25 

9.41 

12.3 

12.4 

15.0 

6.28 

5.57 

544 

921 

12.0 

122 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

Fig,  2 

Fig,  2 

Fig  2 

15 


It  is  of  interest  to  relate  the  (OO)-values  in  Fig.  2  to  the  asymptotic  values 
given  by  Eq.  (17).  Firstly,  in  Fig.  2  the  slope  is  closely  -2/3  at  the  higher 
velocities,  as  in  Eq.  ( 1 7),  but  it  decreases  at  the  lower  velocities.  Secondly, 
substitution  of  Eq.  ( 1 7)  into  Eq.  (A8c)  yields 

Q& O+O)  =  2.74x1 0 r  V /I 000 T™ .  (19) 

At  1000  K,  Eq.  (19)  gives  a  value  larger  than  the  1000-K  value  for  (OO)-collisions 
in  Table  3  by  the  ratio  2.74/1.94  =  1  41  (a  ratio  which  would  be  about  1.21  if  the 
slope  in  Fig.  2  did  not  decrease).  The  lack  of  complete  consistency  between  Eq 
( 1 7)  and  Fig.  2  may  be  due  partly  to  the  pseudo  nature  of  the  cross  sections  in 
Table  3  and  Fig.  2  and  partly  to  the  aforementioned  inclusion  in  Ref.  69  of  excited 
states  with  their  concomitant  variety  of  potential  curves. 

We  emphasize  that  we  are  less  concerned  with  the  absolute  values  in  Fig.  2 
than  the  values  for  (O+N2)-  and  (0+02)-collisions  relative  to  those  for  (OO)- 
collisions.  We  see  that,  viewed  very  simply,  5/4  is  a  factor  by  which  we  can 
convert  from  (OO)-cross-sections  to  approximate  mean  values  for  (ON2)-  and 
(OC>2)-collisions. 

2.2.3  (0++0)-Collisions. 

Since  we  will  later  describe  some  reservations  about  applying  the 
polarization  potential  for  uo;-collisions  with  atmospheric  species,  we  will 
consider,  first,  (0++0)-collisions  for  which  we  feel  there  is  a  relatively  firm 
basis. 

We  now  apply  Eq.  (14)  for  the  (!2,6,4)-potential.  For  0,  <r=0.79  (cf.  Table 
A- 1 );  for  0*,  rm  *1.164x10'®  cm  and  c0  *  6.663  eV  [Ref.  34,  p.  504],  Thus,  from  Eq. 
( 14c)  we  find 

/(0++0)  =  0.6896.  (20) 

In  Fig.  3a  we  plot  -<Pir)  from  Eq.  ( 1 4a)  and  in  Fig.  3b,  the  ratio  ( R)  of  the  A 
term  to  the  sum  of  the  A  and  A  terms: 


In  Fig.  4  we  plot  b0  from  Eq.  ( 1 3c),  r0  from  Eq.  ( 1 3d),  and  the  center-of-mass  energy 
fen  from  Eq.  (2b).  Figures  3  and  4  enable  one  to  consider  the  relative  magnitudes  of 
the  A  and  A  terms,  for  various  velocities,  as  given  in  Table  4.  We  see  that  for 
velocities  below  3  km/sec,  the  A  term  is  at  least  half  of  the  total  of  the  two 

attractive  terms  even  at  the  orbit-radius  ^  =  4/1/2.  For  such  velocities,  use  of  the 
A  term  alone  would  probably  be  satisfactory  in  computing  the  momentum- 
transfer  cross-section,  shown  in  Fig.  5  as  the  solid  line  computed  from  Eq.  (13d), 
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Figure  3b.  Radial  dependence  of  the  ratio  of  the  (r4)-term  to  the  sum  of  the 

(r4)-  and  (r_6)-terms  in  the  (l2,6,4)-potential  for  (0++0)-collisions 
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Figure  4.  Relative  velocity  dependence  of  the  "orbiting  impact  parameter"  4, 
"orbiting  radius”  r0,  and  center-of-mass  energy  for  (0++0)- 
collisions. 


Table  4.  Comparison  of  the  A  and  %  terms  in  the  (12,6,4)-potential 
for  (0+*0)-collisions. 


g 

io-® 

cm 

.  4, 

Ea.(21 ) 

km/sec 

4 

eV 

m) 

/?( 4/V2) 

0.3 

883 

6.25 

3.74  E-03 

*0.96 

0.91 

1 

4.84 

3.42 

4. 1 5  E-02 

0.85 

0.74 

3 

2.79 

1.98 

3.74  E-01 

0.66 

0.50 

5 

2.16 

1.53 

1.04  E+00 

0.54 

0.37 

<7(0* 

+  O]|polpot.  = 

8.1 4x1  O’10 

g~'  =  8.14x10* 

15 

(1301,02 

The  cross  sections  corresponding  to  the  reduced  cross-sections  in  Table  1  (for  the 
column  corresponding  to  Eq.  (20))  are  plotted  as  the  circled  points  in  Fig.  5  for  the 
relative  velocity  range  5.7  s  ^(km/sec)  s  12.86.  The  dashed  line  in  Fig.  5 
extrapolates  these  values  to  about  3  km/sec  where  the  line  intercepts  the  cross 
section  [Eq.  ( 1 3d)]  for  the  pure  polarization  potential. 
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Figure  5.  Velocity-dependent  momentum-transfer  cross-sections  for 

(0++0)-  ->’’,isions,  based  on  ( 1 2,6>4)-potential,  its  extrapolation, 
and  Do'anzation  potential. 


2.3  [(U,L)+,U0+)+(0,N2,02)]-C0LLISI0NS. 

2.3. 1  Introduction. 

In  midcourse  of  our  U-release  studies,  we  paused  to  study  Ba+  collisions 
with  atmospheric  species,  hoping  that  useful  insights  would  be  obtained.  That 
study,  documented  in  Appendices  A,  B,  C,  and  D,  revealed  the  following  important 
facts: 

a.  The  experimental  mobility  of  Ba+  in  N2  (at  room  conditions)  agrees  well 
with  the  theoretical  value  based  on  the  polarization  (r"4)-potentiai 

b.  The  (^-contribution  dominates  the  (r_4)-contribution  in  the  ( 1 2,6,4)- 
potential  for  (Ba++0)-col1isions  (analogous  to  the  situation  for  [(U,U+)+Ol- 
collisions,  as  will  be  shown  later). 

c.  An  effective  value  of  the  parameter  (/‘/Jr„)’/5  -  which  appears  in  the 
asymptotic,  low-energy  momentum-transfer  cross-section  based  on  the 
(r6)-potential  -  is  postulated  to  be  determined  by  requiring  the  equality 
of  mobilities  (at  STP  conditions)  based  on  an  effective  (r^)-potential  and 
on  the  (r_4)-potential.  (At  non-STP  conditions,  the  two  mobilities  will 
differ.) 
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d  Values  of  the  ambipolar  diffusion  coefficient  [/?„  Ef(  (r^)  pot  ]  computed  from 

an  effective  (r~6)-potential  are  -  for  ionospheric  ion  temperatures  ranging 
from  1 135  K  to  750  K  -  about  79%  to  85%  of  the  ambipolar  diffusion 
coefficient  [/}  ^  ^  ]  based  on  the  polarization  potential  An  overall 
comparison  of  Eff  ^.6,.^  with  (sometimes  ambiguous  or  uncertain; 

experimental  values  [Z}(  #xp  ]  derived  from  12  Ba-releases  is  about  the  same 
as  for  L \  ^  ^ ,  so  that  a  strong  preference  cannot  be  established  on  this 
basis 

it  would  be  desirable  to  have  an  experimental  mobility  for  U+  in  Nj  for  ar.v 
gas).  Not  having  that  information,  we  will  make  the  qualitative  extrapolation 
from  (Ba*+N2)-collisions  to  [UMO,N2.02)]-collisions  and  assume  that  the 
mobilities  of  U+  in  atmospheric  species  at  STP  conditions  can  be  predicted  by  use 
of  the  ti — 4)-potent iai  We  will  also  assume  that  effective  values  of  the  parameter 
(/£rc,)'/3  can  be  determined  by  requiring  equality  of  the  mobilities  (at  STP 
conditions)  based  on  an  effective  (r~6)-potential  and  on  the  (r4)-potential.  These 
assumptions  allow  us  to  use  all  the  general  formulas  we  developed  in  the 
Appendixes  for  Ba+  collisions,  and  we  will  refer  to  them  when  convenient.  For 
completeness,  we  will  also  include  the  work  performed  prior  to  our  consideration 
of  Ba*  collisions.  Table  5  collates  numerous  parameters  needed  in  the  subsequent 
discussion 

Table  5.  Parameters  for  UCC  collisions  with  neutral  atmospheric 
species 


K 

a; 

k 

/=u,ir 

uo* 

uo; 

/  =  U,U* 

UO* 

0 

N2 

02 

K 

0  79 
1.76 

1  59 

1499 

25.05 

2821 

[fiiA 

15  05 
25.22 
28  42 

„>rw 

15.10 

25.37 

2861 

0  230 
0.265 
0.237 

[ak  /t(A/k)f2 

0  229 
0.264 

0  236 

/  =U,l / 

uo* 

/  =  U,U+ 

UO* 

1 1  n* 

0 

0.360 

0.360 

3  44 

3  45 

3  45 

N2 

0.453 

0  452 

6.64 

6.66 

6  68 

02 

0414 

0413 

6  70 

6  72 

6  74 
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2.3.2  [U+(0,N2,02)]_Collisions 


2-3.2. 1  (U+Q)-Collisions  To  obtain  the  momentum-transfer  cross-section  for 
(U+0)-co!lisions  based  on  the  (r-6)-potential,  we  adopt  values  of  the  potential 
parameters  available  at  the  beginning  of  this  study,  i.e.,  those  qiven  by  Michels 
[Ref.  57,  p.  185], 


r„*7.76eV,  rm*rf  =  1.84x10"®  cm  (22a, b) 

D- 0.8909/^  =  1.64x1 0‘d  cm.  (22c) 

To  relate  the  relative  velocity  to  tr ,  we  use  Eq.  (3b)  to  get 

/fC4)-l4.99,  (23a) 

which,  with  Eqs.  (22)  and  (4b),  gives 

g  =  9.987x1  OV72  (23b) 


The  momentum-transfer  cross-section  from  Eq.  (5),  based  on  the  Lennard-Jones 
( 12,6)-potential, 

a{g)  -  8.44x1 0~,65<,)(*‘) ,  (24) 

is  given  in  units  of  10*'6  cm2  in  Table  2  and  plotted  in  Fig.  6. 

To  obtain  cross  sections  at  velocities  below  about  3x1 05  cm/sec,  we  use  Eq 
( 1 2a),  with  Eqs.  (22a),  (22c),  and  (23),  to  get 

^(u+  0 1 -.Wpot  =  2. 11  x  1 0- ' V*273  =  9. 79 x  1 0* 15  fig"  cm2  (25a, b) 

For  ^  =  3.1 6x1 05  cm/sec,  Eq.  (25)  gives  45.5xlCr16  cm2,  in  excellent  agreement  with 
the  lowest-energy  value  in  Table  2.  We  will  adopt  Eq.  (25)  as  the  appropriate 
relative-velocity  momentum-transfer  cross-section  for  (U+0)-col  1  isions,  based  on 
the  (r-6)-potential,  for  the  relative  velocity  range  less  than  about  5xio5  cm/sec. 
as  shown  as  the  dashed  line  u  3x1 05  cm/sec)  in  Fig.  6. 

To  obtain  the  cross  section  for  the  effective  (r^)-potential,  we  use  Eq. 
(A20)  and  the  parameters  in  Table  5  to  get 

1.38x1 0-"  =  6.41x1 0-15  cm2 .  (26a, b) 

Equation  (26)  is  also  plotted  in  Fig.  6. 
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Figure  6.  Velocity-dependent  momentum-transfer  cross-sections  for  (U+0)- 
collisions,  based  on  Lennard-Jones  ( 12,6)-potential,  asymptotic 
(r_6)-potential,  and  effective  (r*&)-potential. 


2.3.2  2  [U-t-CNg.O^)]— Collisions.  For  an  interaction  we  will  label  *(r_&)- 
potential  (even  though  it  is  straining  the  concept),  we  appeal  to  results  for 
[0(0,N2,C>2)]-col!isions  and  apply  the  approximate  scaling  factor  inferred  in 
Section  2.2  2  to  relate  [(MN2, (^-collisions  to  (0+0)-collisions,  i.e.,  5/4  By 
applying  this  factor  to  Eq.  (25)  we  have 


<7(U*(NJ,0,)1.,.1Veol  =  |  <r(U+0)|(r.,^  <27a  I ) 

=  264x10‘,V2/3  cm2  (27a2) 

=  1.24x1 cm2.  (27a3) 

For  the  effective  (r6)-potential,  we  use  Eq.  (A20)  and  parameters  in  Table  5  to 
get 

<r(U+  N2)| Eff (r ^  =  1. 7 4 x  1  O’ " g'™  cm2  (27b  1) 

=  8.08x1  O'’5  cm2  (27b2) 

02)|Erf  (r.)_pot  =  l.59x  1 0*'  y 2/3  cm2  (27c  1 ) 

=  7.38x1 0-’5^  cm2  (27c2) 
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2.3.3  [U++(0,N2,02)]-ColllSlons. 


2.3.3. 1  (U* +0)-Co))isions.  For  (U++0)-co!!isions  characterized  by  the 
( 1 2,6,4)-potential  given  in  Eg.  ( 1 4a),  we  use  the  UO+  parameters  [Ref  57,  p  1 961 

4*8.15eV,  rm « 1.843x10^ cm  (28a, D) 

and  find  from  Eq.  (14c)  that 

y(LT+0)  =  0.960.  (28c) 

We  interpret  the  closeness  of  Y  to  1  as  a  strong  hint  that  the  polarization 
potential  by  itself  is  unsatisfactory  for  (U++0)-collisions  in  the  velocity  range  of 
current  interest. 

In  Fig.  7a  we  plot  -<Pir)  from  Eq.  ( 1 4a),  evaluated  for  the  parameters  given 
in  Eq.  (28),  and  in  Fig.  7b,  the  ratio  (/?)  given  by  Eq.  (21 ).  In  Fig  8  we  plot  4  from 
Eq.  ( 1 3c),  r0  from  Eq.  ( 1 3a),  and  the  center-of-mass  energy  4cn  from  Eq.  (2b).  In 
Table  6  we  collate  the  comparison  of  the  7a  and  7s  terms  for  (u++0)-coilisions. 
The  smallness  of  the  Ta  term  relative  to  the  total  of  the  two  attractive  terms, 
even  at  the  lowest  velocity  of  interest  here  (0.3  km/sec),  seems  to  be  another 
strong  hint  that  the  polarization  potential  by  itself  is  unsatisfactory  for  (U++0)- 
collisions.  In  Fig.  9  the  momentum-transfer  cross-section  for  the  polarization 
potential  is  shown  as  the  short-dash  line,  computed  from  Eq.  ( 1 3b)  and 

<7(U++0)|PoipOl=  5.96x1 0-10^'1  cm2  (29a) 

-  5.96x1 0-'Visec  cm2  (29b) 

Table  6.  Comparison  of  the  £*and  Ts  terms  in  the  ( 1 2,6, 4)-potentiai 
for  (U++0)-collisions. 


9 

km/sec 

1  O’8  cm 

4, 

Ed.(21 ) 

4 

e>jV2 

ev 

R(  VV2) 

0.3 

7.55 

5.34 

7.00  E-03 

0.35 

0  21 

l 

4.14 

2.93 

7.78  E-02 

0.14 

0.073 

5 

1.85 

1.31 

1 .94  E+00 

0.031 

0016 

Again,  we  are  interested  in  using  the  f ul  1  ( 1 2,6,4)-potentia!  to  obtain  the 
momentum-transfer  cross-section,  (Owing  to  the  closeness  of  ^(U++0)  to  unity 
per  Eq.  (28),  we  expect  little  difference  from  using  just  the  ( 12,6)-potential )  The 
interpolated  values  of  the  reduced  cross-section  are  given  in  Table  1  for  /  =  0.960 
and  the  actual  cross-sections  are  given  in  Table  7  and  plotted  in  Fig.  9  as  the 
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Figure  8.  Relative  velocity  dependence  of  the  "orbiting  impact  parameter"  4, 
"orbiting  radius"  r0,  and  center-of-mass  energy  for  (U++0)- 
coiiisions. 


10"13 


* 

io-14 


10-'5 

Figure  9.  Velocity-dependent  momentum-transfer  cross-sections  for  (U+*0)- 
collisions,  based  ond^^-potential,  Lennard-Jones  (12,6)- 
potential,  asymptotic  (r6)-potential,  effective  (r6)-potential,  and 
polarization  potential. 
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solid-square  points.  The  velocities  in  Table  7  are  computed  from  Eq.  (4b),  with  r, 
given  by  Eq.  (28a)  and  /zC4)  =  1499  from  Table  5,  giving 


^(km/sec)  =  10.24-/r .  (30) 

Table  7.  Momentum-transfer  cross-sections  based  on  the  (12,6,4)- 
potential  with  /  =  0.960  for  (U++0)-col1isions. 


tr 

.<7,  km /sec 

o,\  O'15  cm2 

0.20271 

4.61 

3.57 

3.80 

0.42279 

6.66 

2.64 

2.82 

0.68414 

8.47 

2.22 

2.37 

0.84229 

9.40 

2.01 

2.14 

1  0296 

10.39 

1.89 

2.02 

Table  1 

Eq.  (30) 

Table  1 

Col.3&Eq.  (28b) 
Fig.  9 

For  completeness  we  also  use  the  ( !2,6)-potential  to  obtain  the  momentum- 
transfer  cross-section  for  (U++0)-collisions.  The  reduced  energies  and  cross 
sections  (in  units  of  nD2)  from  Ref.  31,  with  the  velocity  from  Eq.  (30),  are  given 
in  Table  8,  the  actual  cross-sections  are  plotted  as  the  open-circle  points  in  Fig  9 

Table  8.  Momentum-transfer  cross-sections  based  on  the  (12,6)- 
potential  for  (U++0)-collislons. 


K 

g,  km  /sec 

a  l{n  D2) 

a,  10‘'5cm: 

0.1 

3.24 

5.445 

4.61 

0.2 

4.58 

4.476 

3.79 

0.4 

6.48 

3.368 

2.85 

0.6 

7.93 

2.909 

2.46 

0.8 

9.16 

2.523 

2.14 

0.9 

9.71 

2.505 

2.12 

0.95 

9,98 

2.467 

2.09 

1.0 

10.24 

2.401 

2.03 

Table  1 

Eq.  (30) 

Table  1 

Eq.  ( 1  c) 

Col.  3  & 

Eqs.  (1c),  (28b) 
Fig.  9 

26 


The  long-dash  line  in  Fig.  9  is  computed  from  the  asymptotic  expression  in 
Eq.  (12c),  with  parameters  from  Eq.  (28a, b),  giving 

^+0)1^,.^  =2.15x10-"  g'vl  =  9.98x1  O'15 g^  cm2  (31a, D) 

We  will  adopt  Eq.  (3 1 )  as  the  relative-velocity  momentum-transfer  cross-section 
for  (U++0)-collisions,  based  on  the  (r~&)-potential,  for  the  velocity  range  less 
than  about  5x1 05  cm/sec.  Note  that  this  expression  hardly  differs  from  Eq.  (25) 
for  (U+0)-collisions;  the  uncertainty  in  the  calculation  does  not  warrant 
distinguishing  the  two  formulas. 

To  obtain  the  cross  section  for  the  effective  (r~6)-potential,  we  use 
Eq.  (A20)  and  the  parameters  in  Table  5  to  get 

^(U++0)|Ef((r*)-pot=  1.38x10-' V2/3  =  6.41x10-15  gZ\*  cm2.  (32a, b) 

Equation  (32),  which  is  the  same  as  Eq.  (26)  for  (U+0)-collisions,  is  also  plotted  as 
the  dotted  line  in  Fig.  9. 

2  3.3.2  [U>+(N2,02)]-Collisions.  By  using  Eq.  ( 1 3b)  for  the  polarization 
potential  and  the  polarizabilities  and  reduced  mass  numbers  from  Table  5,  we  nave 


<r(LT +N2)|polpOt=6.86x10"10^’",  =  6.86x1 0"15  cm2  (33al.a2) 
+O2)|poipot=  6.1 4 xIO-'V'- 6.14x1 0',5^ee  cm2.  (33bl,b2) 

For  the  (r~6)-potential  in  Eq.  (27a)  we  used  a  factor  of  5/4  for  a  larger 
target  size.  Similarly,  we  apply  a  factor  of  5/4  to  Eq.  (31 )  and  thereby  suggest  the 
expression 

5  rfir+OlU,,  <34a) 

=  2.69x1  O’"  ^"2/3  cm2  (34b) 

=  1.25x1 0*H  g&l"  cm2  (34c) 


Note  that  this  expression  hardly  differs  from  Eq.  (27a)  for  [U+(N2,02)]-collisions, 
the  uncertainty  in  the  calculation  does  not  warrant  distinguishing  the  two 
formulas. 
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To  obtain  the  cross  sections  for  the  effective  (r_6)-potential,  we  use 
Eq.  (A20)  and  the  parameters  in  Table  5  to  get 

-rfUVN^,,  ,r.)<01  =  1 74  x  I  O'"  =  8  08  x  1  O'15  ^  cm2  ( 35a  1 .32 ) 

»'<u*+0!1lEfi^i-p«.  =  l-59x|0''V'M  =  7'38xl0-,5<^L  cm2  (3561, 62) 
These  two  expressions  are  the  same  as  Eq.  (27b)  for  [U+(N2,02)]-collisions. 


2.3.13  Summary  Formula  for  Polarization  and  Effective  (r'VPotential.  The 
following  expression,  with  the  parameters  in  Table  5,  combines  all  the  formulas 
for  [U++(0,N2,C>2)]-collisions  based  on  the  polarization  and  effective  (r6)- 
potentials: 

T  •  11/2  i 

4*(<?)|p=2.59x1(rg  — - [1.06x1 cm2  (36a  1) 


=  2.59x10 


-IS  ak 


^m/sec 


[0106  /AAJ&hJ  <36a2) 


=>  3.84  x  1 0'1 1  [o'J12  [ A(H>) )"l/3  -  P  =  l /6  (36b  1) 

-  178x1  O'14  \tiiA,k)  <7km/SecTW ,  P  -1/6  (36b2) 


fO,  polarization  potential 
"  1/6,  effective (r-^J-potentlal 

Equation  (36b  1 )  is  the  same  as  Eq.  (A20a2). 

2.3.4  [UO++(0,N2,02)]-Collisions. 

For  [U0M0.N2, (^-collisions,  we  are  unaware  of  any  simplified  procedure 
we  might  use  to  estimate  the  cross  sections.  We  are  inclined  to  discount  possible 
use  of  the  polarization  potential  -  despite  our  presenting  formulas  based  on  it  - 
for  qualitative  reasons  similar  to  those  to  which  we  appealed  in  discounting  it  for 
(U+ ♦Obcollisions,  even  though  we  cannot  provide  numerical  estimates. 


For  the  polarization  potential,  we  have  from  Eq.  ( 1 3b)  and  parameters  in 
Table  5,  formulas  for  U0+  collisions  which  are  virtually  the  same  as  those  for 
U+  collisions  owing  to  little  difference  in  the  reduced  mass  numbers: 
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<r(U0-+0>|polM,*<r(U*+0)| 

polpot 

=>Eq.(29) 


(37a) 


(37b) 

=>Eq.(33a) 

^(U0’+02)|pO)POt  *  ^,(U*+O2)|polpoi.  (370 

=>Eq.(33b) 

For  the  (r6)-potential,  to  account  for  the  larger  size  of  U0+  compared  with 
U+,  we  will  make  the  ad  hoc  assumption  that  the  cross  sections  for  [ir*(0,N2,02)]" 
collisions  should  be  multiplied  by  5/4  to  obtain  the  cross  sections  for 
[UO++(0,N2,02)]-collisions: 

(38ai 

x£q.(3l) 

<r(UO-f(N2,0J)]U..^l  .s<rtU'+(N2.02)l|^wot  (380) 

=*  l  xEq.(34) 

For  the  effective  (i — 6)-potential,  we  use  Eq.  (A  19)  for  the  parameter 
-  since  it  depends  only  on  the  target  species  -  with  Eqs.  ( 1 2c)  and  ( 1 2e)  The 
results  for  the  [UO++(0,N2,02)Kollisions  will  differ  from  those  for 
[ir+(0,N2,02)]-collisions  only  by  a  very  slightly  different  reduced  mass  number, 
certainly  too  little  to  matter.  Hence,  we  write 

^*+°>Lr-W-  <39a) 

=>£q.(32) 

<39b) 

=*Eq.(35a) 

<7(U0  +02)jEff  (r-^j.pot  +02)|eff  (r^pot  (39c) 

=>Eq(35b) 
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2.3  5  Summary  Form  for  [(U,ir,UO+)+(0,N2,02)]-Collisions. 


For  later  use  it  is  convenient  to  write  the  summary  formula 


=  cm2 

$  =  1  km/ sec 

_fl/2,  polarization  potential 

P  1/3,  (r6)-  and  eff.fr^J-pot 


with  the  coefficients  a«p  given  in  Table  9. 

Table  9.  Summary-equation  parameters  for  [(U,U+,U0+)+(0,N2,02)]- 
collisions. 


k 

<,  I0-'5  cm2 

D 

/  =  U 

/SU* 

/  =  U0+ 

Potential 

0 

5.96 

5.96 

pol. 

1/2 

9.79 

9.98 

12.5 

r-6 

1/3 

641 

6.41 

6.41 

eff.  (r-6) 

1/3 

n2 

6.86 

6.86 

pol. 

1/2 

12.4 

12.5 

15.6 

~  r-6 

1/3 

808 

8.08 

8.08 

eff.  (r_&) 

1/3 

02 

6.14 

6.14 

pol 

1/2 

12.4 

12.5 

15.6 

~  r-6 

1/3 

7.38 

7.38 

738 

eff.  (r-6) 

1/3 
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SECTION  3 


STOPPING  POWER,  MEAN  MOMENTUM-TRANSFER  CROSS-SECTIONS, 
AND  RANGE-ENERGY  RELATIONS 


3.1  INTRODUCTION. 

The  goal  of  this  Section  is  to  obtain  simple  range-energy  relations  for  low- 
energy  uranium  species  traversing  the  ambient  atmosphere.  Such  relations  are 
obtained  by  integrating  an  appropriate  stopping-power  formula  involving  mean 
values  of  the  momentum-transfer  cross-sections  developed  in  Section  2. 

To  obtain  a  stopping-power  formula  for  a  single  particle  traversing  a 
Maxwellian  gas,  we  will  adapt  ideas  from  Banks  [Ref.  6]  who  was  interested  in  the 
related  problem  of  elastic  collisions  and  energy  transfer  between  gases  which 
have  separate  Maxwellian  velocity  distributions.  Banks  [Ref.  6]  shows  that  the 
exact  equation,  derived  by  Deslogue  [Ref.  1 7]  and  presented  below  as  Eq  (41 ),  for 
the  energy  transfer  rate  between  gases  of  arbitrary  temperature  and  particle 
mass,  can  be  separated  into  three  fundamental  factors,  each  of  which  depends  on  a 
different  aspect  of  the  collision  process  and  gas  composition,  i.e.,  (Da  ratio  of 
particle  masses,  (2)  the  difference  in  the  gas  thermal  energies,  and  (3)  a 
nonequilibrium  collision  frequency  for  energy  transfer  which  entails  a  mean 
nonequilibrium  momentum-transfer  cross-section. 

Banks  [Ref.  6]  acknowledges  that  his  decomposition  of  Deslogue's  formula 
[Ref.  17]  is  not  rigorous  since  he  did  not  use  the  proper  averaging  techniques  of 
kinetic  theory  needed  to  arrive  at  an  exact  expression.  Whatever  shortcoming  that 
fact  imposes  probably  is  not  very  important  relative  to  (a)  our  additional 
approximations  and  (b)  the  uncertainty  in  the  basic  cross  sections  used  in  our 
applications.  To  provide  insights  concerning  the  numerous  approximations 
required  to  obtain  simple  analytical  expressions,  we  will  first  sketch  Banks' 
development  in  general  terms,  followed  by  more  details  with  equations  and  our 
adapted  formulas. 

3.2  ENERGY  EXCHANGE  BETWEEN  MAXWELLIAN  GASES. 


3.2. 1  Desloge's  Exact  Formula. 


Desloge's  [Ref.  1 7]  result  for  the  average  rate  of  change  of 
energy  of  one  gas  is 


m  -  _Anan  _ 

dt  ~  q  1{rr\^m7)\2xk)^m^T^rr\T2fn 


f  /<7*^)exp  (-Kg2) 

•  0 


the  total  kinetic 
ag  (4ia) 
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where 


and 


4- 


2k  T,  2 kT2 

A  —  - + - 

m]  m2 

q0{g)  =  2/r  JV(<7  ,0)  (1-cos0)sin0d0 


(41b) 
(410 
(41  d) 


U- gas  total  kinetic  energy 
/7-particle  number  density 
/77-particle  mass 
/■-Maxwellian  temperature 
^-Boltzmann's  constant 
^-relative  velocity  between  particles 
<7^)-velocity-dependent  momentum- transfer  cross-section 
(/-particle  velocity  in  laboratory  system 

dV-velocity-space  volume  element 
tf-center-of-mass  scattering  angle 
<7(47,  ^-differential  scattering  cross-section 
/-velocity  distribution  function. 


3.2.2  Banks'  Heuristic  Formula. 


Banks'  development  is  based  initially  on  a  simple  model  of  energy  transfer 
for  a  single  particle  moving  in  a  gas  (exactly  our  case  of  interest!),  but  that 
particle  soon  becomes  a  member  of  a  Maxwellian  gas  interacting  with  the  original 
gas.  By  defining  an  appropriate  collision  frequency  for  energy  transfer  and  a 
momentum-transfer  cross-section,  Banks  derives  a  general  functional  form  for 
the  energy  exchange  rate  between  two  gases.  By  comparing  that  form  with  Eq. 
(41a),  Banks  obtains  specific  equations  for  the  collision  frequency  and  momentum 
transfer  cross-section  which  are  applicable  to  the  problem  of  energy  transfer. 


The  average  loss  of  kinetic  energy  per  collision,  Jr,,  of  a  single  particle  of 
mass  n\  and  kinetic  energy  p  traveling  through  a  gas  composed  of  particles  of 
mass  /?%  and  average  energy  e2  is,  as  shown  (for  isotropic  scattering  in  the  center 
of-rnass  system)  by  Huxley  and  Crompton  [Ref.  35], 

2mxml 


^  (m. 


■mkf 


(*i-f>). 


(42) 


To  describe  the  rate  at  which  the  single  particle  loses  energy,  Banks 
introduces  the  concept  of  the  single-particle  collision  frequency,  given  by 
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V\2  =  fhg<7o,  (43) 

with  Hi  the  ambient  gas  number  density.  To  obtain  the  average  rate  at  which  the 
single  particle  loses  energy,  Banks  combines  Eqs.  (42)  and  (43): 


Jf,_  2/77, 17^ 

At  {rr^+m^2  '  ’ 


7‘  v\7 


(44) 


Next,  instead  of  a  single  particle,  Banks  considers  a  large  number  of  particles 
forming  a  Maxwellian  gas  mixed  with  the  original  gas,  and  approximates  the  total 
average  energy  exchange  rate  by  the  expression 


with 


Irgrr^ 

d<f  "  (/^  +  /77j)2 


/?,  (i",  -  C2)  Vj2 


Uy=r\{m,\^ 


(45) 

(46) 


where  is  now  an  average  energy  which  corresponds  to  the  Maxwellian 
distribution  of  single  particles  permitted  to  become  the  mixed  gas.  Similarly,  yt2 
is  now  an  average  collision  frequency  which,  unlike  Eq.  (43),  must  be  suitably 
defined  to  account  for  the  many  different  relative  velocities  between  the  various 
gas  particles. 


Banks  states  that  Eq.  (45),  while  not  rigorous,  provides  a  functional  form 
which  can  be  used  to  decompose  Eq.  (41 )  into  three  factors:  (Da  ratio  of  masses, 
(2)  a  difference  in  average  particle  energies,  and  (3)  an  energy-transfer  collision 
frequency.  Only  the  third  factor  involves  the  interparticle  forces. 

Banks  argues  that,  in  order  that  the  correct  form  may  be  synthesized  by 
comparing  Eqs.  (41a)  and  (45),  the  functional  form  of  the  average  collision 
frequency  must  be 

V\2~n2gQ0  (47) 

which  allows  for  an  arbitrary  numerical  factor  in  the  final  result  for  vn  The 
quantity  g  is  the  average  relative  velocity  between  the  particles  of  the  two 
Maxwellian  gases, 

|cP^dV2  (48a) 


Q/r1 

1/2 

'RM 

X 

jt\  m2 

while  Q0  is  the  average  momentum-transfer  cross-section  appropriate  for  thermal 
nonequilibrium, 
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°0  =  Kl  f  g5  Q0(g)  exp  {-Kg2)  d  g  (49) 

•*0 

with  K  given  by  Eq.  (41c).  By  comparing  Eqs.  (41a),  (45),  (47),  and  (49),  Banks 
recognizes  the  average  momentum-transfer  collision  frequency, 

^12  =  3^2  g  >  (50a) 

with  g  given  by  Eq.  (48b),  and  Q0  by  Eq.  (49),  as  the  generalization  to  conditions  of 
different  Maxwellian  velocity  distributions  of  two  gases.  In  a  later  reference 
[Ref  7,  Pt.  A,  pp.  190,  191],  Banks  observes  that  this  collision  frequency  is  in  the 
center-of-mass  system  and  related  to  that  in  the  laboratory  system  by 

vin=[m2/{m^m2)\vX2.  (50b) 

Finally,  as  the  goal  of  his  derivation,  Banks  expresses  Eq.  (41  a)  in  terms  of  his 
decomposition  factors  as 


d  u  ram* 

or,  alternatively,  as  given  in  Ref.  7,  Pt.  A,  p.  1 9 1 ,  as 
d£/.  m. 


(5 1  a) 


(51b) 


3.3  PARTICLE  TRAVERSING  A  MAXWELLIAN  GAS. 
3.3.1  Energy  Loss  Rate. 


For  our  applications,  we  want  to  retain  the  identity  of  a  single  particle 
traversing  the  ambient  Maxwellian  gas.  Our  equation  corresponding  to  Eq.  (42)  will 
be  written  as _ 

SE,k(E,)  =  A,k  (£,-£,)  (42*) 

where  the  ambient  energy  per  particle  and  temperature  are  related  by 


and  the  mean  fractional  loss  in  laboratory  energy  is 
2  M,Mk  2/1(4,,) 

'*  *  (M,  +  Vk?  *  Af+Ak 
with  the  reduced  mass  number 

A.  Ak 
fl{A,k)~~A~7Ak 

Our  expressions  corresponding  to  Eqs.  (43)  and  (44)  are 

Kt*  =  4  g  Qotk 

and 


(42a*) 


(42b*) 


(42c*) 


(43*) 


(44*) 
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For  our  use  we  will  replace  v,k  in  Eq.  (44*)  by  y,Awith  a  functional  form 
corresponding  to  Banks'  expression,  Eq.  (50),  for  the  mean  collision  frequency, 
v,^\nkg,kOD>k  (50a*) 

vUk  =  [(/V2 !{t1s  +Mj\v,k  ,  (50b*) 

but  we  will  associate  a  fictitious  temperature  T,ua  with  the  incident  particle  by 
the  expression 

=  §  kT'Kt  (52) 

so  that  our  mean  relative  velocity  corresponding  to  Eq.  (48b)  is 


_  [8*  ]\l/,2  ra 

9tk=  n  Z  k* Mk 


o  f2  . 

■  s  1^1 


±]“\2£L  ILV 

Zn  n  +  ML 


(48b!*) 


(48b2*) 


(48b3*) 


[3*J  Mk\  ■ 

In  Appendix  E  we  show  that,  for  our  applications,  Eq.  (48b*)  for  the  mean  relative 
velocity  compares  very  well  with  an  expression  based  on  the  exact  relative 
velocity  distribution  for  a  particle  traversing  a  Maxwellian  gas. 

Our  mean  momentum-transfer  cross-section  corresponding  to  Eq.  (49)  is 
°D,k  =  Kl\o  /^,^)exp(-/r/)d^  (49*) 

with  the  parameter  K  corresponding  to  Eq.  (41c)  now  given  by 


K  = 


3  ' 


(41 c I*) 


(41c2*) 


2 

Thus,  our  expression  corresponding  to  Banks'  Eq.  (51a),  except  that  it  pertains  to  a 
single  incident  particle  instead  of  an  entire  gas,  is 

(5la<) 

with  the  auxiliary  Eqs.  (50a*),  (48b*),  (49*),  and  (41c*). 

3.3.2  Mean  Momentum-Transfer  Cross-Sections. 

Before  introducing  the  stopping  power,  we  specialize  Eq.  (49*)  to  the 
momentum-transfer  cross-sections  in  Section  2  for  [(U,  U+,  UO+) +  (0,  N2,  O2)]- 
coliisions.  Use  of  <r,k(g),  given  by  Eq.  (40),  for  g0,k{g)  in  Eq.  (49*),  gives 

Qoik  -  ,  (I  O5)20  K*  J. ,  (53) 
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where  Jf 'p  is  defined  in  terms  of  Euler's  integral  for  the  gamma  function  [Ref.  1, 
p.  255], 

J0  =  £  xz'p~'  e~x  dx  ( 54a ) 


=  /'(3-p)  (54P) 

>(8 /  3)  =  [(2 x  5) / 32]  r(2 /3)  =  1.5046,  p  =  1  / 3 
=  (54c) 

/'(5/2)  =  [(1x3)/22]/'(1/2)  =  1.3293,  p  =  l/2 


and  from  Eqs.  (41c*)  and  (42a*), 


(10'Vf  = 


ip 


10 


10 


£? 


IP 


E I  +  E,  k 


with 


(55a) 

(55b) 


E*  =  \M,  (105)2  (56a) 

=  5.1 8x1  O'3 A,,  eV  (56b) 

E,k  =Ea  (fy  !Mk)  =  Ea  ( A,/Ak\  (57) 


3.3.3  Stopping  Power. 

To  obtain  a  stopping-power  formula  we  divide  dE,/dt  by  the  incident 
(laboratory)  velocity  V,  to  obtain  the  energy  transfer  per  unit  path  length  (the 
stopping  power)  and  sum  the  stopping-power  terms  for  the  several  species  in  the 
ambient  atmosphere: 


d*, 


14L..L 

K  M  K 


EEf^E, )  3  ,<7,* 


(58) 


Use  of  Eqs.  (48b2*),  (54),  and  (56b)  in  Eq.  (58)  gives 


^  -  -fp  ^  5  ^  + *r 

or.  with  use  of  Eqs.  (42*),  (52),  (48b3*),  and  (57), 

<c>p£p,y  YiAZ'ZA 


E? 


ip 


F  LJ. » 

d/?,  "p  E)n 


rp 


i-i 


with 

and  where  the  numerical  factor  is 


(59a) 

(59b) 

(59c) 
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I-P 


(59d1 ) 


We  now  make  two  related  approximations  in  the  factor 

<59e) 

appearing  in  the  stopping-power  expression  in  Eq.  (59b).  First,  owing  to  ( I )  the 
weak  dependence  (for  p  =  1/3)  (there  is  no  dependence  for  p  =  1/2)  of  the  stooping 
power  on  the  factor  F{E,,E,k)  and  (2)  our  desire  to  get  a  relatively  simple  range- 
energy  formula,  we  will  eliminate  the  dependence  on  Ak  in  the  factor  by  using  an 
average  value  of  Ak,  Ak,  with  plausible  but  nonunique  weighting  factors,  defined  as 

<60) 

with 

(4°T  =  £4,k;T'  <6I) 

Then,  from  Eq.  (57),  we  can  write 

<62) 


and  our  stopping-power  formula  as 


-fp  [e;]>  uyr 


(63) 


Before  assessing  the  consequence  of  using  the  first  approximation 
expressed  in  Eqs.  (60)  and  (62)  for  p  =  1/3,  we  introduce  some  necessary 
parameters  for  current  and  later  use.  Table  10  contains  selected  properties 
(provided  by  I.  L.  Kofsky)  at  four  altitudes  in  a  MSIS-83  atmosphere  [Ref  29]  for 
low-latitude  (33  deg),  low  solar  10.7-cm  flux  (78),  and  magnetically-quiet 
conditions  (Ap  =  3).  Table  1 1  provides  mass  parameters  for  [OJ,  U+,  UO+) +  (0,  N2, 
02)]~col  1  isions.  Table  12  provides  the  rms  thermal  velocities  of  U,  U*,  and  U0+  for 
the  temperatures  in  Table  10.  For  the  collision  pairs  in  Table  1 1 ,  stopping-power 
and  range-energy  parameters  are  given  in  Tables  13a  and  13b  for  the  polarization 
potential,  in  Tables  14a,  14b,  and  14c  for  the  (r-S)-potential,  and  in  Tables  1 5a 
and  15b  for  the  effective  (r6)-potentia1,  corresponding  to  the  momentum- 
transfer  cross-sections  summarized  in  Table  9. 
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Table  10.  Selected  properties  for  a  quiet  MSIS-83  atmosphere  [Ref  291 


7,  km 

T,  K  ^"a,eV 

cm'3 

[0] 

[NJ 

[OJ 

200 

752.6  0.0973 

2.30  E+09 

2.36  E+09 

1.22  E+08 

225 

778.7  0.1007 

1  23  E+09 

8.33  E+08 

3.74  E+07 

250 

792.7  0.1025 

6.82  E+08 

3.09  E+08 

1.20  E+07 

275 

800.3  0.1035 

3.89  E+08 

1.18  E+08 

4.00  E+06 

cm-3 

7,  km 

/7(total) 

fio] 

f[N2] 

f[02] 

200 

4.78  E+09 

0.481 

0.494 

0.025 

225 

2. 1 0  E+09 

0.586 

0.397 

0017 

250 

1.00  E+09 

0.680 

0.308 

0.012 

275 

5.1  1  E+08 

0.761 

0.231 

0008 

Table  11  Mass  parameters  for  [(U,Lr,U0"‘) +(0,N2>02)]-col  1  isions 

/  k 

A 

4,  k 

U,  U+ 

238 

0 

16 

14.99 

0.1  18 

N2 

28 

25.05 

0.188 

02 

32 

28.21 

0.209 

UO+ 

254 

0 

16 

15.05 

0.1  1  1 

N2 

28 

25.22 

0.179 

02 

32 

28.42 

0.199 

Table  12. 

Root-mean-square  thermal  velocities  of  U,  U+,  and  UO+  in 
the  Table- 10  atmosphere. 

km/ sec 

/ 

h  =  200 

7  =  225 

7  =  250  7  =  275 

U,  U+ 
UO  + 

0.2809 

0.2719 

0.2858 

0.2766 

0.2883  0.2897 

0.2791  0  2804 

Table  13a.  Stopping-power  and  range-energy  parameters  for 
polarization  potential  in  [LT  +(0,N2,02)]-col  1  isions. 


_ 200  km _ 

_ 225  km _ 

z0,p 

'tfk 

4>  urr1 

j°.p 

/'tk 

aaaW 

i 

k 

10s  cm 

o 

i 

<* 

O 

3 

2. 

10s  cm 

1  O'6  cm"1 

U+ 

0 

0.730 

1.62 

1.36 

0.87 

N2 

0.618 

3.04 

1.75 

1  07 

02 

13.3 

0.16 

43.5 

0.05 

[40'pr’  = 

4.82 

Jx 

O 

TJ 
*  I1 

II 

1.99 

4°‘P  = 

2.07  km 

4>-p  = 

5.03  km 

f°  = 

12.7 

r°  = 

12.2 

250  km 

275  km 

i°,p 

'♦ik 

4*  u,°kpr' 

JO.P 

^ik 

4*  r1 

/ 

k 

10s  cm 

!0'7cm4 

10s  cm 

10’ 7  cm'1 

U+ 

0 

2.46 

4.80 

4.31 

2.74 

N2 

4.72 

3.98 

12.4 

1.52 

02 

136 

0J5 

407 

0.05 

U0'9]'1  = 

8.93 

WT'  = 

4.31 

40,P  = 

1 1.2  km 

4°,P  = 

23.2  km 

f°  = 

12.0 

f°  = 

119 

Table  13b  Stopping-power  and  range-energy  parameters  for 
polarization  potential  in  [UO++(0,N2,02)]-collisions. 


200  Km  _ 225 Jim 

^  rjo-Pi-'  jo,p  A 


»0,p 

*tk 

4*  wrr' 

7°.p 

^ik 

4a  ur  r1 

/ 

* 

hUH 

1 0'6  cm'1 

105cm 

1 0’6  cm’1 

UO* 

0 

0.730 

1.52 

1.36 

082 

N2 

0  618 

2.90 

1.75 

1.02 

02 

13.3 

0.15 

43.5 

0,05 

[4°-pr'  = 

4.57 

[40iPr'  = 

1.89 

40,P  = 

2.19  km 

T> 

II 

5.29  km 

13.5 

13  1 
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Table  1 3b.  Stopping-power  and  range-energy  parameters  (Cont'd) 


250  km _  _ 275  km 


JO,  P 

4*  LftT 

jOk 

Aik 

4,  Uf  ]*' 

/ 

K 

105  cm 

10* 7  cm*' 

10s  cm 

10* 7  cm'1 

uo+ 

0 

2.46 

4.51 

4.31 

2.58 

N2 

4.72 

3.79 

12.4 

1.44 

02 

136 

0.15 

407 

0.05 

M°'pr'  = 

8.45 

[40>pr'  = 

4.07 

% 

•o 

II 

1 1.8  km 

© 

*o 

II 

24.6  km 

r°  = 

12.8 

r°  = 

12.7 

Table  14a.  Stopping-power  and  range-energy  parameters  for 
(r~6)-  potential  in  [U+<0,N2,02)]-collisions. 


JO.P 

A\k 

a,  wf  r 

7°.p 

A\k 

Afk  urr1 

/ 

k 

10s  cm 

1 0’6  cm'' 

1 0*6  cm-1 

U 

0 

0.444 

2.66 

0.830 

1.42 

N2 

0.342 

5.50 

0.968 

1  94 

02 

6.61 

0,32 

21.6 

0.10 

u°'pr'  = 

8.48 

u^r1 = 

3.46 

40,P  = 

1.18km 

T) 

II 

2.89  km 

A,/Ak  = 

9.76 

Ar/Ak  = 

10.3 

_ 250  Km _  _ 275  km 

&  4>urr'  w  4A^r 

L  k  105cm  10~7cm~'  10s  cm  10'7cm~' 

U  0 

N2 
02 


1.50  7.87 

2.61  7.20 

67.2  0.31 

[4°'p]-'=  15.38 

40,P  =  6.50  km 

A,/Ak  =  10.8 


2.63 

4.49 

6.83 

2.75 

202 

0.10 

[4°'P]"  = 

7.34 

A0,P  _ 

4  = 

13.6  km 

> 

ii 

1 1.5 
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Table  1 4b  Stopping-power  and  range-energy  parameters  for 
(r-6)-  potential  in  [U++(0,N2)02)]-collisions. 


_ 225  km _ 

3°.p 

'tlk 

4>  urr’ 

7°'p 

'tfk 

4,  urr1 

/ 

k 

10s  cm 

1  O'6  cm*1 

10s  cm 

1 0'6  cm'1 

u+ 

0 

0.436 

2.71 

0.815 

1.45 

N2 

0.339 

5.55 

0.960 

1.96 

02 

6.56 

0.32 

21.4 

0.LQ 

8,58 

i4#*pr'  = 

3.51 

T3 

It 

1.17  km 

40,P  = 

2.85  km 

A, /At- 

9.77 

4/4  = 

10.3 

250  km _ 

_ 275  km _ 

»°.p 

'Hk 

4>  urr1 

JO.P 

A\k 

4*  u,f  r1 

/ 

k 

10s  cm 

o 

1 

-4 

o 

3 

10s  cm 

1 0'7  cm'1 

u+ 

0 

1.47 

8.03 

2.58 

4.57 

N2 

259 

726 

6.78 

2.77 

02 

66.7 

0.31 

200 

0.10 

U'^r1 = 

15.60 

© 

~o 
'  i* 

II 

7.44 

4°,P  = 

6.41  km 

4> 

x> 

II 

13.4  km 

4/4  = 

10.9 

4/4  = 

11.5 

II 

II 

II 

II 

II 

II 

II 

ii 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

Table  14c.  Stopping-power  and  range-energy  parameters  for 
(r~6)-  potential  in  [UO%(0;N2,02)]-collisions. 


200  km  _ 225  Km. 

^  f  JO.PVI  }0,P  . 


30. p 

''Mk 

4*  itfr 

3°>p 

•*lk 

4>  wsr 

/ 

k 

Rr  cm 

1 0*6  cm'' 

10s  cm 

1 0'6  cm'1 

U0+ 

0 

0.348 

3.19 

0.650 

1.71 

N2 

0.272 

6.58 

0.770 

2.32 

O2 

5.25 

0.38 

17.1 

0.12 

[4°-pr'  = 

10.15 

uo,p  r' = 

4.15 

40,P  = 

0,985  km 

40,P  = 

2.41  km 

4/4  = 

10.4 

4/4  = 

1 1.0 
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Table  14c.  Stopping-power  and  range-energy  parameters.  (Cont'd) 


250  km _  _ 275  km 


?0.p 

'Hk 

4*  urr1 

*ik 

4k  urr1 

/ 

k 

10s  cm 

10" 7  cm'1 

10s  cm 

0 

» 

-si 

0 

3 

uo+ 

0 

1.17 

9.49 

2.06 

5.39 

N2 

2.07 

8.65 

5.43 

3.30 

02 

53.4 

0.37 

160 

0.12 

[4‘PT'  = 

18.51 

[40‘pr'  = 

8.81 

40,P  = 

5.40  km 

jO,P  _ 

A  = 

1 1.4  km 

> 

n 

1 1.6 

4/4  = 

12.3 

Table  15a.  Stopping-power  and  range-energy  parameters  for 

effective  (r6)-  potential  in  [(U,Lf)+(0,N2,02)]-collisions, 


200  km _  _ 225  km 


jO,P 

"ik 

4k  urr' 

jo.p 

*ik 

4k  [4°k'P]"' 

/ 

k 

105  cm 

1 0‘6  cm'1 

105  cm 

1 0"6  cm'1 

U,  U+ 

0 

0.678 

1.74 

1.27 

0.93 

N2 

0.524 

3.59 

1.49 

1.26 

O2 

1  1.1 

0.19 

362 

0.06 

uo'pr'  = 

5.52 

[40iPr'  = 

2.25 

40,P  = 

1.81  km 

40,P  = 

4.44  km 

11 

K* 

s 

9.77 

4/4  = 

10.3 

_ 250  km _ 

275  km 

4k P 

4k  u?r' 

4k P 

4k  ufr’ 

/ 

k 

10" 7  cm'1 

ft--'  ISi 

10" 7  cm'1 

u,  u+ 

0 

2.29 

5.15 

4.01 

2.94 

N2 

401 

4.69 

10.5 

1.79 

O2 

1 13 

0.18 

339 

0.06 

[4°-pr'  = 

10.02 

[4°’pr'  = 

4.79 

40,P  = 

9.98  km 

40,P  = 

20.9  km 

4/4  = 

10.9 

> 

n 

1 1.5 

Table  1 5b.  Stopping-power  and  range-energy  parameters  for 

effective  (r6)-  potential  in  [UO++(0,N2,02)]-col1isions. 


200  km 

225  km 

jo,p 

'tik 

4>  urr' 

Ak  urr' 

/' 

k 

I05  cm 

1 O*6  cm'1 

l  O5  cm  1  O'6  cm'1 

uo+ 

0 

0.678 

1.64 

1.27  0.88 

N2 

0,524 

3.42 

1.49  1.20 

02 

111 

0.18 

36.2  005 

[4°'p  r1  = 

5.24 

u°-pr'  =  213 

40,P  = 

1.91  km 

4°'p  =  4.69  km 

A,/Ak  = 

10.4 

o 

ii 

s. 

250  km _ 

275  km 

i°,p 

<*ik 

Ak  urr 

Ak  u,y  r1 

/ 

k 

105  cm 

1  O’7  cm-1 

1 0s  cm  1 0'7  cnrf 1 

uo+ 

0 

2.29 

4.85 

4.01  2.27 

N2 

4.01 

4.46 

10.5  170 

02 

1  13 

0.18  . 

339  0  06 

= 

9.49 

[4°'p]-'  =  4.53 

4°'P  = 

10.5  km 

40lP  =  22  1  km 

A,/Ak  = 

1 1.6 

A,/Ak  =  12.3 

To  proceed  now  with  the  assessment  of  the  first  approximation,  we  use  the 
information  in  Table  15a  to  evaluate  the  following  two  expressions  for  [(u,  U+)  * 
(0,  N2,  02)]-collisions,  based  on  the  effective  (r_6)-potential,  at  250-km  altitude 
in  the  atmosphere  given  in  Table  10. 

(exact)  £  [e+  A,  /  AkV*  =  (64a) 

(approx. )[£■+  A, /  Akr  £  4* [^pf '  =  1.588x1  O'6  (64b) 

with  the  reduced  energy 

c*£,/£a.  (64c) 

For  r=  5,  the  error  is  only  0.69%,  which  is  certainly  acceptable. 
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The  second  approximation  is  to  replace  the  slowly-varying  factor  E(E,,E/k) 
by  a  constant  factor,  F{£itE,k),  where  the  average  energy  E,  is  defined  somewhat 
arbitrarily  as 


E^'^+E)") 

(65a) 

Ejs)  =  starting  energy 

(65b) 

?Pm finishing  energy. 

(65c) 

In  support  of  the  second  approximation,  note  that  the  variation  of  the  factor 
f{E,,E,k)  is  less  than  about  10  percent  over  the  energy  range  of  principal  interest, 
as  seen  in  Table  16  where  we  use  a  nominal  value  of  10  for  the  ratio  A,  !Ak. 


Table  16.  Illustration  of  a  slowly-varying  factor  in  the 
stopping-power  formula  for  p  =  1/3. 


1+10 
€  +  10 


r  = 

-i  1/6 


3  5  10 

0.973  0.950  0.905 


With  this  second  approximation,  our  stopping-power  formula  in  Eq.  (63) 
becomes 


*Ln  r)  Ll L 

d  R,  £-*  EY1 


eV/km 


with 


^E.i.p  = 


f.  t  vw+w 


(66a) 


(66b  1 ) 


Fp  [f°  ]p  [t  +A,  /4j*'p 

*-F,/E. 

=  \{es+cf) 

S-C/E, 

-5.1 8x1  O'3  A,  /f„(eV) 
cs  =  e°[y,{  km/ sec)]* 

Values  of  e°  are  given  in  Tables  1 3a  and  1 3b. 


(66b2) 

(66c 1) 
(66c2) 
(66d1 ) 
(66d2) 
(66e) 
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3.3.4  Range-Energy  Relation. 


If  we  Integrate  Eq.  (66a)  from  start  (s)  to  finish  (f)  (and  remember  that  /? 
is  the  distance  traversed,  not  the  residual  range),  we  have 


£-1/2 

R,  8icf,cs)  km 
with  the  reduced  energy  given  by  Eq.  (64c)  and  the  coefficient  CtSp  by 

0,u>  =  E*2  Ctw  km,  (67c) 

collated  in  Table  17  for  U,  U+,  and  UO+  at  the  altitudes  and  atmosphere  in 
Table  10. 

The  integral  in  Eq.  (67)  may  be  expressed  and  evaluated  [Ref.  13,  p.  63, 
No.  92]  as 

f,  /»V2  _ 

G(ef,es)*\  — :  tic  (67dl) 

==7[+/r2  (67d2) 


(67a) 

(67b) 


(6703) 


4  ■  ln| 


r±-\  eu,2+\ 


42+ir;2-i 


1/2 


for  cf  >  1 


(67d4) 


If  one  ignored  the  ambient  energy  in  the  denominator  of  the  integrand,  the 
integrand  would  be  rU7  and  only  the  term  4  would  obtain.  But  with  retention  of 
the  ambient  energy  in  the  denominator,  we  have  a  logarithmical  ly-divergent 
integral  for  cf-\.  This  means,  of  course,  that  our  formulation  of  the  energy-loss 
process  is  inappropriate  for  obtaining  a  range-energy  relation  if  one  insists  on 
slowing  the  particle  all  the  way  to  ambient  energy.  We  shall  not  so  insist!  We 
recognize  that  nature  knows  how  to  thermalize  an  energetic  particle.  We  expect 
that  as  the  particle  slows  to  the  point  where  it  is  a  viable  candidate  for  inclusion 
as  a  member  of  the  thermal  distribution,  one  should  probably  be  concerned  with 
some  sort  of  a  random  walk  or  diffusion  problem  (for  which  projected-motion 
questions  could  be  investigated)  instead  of  a  (linear)  range-energy  problem  But 
investigating  that  aspect  of  the  problem  is  beyond  the  scope  of  the  current 
considerations. 
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Table  1 7.  Range-energy  coefficients  Cr jp  for  U,  U+,  and  UO+  at  several 
altitudes  in  a  quiet  MSIS-83  atmosphere. 


C t,\f>  i  km 


Potential 

/ 

h  =  20Q_ 

h  =  225 

h  =  250 

h  =  275 

Polariz. 

U+ 

0.581 

1.44 

3.23 

6.73 

uo+ 

1.46 

3.30 

690 

(u*.  uo-) 

0.589 

1.45 

3.27 

6.82 

(r-6) 

u,  u* 

0.315 

0.776 

1.74 

3  63 

c 

o 

+ 

0.255 

0.625 

1.40 

2.95 

Eff.  (r-6) 

u,  u+ 

0.482 

1.19 

2.68 

558 

uo+ 

0.492 

1.22 

zn 

5J1 

(u.u'.uo*) 

0.487 

1.21 

2.71 

5.65 

Notes: 

1 .  For  the  polarization  potential,  the  difference  between  the  U+  and  UO+ 
parameters  here  (and  in  Tables  13a  and  13b)  is  due  to  slightly  different 
reduced  masses,  maintenance  of  which  is  not  justified  since  it  was 
ignored  in  the  summary  cross-section  parameters  in  Table  9  in  Section  2 

2.  To  obtain  numerical  values  for  the  (r~&)-  and  effective  (r-&)-potentials, 
we  have  evaluated  the  weakly  energy-dependent  factor  involving  e  [Eq. 
(66c2)l  by  using  a  value  of  0.57  km/sec  (the  rms  velocity  for  a  3 1 00-K  U 
source)  for  V,  in  Eq.  (66e)  for  es  and  a  value  of  1 .0  for  ef  in  Eq.  (66c2). 
Some  applications  may  require  re-setting  these  values. 

3.  For  the  (r~6)-potential, 

(a)  the  formal  distinction  between  U  and  U+  parameters  in  Tables  1 4a 
and  14b  is  ignored  here  by  using  only  the  U  parameters  from  Table 
1 4a,  and 

(b)  the  difference  between  the  UO+  and  U+  results  is  due  almost 
entirely  to  the  ad  hoc  factor  of  5/4  introduced  in  Eq.  (38)  and  only 
slightly  to  the  reduced-mass  difference. 

4.  For  the  effective  (i — 6)-potential,  the  difference  between  the  (U,  u+)  and 
UO+  parameters  here  (and  in  Tables  1 5a  and  1 5b)  is  due  to  slightly 
different  reduced  masses,  maintenance  of  which  is  not  justified  since  it 
was  ignored  in  the  summary  cross-section  parameters  in  Table  9  in 
Section  2. 
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To  gain  some  insight  into  the  behavior  of  the  integral  in  Eq.  (67di ),  we  nave 
recorded  In  Table  18  the  results  of  evaluating  the  terms  4,  T2>  the  sum  4+4,  and 
the  ratio  (4+42)/7;  for  starting  reduced-energies  es  of  3,  5,  and  10.  For  some 
applications  where  a  very  simple  expression  for  the  integral  is  desired  -  and  one 
can  accept  the  corresponding  approximation  -  one  may  choose  to  use  only  the  term 
T}  with  €f»  1.  In  this  regard  it  is  of  interest  to  note  from  Table  18  that 

7lef-\]~{T\+r2)[ef-2].  (68) 

The  significance  of  cf  equalling  2  (i.e.,  that  it  corresponds  approximately  to  the 

89th  percentile  of  the  integral  Maxwellian  energy  distribution)  can  be  seen  from 
Table  19  which  provides  a  selected  set  of  fractional  values  of  the  integral 
Maxwellian  energy  distribution,  F(ef),  with  energies  not  exceeding  ef,  as  derived 
from  Appendix  F  giving  equations  and  graphs  for  the  differential  and  integral 
distribution  functions.  Table  19  may  also  aid  one  in  determining  where  the  cutoff 
might  be  set  for  the  integral  in  Eq.  (67). 

Table  18.  Comparison  of  approximate  reduced  range-energy  integrals. 


cs  =  3  es  =  5 _  _ £>  =  10 


4+7?  4+  T2  7\+ 4 


Cf 

4 

t2  4+4 

_4_  . 

4 

4 

4+42 

4  4 

4 

4+4 

;  4 

1.0 

1.46 

2.47 

4  32 

1.1 

1  37 

2.42 

3.79 

2.77 

2.37 

2.77 

5.15 

2.17 

4.23 

3.08 

7  31 

1  73 

1.2 

1.27 

1.77 

3.05 

2.39 

2.28 

2.13 

4.41 

1.93 

4.13 

2.43 

6.57 

1.59 

1.3 

1.18 

1.41 

2.59 

2.19 

2.19 

1.76 

3.96 

1.80 

4.04 

2.07 

6.12 

1  51 

1.4 

1.10 

1.16 

2.26 

2.06 

2.1 1 

1.52 

3.62 

1.72 

3.96 

1.82 

5.78 

1.46 

1.5 

1.01 

0.98 

1.99 

1.96 

2.02 

1.33 

3.35 

1.66 

3.88 

1.64 

5.51 

1.42 

1.6 

0.93 

0.83 

1.76 

1.89 

1.94 

1.18 

3.13 

1.61 

3.79 

1.49 

5.29 

1.39 

1.7 

0.86 

0.71 

1.57 

1.83 

1.86 

1.06 

2.93 

1.57 

3.72 

1.37 

5.09 

1.37 

1.8 

0.78 

0  61 

1.39 

1.78 

1.79 

0.96 

2.75 

1.54 

3.64 

1.27 

4.91 

1  35 

1.9 

0.71 

0.52 

1.23 

1.74 

1.72 

0.88 

2,59 

1.51 

3.57 

1.18 

4.75 

1.33 

2.0 

0.64 

0.45 

1.08 

1.70 

1.64 

0.80 

2.44 

1.49 

3.50 

1.1  1 

4.60 

1.32 

2.1 

0.57 

0.38 

0.95 

1.67 

1.57 

0.73 

2.31 

1.47 

343 

1.04 

4  47 

1  30 

2.2 

0.50 

0.32 

0.82 

1.64 

1.51 

0.67 

2.18 

1.45 

3.36 

0.98 

4.34 

1.29 

2.3 

0.43 

0.27 

070 

1.62 

1.44 

0.62 

2.06 

1.43 

3.29 

0.93 

4.22 

1.28 

2.4 

0.37 

0.22 

0.58 

1.60 

1.37 

0.57 

1.95 

1.42 

3.23 

0.88 

4.1  1 

1.27 

2.5 

0.30 

0.17 

0.48 

1.58 

1.31 

0.53 

1.84 

1.40 

3.16 

0.84 

4.00 

1.26 

2.6 

0.24 

0.13 

0.37 

1.56 

1.25 

0.49 

1.74 

1.39 

3.10 

0.80 

3  90 

1.26 

2.7 

0.18 

0.10 

0.27 

1.54 

1.19 

0.45 

1.64 

1.38 

3.04 

0.76 

3.80 

1.25 
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Table  18.  Comparison  of  approximate  reduced  range-energy  integrals.  (Cont'd) 


e,  =  3  es  =  5  JV  =  1 0 


€f 

4 

4  . 

4+4 

4+4 

4 

4 

4 

4+7* 

4+4  4  4 

r2 

4+4 

4+4 

4 

2.8 

0.12 

0.06 

0.18 

1,53 

1.13 

0.42 

1.54 

1.37 

2.98 

0.72 

3.70 

1.24 

2.9 

0.06 

0.03 

0.09 

1.51 

1.07 

0.38 

1.45 

1.36 

2.92 

0.69 

3.61 

1.24 

3.0 

0.00 

0,00 

0.00 

1.01 

0.35 

1.36 

1.35 

2.86 

0.66 

3.52 

1.23 

3.2 

0.89 

0.30 

1.19 

1.34 

2.75 

0.61 

3.35 

1.22 

3.4 

0.78 

0.25 

1.04 

1.32 

2.64 

0.56 

3  20 

1.21 

3.6 

0.68 

0.21 

0.89 

1.31 

2.53 

0.52 

3.05 

1.20 

3.8 

0.57 

0.17 

0.74 

1.30 

2.43 

0.48 

2.90 

1  20 

4.0 

0.47 

0.14 

0.61 

1.29 

2.32 

0.44 

277 

1  19 

4.5 

0.23 

0.06 

0.29 

1.27 

2.08 

0.37 

2.45 

1.18 

5.0 

0.00 

0.00 

0.00 

1.85 

0.3 1 

2.16 

1.17 

5.5 

1.63 

0.26 

1.89 

1.16 

6.0 

1.43 

0.21 

1  64 

1  15 

6.5 

1.23 

0.17 

1.40 

1  14 

7.0 

1.03 

0.14 

1.17 

1.14 

75 

0.85 

0.1  1 

096 

1  03 

8.0 

0.67 

0.08 

0.75 

1.03 

8.5 

0.49 

0.06 

0.55 

1.12 

9.0 

0.32 

0.04 

0  36 

1.12 

9.5 

0.16 

0.02 

0.18 

1.1  1 

10.0 

0.00 

0.00 

000 

Table  19.  Integral  fraction  of  Maxwellian  energy  distribution. 


€f 

Fief) 

€f 

Fief) 

cf 

Fie,) 

€f 

Fief) 

1.0 

0.608 

1.5 

0.788 

2.0 

0.888 

2.5 

0.942 

1.1 

0.652 

1.6 

0.813 

2.1 

0.902 

2.6 

0.950 

1.2 

0.692 

1.7 

0.835 

2.2 

0.914 

2.7 

0.956 

1  3 

0.728 

1.8 

0855 

2.3 

0925 

2.8 

0.962 

14 

0.759 

1.9 

0.873 

2.4 

0.934 

2.9 

0.966 

3.0  0.971 
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3.3.5  Energy-Time  Relation. 


From  the  basic  relationship  between  energy-loss  rate  and  stopping  power, 
Eq.  (58),  we  can  write 


tit  = 


tiE, 


K  tiE,-/  tifi, 
7T  tiE, 


1 2  E,  tiE'ltiR, 

By  using  Eq.  (66a)  for  the  stopping  power,  we  have 

d£, 


tit  =  - 1 


with  the  scale  time 


r  •  Ctip  12  sec 

-C,u4ff,H  2£„) 

=o.,,i^r 


(69a) 

(69b) 


(69c) 

(69d) 

(69e 1 ) 
(69e2) 

(69e3) 


Values  of  {  )  ,  the  rms  thermal  velocity  of  U,  U+,  and  U0+,  are  given  in 

Table  12  for  the  temperatures  and  altitudes  in  Table  10  Values  of  the  scale  time 
t°,  corresponding  to  the  conditions  in  Table  17,  are  given  in  Table  20. 


Table  20.  Scale  time  t°  for  U,  U+,  and  U0+  in  energy-time  and  range-time 
relations  at  several  altitudes  in  a  quiet  MSI5-83  atmosphere. 


^,sec 


Potential 

/ 

h  =  200 

h  =  225 

h  =  250 

h  -  275 

Polariz. 

U+ 

2.07 

5.04 

11. 2 

23.2 

uo+ 

2.19 

323 

JLLS 

246 

(u*.  U0‘) 

2.13 

5.16 

1 1.5 

23.9 

(r6) 

u,  u* 

1.12 

2.72 

6.04 

12  5 

uo+ 

0.938 

2.26 

5.02 

10.5 

Eff.  (r-6) 

u,  u+ 

1.72 

4.16 

9.30 

19.3 

uo+ 

LSI 

4.41 

978 

204 

(uu'.uo') 

1.77 

4.29 

9.54 

19.9 
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Integrating  Eq.  (69d),  we  have 


or 


(70a) 

(70b) 


jiy-Mfj-Oexp 


(71) 


Equation  (71 ),  which  is  exact  relative  to  the  simplified  stopping  power  expressed 
by  Eq.  (66a),  has  been  used  to  provide  in  Table  2 1 ,  for  several  values  of  the 
reduced  starting  energy  (cs),  values  of  the  reduced  finishing  energy  ief)  as  a 
function  of  the  reduced  time  The  fact  that  infinite  time  is  required  for  c 

to  be  reduced  to  1.0  corresponds  to  the  divergent  integrand  in  Eq.  (67d  1 ). 


If  one  ignored  the  ambient  energy  in  the  denominator  of  the  integrand  of  Eq. 
(70a),  the  integrand  would  be  c'  and  the  integral  would  be 


or 


t, 


In 


£,  *1 


exp  (-tf/t°),  £f  . 


(72a) 

(72b) 


Equation  (72b)  has  been  used  to  provide  in  Table  22,  for  several  values  of  the 
reduced  starting  energy  (e3),  values  of  the  reduced  finishing  energy  ief  i  I )  as  a 
function  of  the  reduced  time 


3.3.6  Range-Time  Relation. 

Three  range-time  relations  will  be  presented. 

First,  to  obtain  range-time  information,  one  can  use  the  energy-time 
relation  given  by  Eq.  (71)  with  the  range-energy  relation,  Eq.  (67),  with  a 
reasonable  value  for  the  cutoff.  To  gain  further  insight  into  the  behavior  of  the 
integral  in  Eq.  (67dl )  when  it  is  used  for  range-time  relations  instead  of  range- 
energy  relations,  we  have  recorded  in  Table  23  the  results  of  evaluating  the  terms 
Tu  T2,  the  sum  f}+T2,  and  the  ratio  {T^T2)/Ty  as  a  function  of  the  reduced  time 
( tf/t° ),  for  starting  reduced-energies  es  of  3,  5,  and  10.  (In  this  implementation, 
the  result  of  evaluating  ef  from  Eq.  (71 )  was  used  in  Eqs.  (67d3)  and  (67d4).) 


50 


Table  21.  Reduced  finishing  energy  versus  reduced  time  for 

several  values  of  reduced  starting  energy,  per  Eq.  (71 ). 


t,  /t° 


0.0 

0.1 

0.2 

0.3 

04 


£>  =  3 


3.00 

2.81 

2.64 

2.48 

2.34 


05 

0.6 

0.7 

0.8 

0.9 


221 
2.10 
1.99 
1.90 
:  3i 


10  1.74 

1.1  1.67 

1.2  160 

1.3  1.55 

1.4  149 


1.5 

1.6 

1.7 

1.8 
1  9 


1.45 
1.40 
1.37 
1.33 
1  30 


2.0 

2.1 

2.2 

2.3 

2.4 


1  27 
1.24 
1.22 
1.20 
1  18 


2.5 

2.6 
2  7 
2.8 
2.9 
3.0 


116 
1.15 
113 
1.12 
1.1  1 
1.10 
m 


€f 


£>  =  5  £>=4  0 


5.00 

10  00 

4.62 

9.14 

4.27 

8.37 

396 

7.67 

3.68 

7  03 

342 

6.46 

3.20 

5.94 

299 

5  47 

2.80 

5.04 

2.63 

4.66 

2.47 

431 

2.33 

4.00 

2.20 

3.71 

2.09 

3.45 

1  99 

3  22 

1.89 

3.01 

1.81 

2.82 

1.73 

2.64 

1.66 

2.49 

1  60 

2.35 

1  54 

222 

1.49 

2.10 

1.44 

2.00 

1.40 

1.90 

1  36 

1  82 

1  33 

1.74 

1.30 

1.67 

1.27 

1  60 

1.24 

1.55 

1.22 

1.50 

1.20 

1.45 

mtMpwE 

mBBm 

oo 


Table  22.  Reduced  finishing  energy  versus  reduced  time  for 

several  values  of  reduced  starting  energy,  per  Eq.  (72b) 


tftt 0 

*V 

*V  =  3 

*V  =  5 

ii 

o 

0.0 

3.00 

5.00 

10.00 

0.1 

2.71 

4.52 

9.05 

0.2 

2.46 

4.09 

8.19 

0.3 

2.22 

3.70 

7.41 

0.4 

2.01 

3.35 

6.70 

0.5 

1.82 

3.03 

6.07 

0.6 

1.65 

2.74 

5.49 

07 

1  49 

2.48 

4.97 

0.8 

1.35 

2.25 

4.49 

0.9 

1.22 

2.03 

4.07 

1.0 

1.10 

1.84 

368 

1.1 

1.00 

1  66 

3.33 

1.2 

1.51 

3.01 

1.3 

1.36 

2.73 

1.4 

1.23 

2  47 

1.5 

1.12 

2.23 

1.6 

1.01 

2.02 

1  7 

1  83 

1.8 

1.65 

1.9 

1.50 

2.0 

1.35 

2.1 

1.22 

2.2 

1.11 

2.3 

1.00 
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Table  23;  Comparison  of  approximate  reduced  range-time  integrals 


C* 

=  3 

Cs 

=  5 

10 

t, 

4+7? 

4+^? 

w, 

r 

4 

r2 

4+7? 

4 

4 

t2 

4+7? 

4 

4 

r2 

4+T? 

4 

0.0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.1 

0.00 

0.06 

0.17 

1.53 

0  17 

0.05 

0.22 

1  26 

0.28 

0.03 

031 

1  12 

0.2 

0.22 

0.12 

0.34 

1.55 

0.34 

0.09 

0.43 

1.28 

0.54 

0.07 

0.6! 

1.12 

0.3 

0.31 

0.18 

050 

1,58 

0.49 

0.14 

0.63 

1.29 

0.79 

0  10 

0  89 

1.13 

0.4 

0.40 

0.25 

0.65 

1.61 

0.63 

0.19 

0.83 

1.30 

1.02 

0.14 

1  16 

1  M 

0.5 

0.49 

0.31 

0.80 

1.64 

0.77 

0.25 

1.02 

1.32 

1.24 

0.18 

1  42 

1.14 

0.6 

0.57 

0.38 

0.95 

1.67 

0.90 

0.30 

1.20 

1.34 

1.45 

0.22 

1  67 

1  15 

0.7 

0.64 

0.45 

1.09 

1.70 

1.02 

0.36 

1.37 

1.35 

1.65 

0.26 

191 

1  16 

0.8 

071 

0.52 

1.23 

1  74 

113 

0.42 

1.54 

1.37 

1.83 

0.30 

2.14 

t  i  n 
i  i  / 

0.9 

0.77 

0.60 

1.37 

1.77 

1.23 

0.48 

1.71 

1.39 

2.01 

0.35 

2.36 

1.17 

1.0 

0.83 

0.67 

1.50 

1.81 

1.33 

0.54 

1.87 

1.41 

2.17 

0.40 

2.57 

1.18 

1.1 

088 

0.75 

1.63 

1.85 

1.42 

0.61 

2.02 

1.43 

2.33 

0  44 

277 

1  19 

1.2 

0.93 

0.83 

1.76 

1.89 

1.50 

0.67 

2.17 

1.45 

2.47 

0.50 

2.97 

1  20 

1.3 

0.98 

0.91 

1.88 

1.93 

1.58 

0.74 

2.32 

1.47 

2.61 

0.55 

3.16 

1.21 

1.4 

1.02 

0.99 

2.01 

1.97 

1.65 

0.81 

2.46 

1.49 

2.74 

0.60 

3.34 

1.22 

1.5 

1.06 

1.07 

2.13 

2.01 

1.72 

0.88 

2.60 

1.51 

2.86 

0.66 

3.52 

1.23 

1.6 

1.09 

1.15 

2.25 

2.05 

1.78 

0.96 

2.74 

1.54 

2.97 

0.72 

3.69 

1.24 

1.7 

1.13 

1.24 

2.37 

2.10 

1.84 

1.03 

2.87 

1.56 

3.07 

0.78 

3.85 

1.25 

1.8 

1.16 

1.32 

2.48 

2.14 

1.89 

1.1  1 

3.00 

1.58 

3.17 

0.84 

401 

1  27 

1.9 

1.18 

1.41 

2.60 

2.19 

1.94 

1.19 

3.13 

1.61 

3.26 

0.91 

4.17 

1.28 

2.0 

1.21 

1.50 

2.71 

2.24 

1.99 

1.27 

3.25 

1.64 

3.35 

0.97 

4.32 

1.29 

2.1 

1.23 

1.59 

2.82 

2.29 

2.23 

1.35 

3.38 

1.66 

3  42 

1.04 

4.46 

1  30 

2.2 

1.25 

1.68 

2.93 

1.34 

2.07 

1.43 

3.50 

1.69 

3.50 

1.11 

4.61 

1.32 

2.3 

1.27 

1.77 

3.04 

2.39 

2.10 

1.51 

3.62 

1.72 

3.57 

1.18 

4.75 

1  33 

2.4 

1.29 

1.86 

3.15 

2.44 

2.14 

1.60 

3.74 

1.75 

3.63 

1.25 

4.88 

1.35 

2.5 

1.31 

1.95 

3.26 

2.50 

2.17 

1.68 

3,85 

1.78 

3.69 

1.33 

5.02 

1  36 

2.6 

1.32 

2.05 

3.37 

2.55 

2.19 

1.77 

3.97 

1.81 

3.74 

1.41 

5.15 

1  38 

2.7 

1.33 

2.14 

3.47 

2.60 

2.22 

1.86 

4.08 

1.84 

3.79 

1.48 

5.28 

1.39 

2,8 

1.35 

2.23 

3.58 

2.66 

2.24 

1.95 

4.19 

1.87 

3.84 

1  56 

5.40 

1  4! 

2.9 

1.36 

2.33 

3.69 

2.72 

2.26 

2.04 

4.30 

1.90 

3.88 

1.65 

5.52 

1.42 

3.0 

1.37 

2.42 

3.79 

2.77 

2.28 

2.13 

4.41 

1.93 

3.92 

1.73 

5.65 

1  44 

oo 

1.46 

2,47 

432 
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In  those  special  cases  where  the  ambient  energy  may  be  ignored  in  obtaining 
an  approximate  range-energy  relation,  such  as  is  given  by  only  term  4  in  Eq 
(67d2),  two  approximate  range-time  expressions  can  be  obtained.  Use  of  Eq  (7 1 ) 
for  cf  gives  the  first  of  these  two  approximate  expressions, 

R,  *  2  6^ ,  p  [  4/2  -  (I  +  ( ^  - 1 )  exp  ( -^  /  )  ]l/2  ].  (73) 

Equation  (73)  has  a  hybrid  character  because  Eq.  (71 )  for  eft  as  noted  above,  is  an 
exact  expression  relative  to  the  stopping  power  expressed  by  Eq.  (66a)  and 
includes  the  effect  of  the  ambient  energy. 

The  second  of  these  two  approximate  expressions  is  obtained  by  use  of  Eq 
(72b)  for  ef, 


tf,*2f,,p4'J[l-expK/2n|, 


(74a) 


provided 


cs  2  exp  it,  It) 


(74b) 
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SECTION  4 


DIFFUSION  COEFFICIENTS 
4  I  ORDINARY  DIFFUSION  COEFFICIENT 

In  Appendix  A,  motivated  in  part  by  subsequent  applications  to  Ba  release-: 
in  Appendix  D,  we  develop  ordinary  diffusion  coefficients  (D,)  based  on  the 
interaction  potentials  discussed  in  Section  2  for  momentum-transfer  cross- 
sections.  Here,  we  present  those  formulas  for  application  to  uranium  oxide  ions 
Equation  (A32),  based  on  the  polarization  potential  and  effective  (r-cu-potennm 
is 

,  3.2  U I  Of  T  krrf  „e., 

T  1  sec 

k 

fO,  polarization  potential 
.1/6,  effective  (/”6)-potentiai 

values  of  the  collision  parameters  (involving  the  reduced  mass  and 
polarizabilities)  are  given  in  Table  5  in  Section  2.  Values  of  the  atmospheric 
parameters  are  given  in  Table  10  in  Section  3  for  four  altitudes  in  a  quiet  MSIS-S: 
atmosphere  There  is  little  difference  between  the  reduced  masses  of  U*,  DC*  and 
UOt  for  collisions  with  atmospheric  species  Hence,  we  will  evaluate  Eq  (75> 
specifically  for  uo*,  but  apply  the  results  to  U+,  UO*,  and  uo:,  as  recorded  m  Taoie 
24 


Table  24  u*,  UO*,  and  uo;  diffusion  coefficients  at  several 
altitudes  in  a  quiet  MSIS-83  atmosphere 


n 

km 

D,,  km2 /sec 

UTJT, 

(\,  km2 /sec 

Polariz. 

Eff.  (/“-*) 

Polariz. 

Eff.  (r -) 

0.099 

0.083 

2.5 

0.25 

0  21 

225 

0  25 

0  21 

2.5 

0  62 

0  52 

250 

0.57 

0  48 

2.5 

1.4 

l  2 

275 

1  2 

1  0 

25 

3  0 

2  5 
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4.2  AMBI POLAR  DIFFUSION  COEFFICIENT. 

The  formula  normally  applied  as  the  coefficient  for  ambipolar  diffusion 
along  the  geomagnetic  field  (/})  is  that  given  for  the  ambipolar  diffusion 
coefficient  by  Eq.  (D4b)  in  Appendix  D, 

a**0([\  +  T*/T'  I.  (76) 

where  Df  is  the  ordinary  diffusion  coefficient  in  Eq.  (75)  and  Te  and  T,  respectively 
are  the  electron  and  ion  temperatures  in  the  atmosphere.  Such  temperatures,  not 
easily  predicted  for  a  given  model  atmosphere,  are  discussed  in  Appendix  D  with 
respect  to  applications  to  Ba  releases.  Here,  we  shall  simply  adopt  a  nominal 
value  of  1.5  for  the  ratio  Te/T( .  Values  of  the  ambipolar  diffusion  coefficient  are 
given  in  Table  24  for  four  altitudes  in  a  quiet  MSIS-83  atmosphere. 
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SECTION  5 


COMPARISONS  OF  MOMENTUM-TRANSFER  CROSS-SECTIONS 
WITH  SELECTED  URANIUM  OXIDE  REACTION  CROSS-SECTIONS, 
RATE  COEFFICIENTS  FOR  SELECTED  URANIUM  OXIDE  REACTIONS 


5.1  INTRODUCTION. 

It  Is  Interesting  to  compare  some  of  our  momentum-transfer  cross-sections 
developed  in  Section  2  with  cross  sections  for  certain  reactions,  specially  atom- 
transfer  reactions.  (The  writer  expects  momentum-transfer  cross-sections 
should  not  be  less  than  atom-transfer  cross-sections.)  Also,  as  an  auxiliary 
product  of  these  considerations,  rate  coefficients  are  derived  for  selected 
uranium  oxide  reactions  for  which  velocity-dependent  cross-sections  are  either 
known  or  assumed,  based  on  extrapolations  and/or  averaged  values  of 
measurements. 

Before  considering  specific  reactions,  we  develop  a  needed  relation  between 
a  power-law  reaction  cross-section  in  the  form 

ln 

a(g)=a0(g0)  &  (77) 

and  the  corresponding  rate  coefficient,  k.  We  start  with  the  general  relation 
be  ween  the  rate  coefficient  and  velocity-dependent  cross-section  for  a 
Maxwellian  distribution  of  velocities,  given  as,  e.g.,  by  Light,  Ross,  and  Shuler  [Ref 
42,  Eg.  (30)], 


2  jtjcT 


4*Jq  g*(7(g)ex P-f^  *9 


By  using  Eq.  (77)  and  letting  g0  be  the  mean  value  of  the  relative  velocity 
distribution,  i  e., 

r  -r-  -rf/2 


?o=g  = 


8  *T 


w  •  get 


k{T)  =  <70{g)d{T)Gn 


(79a) 


(79b) 


with  Gn,  proportional  to  Euler's  integral  for  the  gamma  function  [Ref.  I ,  p.  255], 
Gn  =[4/jrJ"/2  /’(2-/:72)  (790 

/"(2-  /?/2)=  f  x'~n'2e~*  6x,  (79d) 

Jo 


tabulated  in  Table  25  for  several  values  of  n. 
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Table  25.  Values  of  the  parameter  4  relating  rate  coefficient  and 
power-law  reaction  cross-section. 


n 

r(2-n/2) 

[4/*]"'2 

(?n  =  [4/^Jn/2  A2- n/2) 

0 

n  2)  =  i  .ooooo 

1.0000 

1 .0000 

1/2 

nil  A)  =  0  91906 

1.0622 

0.9763 

1 

n 3/2)  =  0.88623 

1.1284 

1.000 

3/2 

A5/4)  =  0.90640 

1.1986 

1.086 

5.2  ir+o2-»uo++o. 


5.2. 1  Results  of  Armentrout  and  Beauchamp  [Ref.  4], 


The  cross  section  measured  by  Armentrout  and  Beauchamp  [Ref  4]  for  the 
reaction 

ir+o2-*ucr+o  (80) 

is  well  represented  by  the  expression 

[50X1  O''6  AO/1  /-  -7  U 

,  0.24</fo1s3  43eV 


<t[U+(02,0)U0+]  = 


rw 
-c n 


(81a) 


21x1 0",e 


,  3  43 sl8eV 


CM 


The  exponents  of  the  energy  dependence  in  Eq.  (81a)  are  given  in  Ref.  4,  but  the 
coefficients  are  inferred  by  the  writer  from  reading  the  small  graphs  in  Ref  4  To 
express  the  lower-energy  portion  of  the  cross  section  in  terms  of  the  relative 
velocity,  we  first  write 

I.66xI0-2VK*)^2  „  ? 

^  "  y  =1.46x10-"/ 


Ecn  - 


2x1.6xl0-12 
=  0.146^ 

so  that  substitution  of  Eq.  (81  b2)  in  Eq.  (8 1  a  1 )  gives 


(8 1 b 1 ) 
(81b2) 


WU+(02,0)U0+]  = 


8.08x10- 

^km/sec 


,-15 


cm 


>•28^^^4.841 


(810 


The  cross  section  in  Eq.  (81c)  is  shown  in  Fig.  10  as  the  solid  line  and  extrapolated 
to  lower  velocities  as  the  short-dash  line. 


Armentrout  and  Beauchamp  [Ref.  4]  also  quote  a  rate  coefficient  at  300  K  for 
Reaction  (80), 
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Figure  10.  Comparison  of  momentum-transfer  cross-sections,  based  on  various 
interaction  potentials,  for  (ir+02)-col1isions  with  the  cross  sections 
for  the  atom-exchange  reaction  ->U0*+0. 


X'[Lr(02,0)UO^]  =  5.6xlO-,° 


cm3 

molecule  sec ' 


(82) 


obtained  by  extrapolating  their  data  to  lower  velocities  and  averaging  over  a 
Maxwellian  distribution.  We  now  show  the  essential  consistency  between  Eq  (82) 
and  our  formalism.  First,  we  rewrite  Eq.  (81c)  in  the  form  of  Eq.  (77),  i.e., 


<r[Lr(02,0)U0*]  =  1.1 8x1 0 


4.74  x  IQ4 

^cnvsee 


1/2 


so  that 


<r0  =  I.18x10'm  cm2,  g0  =  4.74x1 04  cm/ sec,  n  =  l/2 

Then,  from  Eq.  (79b),  we  have 

X'lUt(O2,0)UO+]  =  5.44xl0*,° 


crrr 


molecule  sec ' 


(83a) 
(83b, c,d) 

(83e) 


in  essential  agreement  with  the  value  quoted  by  Ref.  4,  Eq.  (82).  The  2.9% 
difference  is  probably  due  to  our  error  in  reading  the  small  graph  in  Ref.  4 
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5.2.2  Results  of  Johnsen  and  Biondi  [Ref.  37], 


Johnsen  and  Biondi  [Ref.  37]  measured  the  rate  coefficient  for  Reaction  (80) 
at  a  temperature  of  300  K  and  report 


/('[U+(02,0)U0+]  =  8.5 


+4 


xIO 


,.io 


cnf 


molecule  sec 


(84) 


We  now  infer  the  constant  cr0{g0)  [cf.  Eq.  (77)]  implied  by  the  rate  coefficient  in  Eq. 
(84).  If  we  assume  the  exponent  n  is  the  same  as  in  the  measurements  reported  by 
Ref.  4,  namely,  n  =  1/2,  then  from  Eq.  (79b)  we  have 


= 


K 

0.9763  g 


(85a) 


<7  = 


QjcT 

L  '*’/*  J 


1/2 


=  4.74x10''  cm/ sec  =  0.474  km/ sec 


8.5 


^(g)  = 


+4 


xio- 


'10 


0.9763x4.74x10'’ 


k(^)  =  l.84 


-0.22 

+0.86 


xlO'1*  cm2 


(85b) 

(85c) 

(85d) 


In  Fig.  10,  the  cross-section  value  in  Eq.  (85d)  is  plotted  at  the  relative  velocity  in 
Eq.  (85b)  and  labeled  as  Ref.  37.  The  central  point  of  this  inferred  cross-section 
value  is  larger  than  the  extrapolated  cross-section  of  Ref.  4  by  a  factor  of  *  1  6 

5.2.3  Momentum-Transfer  Cross-Sections  for  LF+Oj. 

From  Eq.  (40)  and  Table  9  in  Section  2  (or  equivalently,  Eqs.  (33b2),  (34c), 
and  (35b2)),  we  have  three  momentum-transfer  cross-sections  for  the  collision 
pair  IT+Oj, 


«r(U*  + 02  )j poipot  =  6. 1 4  x  1 0" 1 5 cm2 

(86a) 

<r[U++02^rl).pet=!.25xl0-H  g^  cm2 

(86b) 

^+°2^ff(r-‘H>ot  =738xlO-'5<^L  cm2. 

(860 

These  cross  sections  are  plotted  and  respectively  labeled  in  Fig.  10  as  "Pol.  pot  ", 
"~(r'&)-pot.",  and  "Eff.  (r-6)-pot."  The  slope  of  the  cross-section  based  on  the 
polarization  potential  certainly  disagrees  with  that  measured  by  Ref.  4;  the  slope 
of  the  other  two  momentum-transfer  cross-sections  disagrees  less  strongly.  At 
about  0  47  km/sec,  Eq.  (86b)  is  certainly  consistent  with  the  value  inferred  from 
the  Ref.  37  measurement  of  rate  coefficient;  Eq.  (86c)  disagrees  only  mildly  with 
the  extrapolation  of  the  Ref  4  data 
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5.2.4  Comments  on  Literature  Rate  Coefficients. 


The  following  comments  are  peripheral  to  the  main  thrust  of  this  section, 
but  are  relevant  to  the  aspect  of  rate  coefficient  for  Reaction  (80).  Gilmore's 
suggested  rate  coefficient  [Ref.  24],  7xl0‘10  cm3/(molecule  sec),  for  Reaction 
(80)  is  the  average  of  the  value  determined  (by  extrapolation)  in  Ref.  4  and  the  Ref. 
37  value.  Since  an  energy-dependent  cross-section  was  assumed  in  the  Ref.  4 
extrapolation,  we  suggest  that  a  temperature  dependence  should  appear  in  such  an 
averaged  rate  coefficient.  To  express  the  Ref.  4  results  in  a  temperature- 
dependent  rate-coefficient  form,  we  first  use  Eq.  (81b)  to  write 


2,56x10-'2 

‘7°<‘7)  =  p - jw 

l^cm/sec  J 


(87a) 


so  that  we  have  the  product 


cr0(g)  g{T)  =  2.56x  1 0"'2  [  g(T)}'/2 


(87b) 


or,  with  Eq.  (79a), 


er0(g)  <?(T)  =  256x10" 


'8 *rV 

n  /i 


=  5.58x10-3^  . 


(87c) 


(87d) 


Substitution  of  Eq.  (87d),  and  the  value  of  from  Table  25,  in  Eq.  (79b)  gives 

■  j.  -|l/4  3 

XWO^UO-i-S^xlO-  [355]  mole^c  (87e) 

Thus,  the  temperature-dependent  factor  is  [77  300fT  Hence,  we  suggest 
modifying  the  Ref.  24  rate  coefficient  to  become 


X'[U*(02,0)U0+]  =  5.44x10" 


’  T  T*  3 

X'[U*(p2fO)UO*  ]  =  7 x,  o-10  ^ 


The  essence  of  the  above  comment  on  temperature  dependence  also  applies 
to  the  rate  coefficient  for  Reaction  (80)  given  in  Table  I  of  Ref.  5,  where  the 
extrapolation-based  value  from  Ref.  4  is  given,  i.e.,  5.6  x  1 0‘ 10  We  suggest  that, 
in  generalizing  that  value,  one  should  write 


*[im,  0)ixr]  =  7xl0-'0 


molecule  sec 
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5.3  uo++o2->uo;+o. 

5.3. 1  Results  of  Armentrout  and  Beauchamp  [Ref.  4). 


The  cross  section  measured  by  Armentrout  and  Beauchamp  [Ref.  4]  for  the 
reaction 

U0*+02-»U0£+0  (90) 

is  well  represented  by  the  expression 


dU0*(02,0)U0+2]  = 


40x1  O'16 


06 

L^cmJ 


,a55 


7.4x10-'6 


8 


-CM 


lL5 


0.42  *£ms  6.07  eV 


6.07  s  ^  s  22  eV 


(91a) 


The  exponents  of  the  energy  dependence  in  Eq.  (91a)  are  given  in  Ref.  4,  but  the 
coefficients  are  inferred  by  the  writer  from  reading  the  small  graphs  in  Ref  4.  To 
express  the  lower-energy  portion  of  the  cross  section  in  terms  of  the  relative 
velocity,  we  first  write 

^--S^-l47x'0-V 

-0.147^  (92b) 

so  that  substitution  of  Eq.  (92b)  in  Eq.  (9 1  a  1 )  gives 


WU0+(02,0)U0+2]  = 


8.67x10 

Ck 


-15 


cm2,  1.75  6.42  (91b) 


The  cross  section  in  Eq.  (91b)  is  shown  in  Fig.  1 1  as  the  solid  line  and  extrapolated 
to  lower  velocities  as  the  short-dash  line. 


To  develop  the  rate  coefficient  at  temperature  T  corresponding  to  the 
(downward  extrapolated,  where  necessary)  cross  section  in  Eq.  (91b),  we  first 
rewrite  Eq.  (91b)  in  the  form  of  Eq.  (77),  i.e., 


so  that 


274x10' 


0  ' 

Then,  from  Eq.  (79b),  we  have 


2.74x1  O'9  [  g{T)  1 

(93a) 

l£(niu  Um/McJ 

n=\.\ . 

(93b, c,d) 

274x1  O'9 
'  [ giTT  *na" 

(93e) 
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Figure  1 1  Comparison  of  momentum-transfer  cross-sections,  based  on  various 
interaction  potentials,  for  (UO++02)-col1isions  with  the  cross 
sections  for  the  atom-exchange  reaction  UOv+o2  ~»uo;+0. 


rtunpa  ,0)U02  ]  =  2.74x1  O'9  [^(D]-0'0 

■  1-0105 

1.01 


=  274x1  O'9 
=  9.44x1  O'10 


8  kT 


j 
T 


1-005 


300 


cm 


moleculesec 


(93f) 

(93g) 

(93h) 


5.3.2  Results  of  Fite  and  Lo  [Ref.  19]. 

Fite  and  Lo  [Ref.  19]  measured  the  rate  coefficient  for  Reaction  (90)  at 
thermal  energies  and  reported 

*(ur  (02 ,0)U02+  ]  =  (203  ±  0.44)  x  1 0'9  — ^ -  (94) 

2  2  moleculesec 

The  term  "thermal  energies"  is  not  further  defined.  This  rate  coefficient  is 
typically  quoted  and  used  as  being  implicitly  temperature  independent  le.g.,  Ref.  3, 
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p.  48,  Ref.  5,  Table  !,  Ref.  24],  Such  usage  implies  an  assumed  ^''-dependence  for 

the  cross  section,  which  is  certainly  close  to  that  found  by  Ref.  4  at  higher 
energies. 

Similar  to  our  treatment  of  the  Ref,  37  rate  coefficient  in  Eq.  (84),  we  infer 
the  constant  a0  [cf.  Eq.  (77)]  implied  by  the  rate  coefficient  in  Eq.  (94)  if  we 
assume  the  exponent  n  is  the  same  as  in  the  measurements  reported  by  Ref  4, 
namely,  n  *  1.1  (although  use  of  n  =  1.0  would  not  make  a  significant  difference), 
then  from  Eq.  (79b),  for  T=  300  K,  we  have 


(203±0.44)xl0~9 
4.73x104x1.01  1 


=  (4.24±0.92)x10'Hcm2 . 


(95) 


in  Fig.  1 1  this  cross  section  is  plotted  at  a  relative  velocity  of  0.473  km/sec  and 
labeled  as  Ref.  19.  The  central  point  of  this  inferred  cross-section  is  larger  than 
the  extrapolated  cross-section  of  Ref.  4  by  a  factor  of  *2.2. 

Before  leaving  the  rate  coefficient  reported  in  Ref.  19,  it  is  worth  noting 
that  Eq.  (10)  in  Ref.  19  used  to  evaluate  k  for  Reaction  (90)  is  proportional  to  the 
ratio  of  two  rate  coefficients,  each  of  which  is  corrected  upward  by  the  factor 
2.39  reported  by  Halle,  Lo,  and  Fite  (Ref.  27].  Thus,  the  rate  coefficient  given  by 
Eq.  (94)  is  unaffected  by  the  correction  factor  from  Ref.  27. 


5.3.3  Momentum-Transfer  Cross-Sections  for  UCT+02. 

From  Eq.  (40)  and  Table  9  in  Section  2  (or  equivalently,  Eqs.  (37c,  (38b),  and 
(39c)),  we  have  three  momentum-transfer  cross-sections  for  the  collision  pair 
U0%02, 

<r(UO-+0,)|Klw«6.l4xlO-,5<^  cm2  (96a) 

*  1.56x1 0'H  ^  cm 2  (96b) 

^UOVO^^^SxlO-'5^  cm2.  (960 

These  cross  sections  are  plotted  and  respectively  labeled  in  Fig,  1 1  as  "Pol  pot  ", 
"*(r-6)-pot",  and  "Eff.  (r6)-pot.”  The  slope  of  the  polarization-potential  based 
cross-section  closely  agrees  with  that  measured  by  Ref.  4;  the  slope  of  the  other 
two  momentum-transfer  cross-sections  disagrees  more  strongly.  At  about  0  47 
km/sec,  Eq.  (96b)  lies  between  the  value  inferred  from  the  Ref.  19  measurement  of 
rate  coefficient  and  the  extrapolation  of  the  Ref.  4  data. 
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5.4  U+02-»U0;+e. 

5.4. 1  Results  of  Fite,  Lo,  and  Irving  [Ref.  20]  and  of  Halle,  Lo,  and  Fite  [Ref  27] 

The  effective  cross-section  for  the  associative-ionization  reaction 
U+02^U02+e  (97) 

measured  by  Fite,  Lo,  and  Irving  [Ref.  203  is 

a{ U(02,e)U02  r=(!68±0.27)x10-17  cm2  (98a) 

at  a  relative  velocity,  e  g.,  of  6.30  xlO4  cm/sec  [Ref.  20],  Fite  et  ai.  [Ref  20]  also 
explain  that  the  "true"  cross-section  -  instead  of  the  "effective"  cross-section 
given  above  -  is  smaller  by  the  factor  1/1.42  under  the  experimentally-based 
assumption  that  the  cross  section  is  independent  of  velocity  [but  the  velocity 
range  is  not  given].  Thus, 

a\  U(02  ,e)UO+2  T*  =  (t  42)* '  x  a[  U02 ,  e)UO+2  ]e" 

=  (1.18±0.19)xl0*17  cm2.  (98b) 


This  correction  factor  apparently  takes  account  of  the  finite  thermal  velocities  of 
the  02  targets,  as  explained  in  detail  in  Ref.  20,  Appendix  A. 


These  cross  sections  were  revised  upward  by  Halle,  Lo,  and  Fite  [Ref  27] 
who  state  the  revised  effective  cross-section  is 


dtfQ2,e)U0;KL<- (4.0 1±0.55)x  1 O*'7  cm2 ,  (99a) 

larger  by  a  factor  of  2.39  than  the  value  determined  in  Ref  20  Correspondingly, 
we  have 


dU(02,e)UO;]^  =  (1.42)*'  xdU(02,e)U0$  ^ 
<7[UQ2,e)UO2]^9e<3  =  (2.83±0.39)yl0~17  cm2 . 


(99b) 


The  rate  coefficient  corresponding  to  Eg.  (99b),  from  Eg.  (79b)  and  for  n  =  0, 


*[U(02,e)U0;  =  a[  U(02,e)U02  g(T)  <?„ 

'8 


17 


=  2.82x10 
=  1.34x1 0"'2 


J 

z/2 

300 


cm 


moleculesec 


(100) 


This  rate  coefficient  differs  from  that  in  Table  1  oi  Ref  5  in  that  we  have  (a)  used 
instead  of  cr^,^  and  (b)  explicitly  written  the  temperature  dependence. 
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5.42  Results  of  Young,  Dehmer,  Cohen,  Pobo,  and  Wexler  [Ref.  68] 

Reaction  (97)  was  measured  by  Young  et  al.  [Ref.  68]  who  report  results 
which  can  be  expressed  as 


WU02,e)U0+2]  =  1.58x10 


-17 


0.63 


&  i 


.rrvsfc  j 


cm2,  1.2*#^  *6.5 


(101) 


Young  et  al.  [Ref.  68]  state  that 

"...  extrapolation  of  the  deconvoluted  U02  cross  section  to 
vr '  6.3  xlO4  cm/sec  (Erel  =  0.058  eV)  with  the  use  of  the  functional 
form  derived  in  Sec.  Ill  [our  Eq.  ( 101 )]  leads  to  a  value  of  1.58  xlO-17 
cm2,  which  is  to  be  compared  with  the  "Q-true"  of  1.19  x  1 0~ 1 7  cm2 
reported  by  the  authors  of  Ref.  20  at  the  same  average  relative 
velocity  in  their  Maxwellian  beam-gas  experiment." 

"While  this  level  of  agreement  is  encouraging,  it  is  to  be  emphasized 
that  the  conditions  of  the  present  experiment  were  not  optimized  for 
the  measurement  of  absolute  cross  sections, ...  ." 


Figure  i2.  Cross  sections  for  the  associative  ionization  reaction  U+02  -*U02+e 
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in  Fig.  12  we  show  the  results  from  Young  et  al.  [Ref.  68],  from  Fite  et  ai 
[Ref.  20],  and  from  Halle  et  al.  [Ref.  27].  It  is  unfortunate  that  the  work  of  Young  et 
al.  [Ref.  68]  was  done  before  Halle  et  al.  [Ref.  27]  reported  the  upward  correction 
factor  to  their  earlier  work  in  Ref.  20  to  which  Ref.  68  refers.  How  would  Young  et 
al.  have  regarded  the  comparison  of  their  results  with  that  in  Ref  27?  To  this 
writer,  the  level  of  agreement  is  certainly  far  less  encouraging  than  earlier 

5.4.3  Momentum-Transfer  Cross-Sections  for  U+02. 

For  completeness  we  present  momentum-transfer  cross-sections  for  the 
collision  pair  U+02,  but  we  do  not  expect  such  cross  sections  should  be  comparable 
to  the  cross-section  for  Reaction  (97).  From  Eq.  (40)  and  Table  9  in  Section  2  (or 
equivalently,  Eqs.  (27a3)  and  (27c2)),  we  have  two  momentum-transfer  cross- 
sections  for  the  collision  pair  U+02, 

ff[U*0J]|„r<,.|>M  .  124x10 cm2  { 1 02a) 

cr(U+02)[E)t  (r*,_pot  =  7.38x1 0‘IS  ^  cm1.  (1020) 

These  two  momentum-transfer  cross-sections  are  several  orders-of-magmtude 
larger  than  the  cross  section  for  Reaction  (97).  If  we  consider  a  value  of  0.63 
km/sec  for  the  relative  velocity,  then  the  average  reaction  cross-section  - 
derived  from  ( 1 )  the  extrapolated  value  ( 1 .58  x  1 0~ 1 7  cm2)  from  Ref.  68  and  (2)  the 
value  (true,  revised:  (1.68/1.42)  x  10" 17  x  2.39  =  2.83  xl0"i7  cm2)  from  Ref  27  - 
is  0.5(  1.58  +  2.83)  x  10" 17  =  2.21  x  1 0' 1  /  cm2,  smaller  than  the  two  momentum- 
transfer  cross-sections  by  respective  factors  of  763  and  454. 

5.5  U+02-*U0+0. 

5.5. 1  Results  of  Fite,  Lo,  and  Irving  [Ref.  20]  and  of  Halle,  Lo,  and  Fite  [Ref  27], 

It  appears  that  the  reaction  of  U+02,  at  low  energy,  has  two  principal 
branches,  one  described  by  Reaction  (97)  for  associative  ionization  and  the  other 
with  neutral  products  -  with  a  much  smaller  cross-section  -  by 

U+02-*U0+0  (103) 

The  total  effective  cross-section  for  the  U+02  reaction  is  [Ref.  20] 

<r[U02,X)YjSa(  =  (1.78 ±0.45) xiO"15  cm2 .  (104) 

By  subtracting  the  associative  ionization  cross-section  from  the  total  cross- 
section,  we  get  that  for  the  neutral  products, 


dU02,0)U0r  =^[U02,X)Y£,lal  -  a{  U(02,e)U0+2Cfvlsed 


=  (1.78±0.45)x1 0"15  - (4.0 l±0.55)x) O’17 
=  (1.74±0.45)xl0'15  cm2 . 


( 105a) 
(105b) 
(105c) 
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For  the  "true"  cross-section,  we  have 


(106a) 

(1060) 


<7[U02,0)U0]tnj*  =  (142)'lx  g[U02,0)U0r 
<r[U(02,0)U0]tnj*  =  (1.22±0.32)xl0"15  cm2. 


n-0, 


The  rate  coefficient  corresponding  to  Eq.  (106b),  from  Eq.  (79b)  3hd  for 
is 


*[U(02,0)U0r*  *  «r[U(02,0)U0P*  g{T)  60 


This  rate  coefficient  differs  from  that  in  Table  I  of  Ref.  5  principally  in  that  we 
have  (a)  used  <rtn"  instead  of  o'"  and  (b)  explicitly  written  the  temperature 
dependence.  Augmentation  of  k by  the  factor  1.42  to  give  ke1'  results  in 
k*n  =  8.22x10"",  somewhat  smaller  than  the  9x10""  in  Table  I  of  Ref.  5,  for  an 
unknown  reason. 


5.5.2  Momentum-Transfer  Cross-Sections  for  U+02. 

The  two  momentum-transfer  cross-sections  given  by  Eqs.  (102a)  and  ( 1 02b) 
are  approximately  an  order-of-nagnitude  larger  than  the  cross  section  for 
Reaction  (103).  At  a  relative  velocity  of  0.63  km/sec,  the  reaction  cross-section 
is  smaller  than  the  two  momentum-transfer  cross-sections  by  respective  factors 
of  14  and  8.2. 

5.6  U+0-»U0++e. 

5.6. 1  Results  of  Fite,  Lo,  and  Irving  [Ref.  20]  and  of  Halle,  Lo,  and  Fite  [Ref.  27], 

The  effective  cross-section  measured  by  Fite  et  al.  [Ref.  20,  pp.  1 236,  1 239] 
for  the  associative  ionization  reaction 

U+0-*U0++e  (108) 

is 

dU(0,e)U0+  r  =  (97  ±20)  x  ^  ( ,  oga) 

=  (!.63±0.43)x!0"15  cm2  (109b) 

where 

^-<r[U(O2,e)UO+2r  =  (1.68±0.27)xl0",7cm2.  (110) 
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Similarly,  by  using  the  stated  [Ref.  20,  pp.  1240,  1248]  ratio  of  true  cross- 
sections  Wg02)  and  reducing  the  effective  cross-section  Q±  Dy  the  factor  i  42  to 


obtain  the  true  cross-section  <%,  we  have 

C\  «  d'JCO.eiUO*  l1™  =  7T  (Ilia) 

u2 

=  (l.21±0.30)xi0',scm2  (111b) 

where 

77-  =  003  ±20)  (111c) 

6/2 

6?2=<T(U02,e)UO;r"=(1.l8±0.l9)x10*,7cm2.  (Hid) 

To  obtain  the  revised  effective  cross-section  we  write 

6^’rev  =  <j[U(0,e)U0+  C  =  (97±20)x6?2rev  (112a) 

=  (3.89  ±  0.96)  xIO'15  cm2  (1 12b) 


where  the  revised  effective  cross-section  Q from  Ref.  27  is 


67i,r,v-dU02,e)UO+2fi!,  =  (4.0i±0.55)x10-'7cm2.  (112c) 


Similarly,  to  obtain  the  revised  true  cross-section  we  write 

^.<r|U(0,e)U0*KT  =  f 

#>rev  =  (2.90±0.69)xl0''5  cm2 

where 

=^[U02,e)U0*2]^  =  (282±0.39)x10',7cm2 . 


(113a) 
(1 13b) 

(1 13c) 


Fite  et  al.  [Ref.  20,  Appendix  A,  p.  1247]  presume  that  the  same  velocity 
independence  of  cross  section  that  applies  to  associative  ionization  for  the  U+02 
reaction  also  applies  to  associative  ionization  for  the  U+0  reaction.  However, 
those  authors  do  not  state  the  velocity  range  for  which  this  presumption  may  be 
valid. 


The  rate  coefficient  corresponding  to  Eq.  (1 13b),  from  Eq.  (79b)  and  for 
n  =  0,  is 

/F[U(0,e)U0*]^  =  ^[LX0,e)U0+Cr  g{T)6^ 

ll/2 


=  290x1  O'15 


8  kT 
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(114) 


k[U(0,e)U0*E?-  1.89x10 


VI0 


T 


il/2 


300 


cm 


moleculesec 


This  rate  coefficient  differs  from  that  in  Table  I  of  Ref.  5  in  that  we  have  (a)  used 
a™  instead  of  a £?  and  (b)  explicitly  written  the  temperature  dependence 

5.6.2  Momentum-T ransf er  Cross-Sections  for  U+  0. 

From  Eq.  (40)  and  Table  9  in  Section  2  (or  equivalently,  Eqs  (25b)  and 
(26b)),  we  have  two  momentum-transfer  cross-sections  for  the  collision  pair 
U+0, 

rtV* 979x1 0-ls cm2  HI 5a) 

(ll5b) 

These  two  momentum-transfer  cross-sections  are  about  a  half  order-of- 
magnitude  larger  than  the  cross  section  for  Reaction  (97).  If  we  consider  a 
somewhat  arbitrary  value  of  0.6  km/sec  for  the  relative  velocity,  then  the 
reaction  cross-section  <r[U(0,e)U0+  given  by  Eq.  ( 1 1 3b)  is  smal ler  than  the  two 
momentum-transfer  cross-sections  by  respective  factors  of  4.8  and  3. 1 .  These 
factors  become  unity  at  respective  relative  velocities  of  4.8  and  3. 1  km/sec 

5.7  SUMMARY  OF  CROSS  SECTIONS  AND  RATE  COEFFICIENTS  FOR  SELECTED 
URANIUM  OXIDE  REACTIONS. 

In  Table  26  we  have  collated  suggested  cross-sections  and  rate  coefficients 
for  the  more  important  uranium  oxide  reactions,  for  which  the  following  notes 
obtain. 

Notes  for  Table  26:  Suggested  Cross  Sections  &  Rate  Coefficients 

for  Selected  [(UO^+^^Jl-Reactions 


a.  "True”  cross-section  [Ref.  20],  revised  [Ref.  27],  Eq.  (113b).  Reaction  rate 
coefficient  is  given  in  two-temperature  form  in  Table  26  as  well  as  one- 
temperature  form,  cf.  Eq.  (114).  See  Section  5.6  for  details. 

b.  See  Fig.  13.  Velocity  dependence  [Ref.  4];  coefficient  is  determined  by 
average  value  at  g  =  4.74  x  1 04  cm/sec  of  ( 1 )  extrapolated  value  (1.18  x 
10~14  cm2)  from  Ref.  4(cf.  Eq.  (81b)  and  Fig.  10)  and  (2)  value  (1.84  xiO-1 4 
cm2;  cf.  Eq.  (85d))  inferred  from  Ref.  37  rate  coefficient  by  assuming  g-]/ 2 
velocity  dependence  from  Ref  4.  See  Section  5.2  for  details. 
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Notes  for  Table  26:  Suggested  Cross  Sections  &  Rate  Coefficients 

for  Selected  [(UOn)nH+(0,02)]-Reactions.  (Cont'd) 


c.  See  Fig.  14  Velocity  dependence  [Ref.  4];  coefficient  is  determined  by 
average  value  at  g  =  473  x  1 04  cm/sec  of  ( l )  extrapolated  value  ( 1 .97  x 

10~ 14  cm2)  from  Ref.  4  (cf.  Eq.  (91b)  and  Fig.  1 1 )  and  (2)  value  (4.24  x10_M 
cm2,  cf.  Eq.  (95))  inferred  from  Ref.  19  rate  coefficient  by  assuming  g-1  1 
velocity  dependence  from  Ref.  4.  See  Section  5.3  for  details. 

d.  See  Fig.  15.  Velocity  dependence  [Ref.  68],  coefficient  is  determined  by 
average  value  at  g  =  6.3  x  1 04  cm/sec  of  ( 1 )  extrapolated  value  ( 1  58  x  1 0_  1 7 
cm2)  from  Ref.  68  (cf.  Eq.  (101)  and  Fig.  12)  and  (2)  value  (true,  revised 

( 1 .68/ 1 .42)  x  1 0- 1?  x  2.39  =  2.83  x  1 0' 1  ?  cm2,  cf.  Eq.  (99b))  from  Ref.  27. 
See  Section  5.4  for  details. 

e.  Subtract  a[ U(02,e)U0; [Ref.  27]  from  the  total  cross-section 
dU02,X)Y]££,  [Ref.  20],  See  Section  5.5  for  details. 


Figure  1 3.  Comparison  of  our  suggested  velocity-dependent  cross-section  with 
experimental  data  for  the  reaction  IT  +02  ->U0*  +  0 
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Figure  14.  Comparison  of  our  suggested  velocity-dependent  cross-section  with 
experimental  data  for  the  reaction  U0++02  -»U02  +  0 


Figure  15.  Comparison  of  our  suggested  velocity-dependent  cross-section  with 
experimental  data  for  the  reaction  U+Oj  -» U02  +e. 


Table  26  Suggested  cross  sections  and  rate  coefficients  for  selected  [(UOn)m+  +  (0,02)]-reactions 
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bee  previous  g  *  relative  velocity,  cm/sec  4  =  U-temperature 

pages  for  notes  //  =  reduced  mass  t2  =  0- temperature 
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APPENDIX  A 


RELATION  BETWEEN  DIFFERENT  EXPRESSIONS  FOR  THE 
DIFFUSION  COEFFICIENT  AND  MOBILITY 


A.  1  THE  PROBLEM. 

In  the  (weakly- Ionized)  plasma  physics  literature  one  sees  the  (ordinary) 
ion  diffusion  coefficient  D  and  ion  mobility  X  (sometimes  written  as  which 
here  is  reserved  for  reduced  mass)  written  as 


o-J-L 

mv 

<*.  e  _eD 
mv~  kT 


by,  e  g.,  Refs.  33,  21,  and  2  (Book,  p.  266),  or  as 


D.*±t 

m 

er  eD 
m  kT 


(Ala) 

(A2a) 


(Alb) 

(A2b) 


by,  e.g.,  Refs.  49,  46,  and  64. 

On  the  other  hand,  elsewhere,  for  both  binary  and  multicomponent  gases,  one 
sees  the  diffusion  coefficient  written  differently  than  Eqs.(Ala)  and  (Alb) 

Binary  Gas.  In  a  binary  gas  consisting  of  particles  with  masses  mt  and  mt 
and  with  concentrations  n,  and  nk,  the  diffusion  coefficient  [Ref.  7,  Pt.  B,  p.  39,  Eq. 
(15.27);  equivalent  to  Dalgarno's  Eq.  (5),  p.  644  in  Ref.  8]  is 


3  x 

rQ  kTy 

Z200tk(nl+nk) 

ft,  *) 

where  Q0,k  is  an  average  diffusion  (or  momentum-transfer)  cross-section  given  by 


-oo 

=  x*o0ik{x)e-x d* 


(A4a) 


with 


A- 1 


/i,kg2 

*  =  (A4b) 

fi,k  =  n i  mk  /  (/77,  +mk) ,  (A4c) 

and  g  being  the  relative  speed  of  the  two  types  of  particles. 

Multicomponent  Gas.  Dalgarno  [Ref.  8,  p.  649,  referencing  Hirshfelder  et  al. 
(1954)]  states  that  when  the  concentration  of  one  of  the  components  is  very  small 
compared  with  those  of  the  others,  the  simplified  formula  for  the  diffusion 
coefficient  of  the  rare  component  (D/)  is  given  by 

o;'*1rkD-l,  (A5) 

k 

where  fk  is  the  fractional  concentration  of  the  *th  particles.  (Note:  D,k  must  be 

evaluated  for  nk  equal  to  /w  This  fact  is  more  readily  recognized  from  sources 
other  than  Dalgarno,  such  as  Refs.  18  and  66.) 

We  now  reconcile  these  different  expressions  for  the  diffusion  coefficient 

A. 2  RECONCILIATION. 

First,  we  assume  there  is  no  question  about  relating  the  collision  time  r  to 
the  collision  frequency  v  by 


t  =  \/v  (A6) 

(cf.  Ref.  7,  p.  186,  Eq.  9. 14)).  Next,  we  assume  we  are  interested  in  an  average 
collision  frequency,  vUk,  corresponding  to  an  average  momentum-transfer  cross- 
section,  obtained  by  averaging  over  a  Maxwellian  distribution  of  velocities.  (We 
assume  only  one  temperature  obtains,  for  -‘r  >1  icity  -  the  argument  is  extendable 
to  the  two-temperature  case.)  This  vUk  is  to  the  effective  collision  frequency 
in  the  laboratory,  related  to  the  collision  frequency  v,k  in  the  center-of-mass 
system  [Ref.  7,  p.  190,  Eq.  (9.40)]  by 


Vlik  * 


m, 


!l _ tj 


m,  +mk 


’  ik  > 


(A7) 


where  vik  is  given  [Ref.  7,  p.  1 89  (Eqs.  (9,33),  (9.34),  (9.35),  &  (9.36)),  with  a 
reference  to  Ref.  6]  by 
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with 


V>k  ~  3  ^kdilfioi  k 


Qt  k  ~ 


(A8a) 


(A8b) 


=  4g,k 


2  2 kT 
c  = - 


(A8c) 


(A8d) 


Thus,  by  collating  expressions,  we  find  the  diffusion  coefficient  as  given  by  Eq. 
(Ala)  for  a  binary  gas  becomes 


n  kT  kT 

mv  m,vUk 

kT  mf + mk  1 

w>  rn, \ 

=  kT__ l_ 
fhk  Kk 
kT  3 

T'ik  ^nkQik^Oik 

_3  kT  1 

4  f^fk^k  ^kT  J  Opik 
/  \V2 

3*  1  8X-7~ 

32  nkQ0fk  Jt/ifk 


(A9a) 


(A9b) 


(A9c) 


(A9d) 


(A9e) 


(A9f) 


which  is  the  same  expression  we  have  for  Dik  in  Eq.  (A3)  if  /?,  «  nk. 

Thus,  by  this  correspondence  the  proper  expression  for  the  r  in  Eqs.  (Alb) 
and  (A2b)  for  a  binary  gas  is 


1  1  m.  n 

"  K*  ?u‘  kTD’> 

3  m,  *fi,k 

4  Mlk^k^Olk  %kT 


(AlOa) 


(A  !0b) 


A. 3  PARAMETERS  SPECIALIZED  FOR  THREE  POTENTIALS. 

We  now  want  to  specialize  the  expressions  for  Q0,  r,  X,  and  D  for  the 
polarization,  (r&),  and  (hybrid)  effective  (r&)  potentials 

A.3. 1  Q0,k ,  The  Mean  Momentum-Transfer  Cross-Section. 


A.3.1.1  Polarization  Potential.  For  a  polarization-potential  interaction,  the 
velocity-dependent  momentum-transfer  cross-section  for  singly-charged  ions  is 
[Ref.  7,  Pt.  A,  p.  218] 

a otk  =  2.2  //I't&f*  (Alla) 

=  259x1 0'9[<3-; / /liA'JgftY* ,  ( A 1  lb) 


where  a*  is  the  polarizability  of  the  neutral  target  particle,  ak  is  in  units  of 
10“24  cm3,  fiik  is  the  reduced  mass  of  the  ion  and  neutral,  (i(Ajk)  is  in  atomic  mass 
units,  e  is  the  electronic  charge,  and  gik  is  the  relative  speed  of  the  ion  and 
neutral.  By  using  Eq.  (A4b),  we  find  that  Eq.  (Alla)  becomes 

aD,k  =  2  2\n(ake* /2kTln /4x,  (A12) 


so  that  Eq.  (A4a)  becomes 


Qoik  - 


2.2  \x 


a, 


\l/2 


\7kT) 


rd5/  2) 


where  [Ref.  1 ,  p.  255] 

/*(/?)  =  f  x"~'e~*dx 

Jo 

7'(5/2)  =  [(1  x3)/22]  nv  2)  =  1.3293. 

Thus, 

„  i  22U/,(5/2)f^r 

@Oik IpoLpot. =  2  [ 2kT  ^  ■ 

«l33x1<r,5(<r;/nw2. 


(Ai3a) 

(A  1 4a) 
(A  1  -4b) 

(A) 3b) 
(A13C) 


A.3.1.2  (r6)-potential.  The  asymptotic  low-energy  momentum-transfer 
cross-section  for  the  r~6  potential, 

*r)=>-2e0{rm/r)\  (A  15) 

obtained  by  combining  Eqs.  (5),  ( 1  c),  and  ( 1 2)  in  Section  2,  is 
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aig,k)  =  xD2  ^{jc) 


=C 


=  C 


irle^eV) 


\l/3 


M,k  9%  , 
(dcMS) 


C  =  0.4342x(3xl.60xl0‘l2)l/3x4;r 
=  9.204  xlO"4 
r  =  7.773x1  OV 

By  using  Eq.  (A4b),  we  find  Eq.  (A  16b)  becomes 

wJ 

so  that  Eq.  (A4a)  becomes 

cm/Z) 


C 

a0>k  =  -yJl 


Ooik 


<r'e  ** 
mSs. 


2kT) 


where  (Ref.  1,  p  255] 

r(Q  /  3)  =  [(2 x  5)  / (32)]  f(2 / 3)  =  1. 5046 . 

Thus 


-*W*|(r-*)-pot 


» 6.92x1  O'4* 

[2kT 


(A  16a) 
(A  1 6b) 

(Aloe) 

(A16d) 
(A  1 6e) 
(A16f) 

(A  1 6g) 

(A  1 7a) 

(A18) 

(A  1 7b) 


A.3.1.3  Effective  (r~6)-Potential.  As  discussed  in  detail  with  respect  to 
(Ba^+N^-collisions  in  Appendix  C,  we  want  to  introduce  an  effective  (r~&)- 
potential  determined  by  requiring  the  equalities  of  mobilities  (at  STP  conditions) 
based  on  (a)  an  effective  (r_6)-potential  and  on  (b)  the  polarization  potential.  This 
requirement  is  equivalent  to  an  equality  of  mean  momentum-transfer  cross- 
sections  (at  STP  conditions)  based  on  the  same  two  potentials,  since  the  mobility 
is  (inversely)  proportional  to  the  mean  momentum-transfer  cross-section.  To 
determine  this  effective  (r_6)-potential,  we  equate  the  right-hand  members  of 
Eqs.  (A  1 3c)  and  (A  1 7b),  set  T  to  273. 1 5  K  (■  Tm),  and  solve  for  the  parameter 

(^f<1)w|E„,r.,<ol  =  4.94xlO-"i(a’;)w.  (AI9) 

Note  that  determining  the  parameter  {r^ejn  in  this  manner  results  in  its  having 
no  dependence  on  any  property  of  the  incident  ion  /  However,  the  velocity- 
dependent  momentum-transfer  cross-section  still  has  a  reduced  mass  factor  Use 
of  Eq.  (A19)  in  (A!6b,c)  gives 
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(r*)-pot  “LEff.  ^  ^)I/3 

~  to;)1'2 

E”  U(4.)^r 

(A20a1 ) 

(A20a2) 

with 

<fE„  =4.55x1  O’19 

(A20a3) 

a,,.  =  3. 84x1 0‘"  . 

(A20a4) 

Use  of  Eq.  (A19)  In  (A17b)  gives 

0«*Un**  “5  25x1  or'^i^or0 . 

(A20b) 

A.3.1.4  Polarization  Potential  and  Effective  (r6)-Potential.  From  Eqs  (A20) 
and  (A  13c)  we  have 

/  T  v* 


/stp  =  273.15K. 


A 


Z?/*lpol  pot 


(A2!a) 

(A21b) 


Hence,  from  Eqs.  ( A2 1 )  and  (A  1 3c),  we  have 


@Oik  jetf  (r-*)-pot  =1-33x10  f  j 

V'STP 


(A22) 


We  can  combine  Eqs.  (A  13c)  and  (A22)  by  writing 


■l33x,cr'J(^fWcmJ 

(A23a) 

0, 

polarization  potential 

(A23b) 

1/6, 

effective  (z"-6)— potential . 

&  = 


A.3.2  r,  The  Collision  Time. 

Substituting  Eqs.  (A17b)  and  (A23)  into  (A  10b)  gives 
I  1.90x1 0*7  "  ** 


4 


'  r 

im 


*  T 


i  3.88x10®  4  f,STP 


\P 


(A24) 

(A25a) 
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(A25b) 


P- 


'0,  polarization  potential 
1/6,  effective  (/^-potential. 


A.3.3  X,  The  Mobility. 


A.3.3. 1  Polarization  Potential.  From  Eqs.  ( A2),  (A9f ),  and  (A  1 3b),  we  have 


St  I  -Ul  _#*stp| 

Ipolpot.  -  ^  Ipolpc 

<3f»STP|  13.8 


*  Ipolpot. 


cm/ sec 
volt/cm 


nL  =  2.69x1  O'9  cm-3 


(A26a) 

'A26b) 

(A26c) 


Equation  (A26b)  is  equivalent  to  that  given  by  Dalgarno  [Ref.  8,  p.  652  or  Ref.  16], 


J£STP  = 


35.9  cm/ sec 


(A27) 


[/i(A)a]V2  vo It/ cm 

where  a  is  in  atomic  units  (a03)  instead  of  1 0~24  cm3  as  for  our  a *,  the  reduced 


mass  fi(A)  is  in  units  of  the  proton  mass,  and  Xsrp  is  referred  to  a  constant  gas 
density  of  2.69  xlO19  cm-3. 


A.3.3.2  (r6)-potentlal.  From  Eqs.  (A2),  (A9f),  and  (A 1 7b),  we  have 


(r-*)-poL 


tstpi  6.82x1  O'15  cm/sec 

**  "[MJPlGcJ*  volt/cm 


(A23a) 

(A28b) 


A.3.3. 3  Polarization  Potential  and  Effective  (r6)-Potential 
(A9f),  and  (A23)  we  have 


qf-STP  Of-STP  g*-STP 

•*1*  *  pot  pot  "  •nT*  Eff >-*)-pot. 


From  Eqs.  (A2), 


(A29a) 

(A29b) 


P  = 


0,  polarization  potential 
1/6,  effective  (r~*) -potential. 


(A29C) 
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A.3.3.4  Collation  of  Ba*  Theoretical  Mobilities  at  STP  Conditions.  The 
parameters  for  Ba+  theoretical  mobility  (and  diffusion  coefficient,  presented 
below)  are  collated  in  Table  27.  For  values  of  a',  see,  e.g„  Ref.  7,  pt.  a,  p  219 
These  parameters  are  used  in  Eq.  (A26b)  to  obtain  the  Ba+  theoretical  mobilities  in 
Table  28.  As  discussed  in  Appendix  C,  the  theoretical  value  for  the  Ba+  mobility  in 
N2  is  larger  than  the  experimental  value  by  only  about  four  percent. 


Table  27.  Parameters  for  Ba+  theoretical  mobility  and 
diffusion  coefficient  in  atmospheric  species. 


k 

at‘k 

.  /  =  Ba+ 

[/t(4k)ar^r 

0 

0.79 

1433 

3.36 

N2 

1.76 

23.26 

6.40 

02 

1.59 

25.95 

6.42 

Table  28.  Ba+  mobilities  in  atmospheric  species  at  STP  conditions 
(Based  on  Eq.  (A26b)  and  properties  in  Table  27.) 


Species  k 
0 
n2 
02 


STP 

X  (Ba+,/0,  [(cm/sec)/(volt/cm)] 

(either  polarization  or  effective  (r'6)-potential) 

410 

2.16 

2.15 


A.3.4  D ,  The  Diffusion  Coefficient. 


From  Eqs.  (AS)  and  (A2)  we  have 


o;' = lr„  vl 

k 


(A30a) 

(A30b) 
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A.3.4. 1  (r'6)-Potential.  From  Eqs,  (A28)  and  (A30b)  we  have 

T  cm 


1  1 5.8  j 

( T 

'STP 

"tot 

l  T  ) 

\l/6 


—  (A3 1 ) 


A.3.4.2  Polarization  Potential  and  Effective  (r6)-Potentla).  From  Eqs.  (A29) 
and  (A26b)  we  have 


n.  3.21x1  V{T^\ 
u>\>- — ~n - T 

" tot  \  ' 


T 


cm 


sec 


P- 


0,  polarization  potential 
1/6,  effective  (/~6)- potential. 


(A32a) 


(A32b) 


A3. 43  Working  Formula  for  Ba*  Diffusion  Coefficient.  The  parameters 
necessary  to  apply  Eq.  (A32)  to  Ba+  are  collated  in  Table  27  and  applied  to  Eq. 
(A32)  to  obtain  Eq.  (A33).  Results  of  evaluating  Eq.  (A33)  for  a  range  of  altitudes 
are  given  in  Table  34 


/^♦Ip  - 


3.2  lx  1 016 


n, 


tot 


5ie 
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T 


3.36£ +6. 40/^+642/; 


(A33) 
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APPENDIX  B 


ION-NEUTRAL  COLLISIONS  IN  MULTICOMPONENT  GAS 

B.  1  COMMENTARY  ON  SELECTED  LITERATURE  FORMULAS  FOR  COLLISION  TIME 

Linson  &  Workman  [Ref.  49,  Eq.  (2a)]  state  that 

"...  the  ion-neutral  collision  time,  r,  is 

r= 4.5x1 09//y,  ,  sec  (Bi) 

where  Nn  is  the  atmospheric  neutral  density  in  cm-3  ... .  A  typical 
ionospheric  temperature  of  1000  K  has  been  taken  for  the  barium  ions. 
The  value  of  r  was  obtained  from  the  temperature-corrected 
measurement  of  the  mobility 

fi  =  erlM  (B2) 

of  barium  ions  in  nitrogen  by  Mitchell  and  Ridler  [Ref.  58]  and 
reprinted  by  Brown  [Ref.  12]." 

We  believe  there  are  two  shortcomings  in  Eq.  (Bl ),  due  to  improperly  (a) 
correcting  for  the  ionospheric  temperature  of  about  1000  K  and  (b)  treating 
mobilities  in  a  multicomponent  gas.  (This  latter  statement  is  equivalent  to 
objecting  to  treating  O2  (which  doesn't  matter)  and  0  (which  does  matter)  as  if 
they  were  N2.) 

The  LW-70  authors  apparently  became  aware  of  the  first  shortcoming 
because  a  later  paper  by  Linson  and  Baron  [Ref.  45]  states  that  the  ion-neutral 
collision  time  is 

r,  =  9x109//V„  .  (B3) 

An  improved  but  still  slightly  incorrect  formula,  related  to  Eq.  (B3),  is 
presented  by  Linson  &  Baxter  (Ref.  46,  p.  58]  and  reproduced  here: 

"An  expression  for  the  collision  time,  r  =  \/v,  can  be  obtained  from 
measurements  of  the  mobility  of  barium  ions  in  nitrogen  gas  [Ref.  58] 
resulting  in 

r  =  8.6x1 0,s//7„  sec  (4.54)  (B3al) 

where  h„  is  an  effective  neutral  concentration  in  particles/m3, 

hn  =  /7Nt  +  /7Oj+0.8/7o.  (4.55)  (B3a2) 

In  deriving  Eq.  (4.54)  it  has  been  assumed  that  the  collision  cross 
section  varies  inversely  as  the  square  root  of  the  temperature  and 
that  the  collision  cross  section  for  oxygen  molecules  and  atoms  is  the 
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same  as  it  is  for  nitrogen  molecules.  The  fact  that  the  lighter  oxygen 
atoms  are  less  efficient  in  stopping  a  heavier  barium  ion  is  expressed 
by  the  appropriate  Langevin  factor, 

where  /\  and  /%  are  the  masses  of  the  nitrogen  molecule  and  oxygen 
atom  respectively." 

The  proper  way  to  obtain  the  collision  time  in  a  multicomponent  neutral  gas 
mixture  can  be  done  in  at  least  a  couple  of  (equivalent)  ways.  In  Section  B.2  we 
compute  it  in  terms  of  mobilities,  here,  we  show  another  quite  direct  way,  by 
summing  the  collision  frequencies  for  the  ion  in  the  various  neutral  gases. 


Summation  of  the  collision  frequencies,  under  the  assumption  of  the 
polarization-potential  interaction,  has  been  used  frequently,  including  a  recent 
paper  by  Lathuillere,  Wickwar,  and  Kofman  [Ref.  41],  from  which  we  reproduce  the 
relevant  portion.  (In  the  Ref.  41  Eq.  (4)  we  have  corrected  the  typographical  error 
in  an  exponent,  changing  it  from  IO-1 1  to  IO"9.) 


"Neutral  density  and  collision  frequency  are  connected  by 

\t/2 


v,„  =  259x1 0*9/?, 


an 


(1) 


{Banks and Kockarts  (Ref.  7]],  where  vin  is  the  collision  frequency 
for  momentum  transfer  between  ion  species  /  and  neutral  species  n, 
measured  in  the  center  of  mass  frame;  n„  is  the  number  density  of 
neutral  species  n,  /e(/,n)  is  the  ion-neutral  reduced  mass  (in  atomic 
mass  units),  and  ^  is  the  neutral  gas  atomic  polarizability  (in  units 
of  10-24  cm'5). 


"The  ion-neutral  collision  frequency,  vin,  used  in  the  calculation  of 
the  incoherent  scatter  spectrum  Is  the  momentum  transfer  collision 
frequency  measured  in  the  laboratory  frame: 


V,n  = 


m. 


m,  +mn 


(2) 


where  m,  and  m„  are  the  ion  and  neutral  particle  masses  (in  atomic 
mass  units),  respectively.  Combining  ( 1 )  and  (2)  and  summing  over 
the  three  principal  neutral  species  in  the  E  region,  we  obtain 

v,n-KN  (3) 

for  each  ion  species,  where  N  in  the  total  neutral  density, 

K  =  2.59x1  O'9  (1//77/X1.761/2/a(/,N2)'/2%(N2) 

U59w/rt/,02)w%(02)+0.79v2/f(/,0)l/2%(0)]  (4) 
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and  %(n)  represents  the  percentage  of  neutral  species  n  in  the  neutral 
gas." 

To  obtain  explicit  numbers,  we  write 


Tr  =  l/Kr 

-  2  59t TF’  [V'-76A.«.  [N!J+  f02]+  ;o.79a.o  (Olf ' 

-  2  [6. 40 [Nj]  +  6.42 {02]  +  3. 36 (OIF' 

8.29x1 09  ,  .  4  ..  t  ,  ,„,u. 

=  fN2M02],0.525[0]fQr  P0lar,2at10n  P°tential  (B3b > 

which  we  believe  to  be  the  proper  version  of  Eqs.  (4.54)  and  (4.55)  in  Ref.  46  (our 
Eqs.  (B3al)  and  (B3a2)). 

B.2  FORMULAS  FOR  MOBILITY,  COLLISION  TIME,  AND  DIFFUSION  COEFFICIENT  IN  A 
MULTICOMPONENT  GAS. 

Dalgarno  [Ref.  8,  p.  659]  notes,  from  the  theory  of  Chapman  and  Cowling  [Ref 
14],  that  the  mobility  X  (a  symbol  used  here  in  preference  to  yz  as  in  Eq.  (B2), 
since  we  reserve  /z  for  reduced  mass)  of  an  ion  in  a  multicomponent  gas  is  given 
by 

X-' =.  £4  X;'  (B4) 

k 

(known  as  Blanc's  law),  where  fk  is  the  fraction  of  the  /Uh  particles  (f*=nk/nt0t9/) 
and  Xy  is  the  mobility  of  the  Ion  in  the  /Uh  gas.  In  using  Eq.  (B4)  it  is  important 
that  Xy  be  evaluated  for  the  total  number  density.  (See  the  related  comment  with 
respect  to  Eq.  (A5).) 

Thus,  to  compute  an  effective  ion-neutral  collision  time  ( r)  for  an  ion  in  a 
multicomponent  gas,  one  should  use 


r'^-x-'^-YryX;'. 
m  mj  k  * 


(B5) 


We  now  specialize  Eq.  (B5)  for  two  different  interaction  potentials, 
polarization  and  effective  r-6,  By  using  Eq.  (A29)  we  get 
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m  nL 
o  ntol 


P- 


0,  polarization  potential 
1/6,  effective  (/"^-potential. 


(B6a) 


(B6D) 


For  Ba+,  Eq.  (B6)  becomes 
3.83x1 09 


ria  = 


n, 


'tol 


^STP 

T 


where  we  have  used  values  for  from  Table  27. 


(/r0/4.10)+(4j/2.16)+(4/2.15) 

STP  / 


(B7) 


The  effective  ion-neutral  diffusion  coefficient  is  also  properly  defined  by 
[Ref.  7,  Pt.  B,  p.  i  64] 


m,  LK* 

k 


(B8) 


B  3  COMPARISONS  OF  Ba+  COLLISION  TIME. 

The  CIRA-65  Model-6  1800-hr(Ref.  15]  atmosphere  is  of  special  interest  in 
Appendix  D  for  comparisons  being  made  there.  Selected  properties  of  that 
atmosphere  for  the  140-  to  300-km  altitude  range  are  given  in  Table  31  and  have 
been  used  in  evaluating  the  Ba+  collision  time  according  to  several  formulas 
presented  here,  with  results  in  Table  29.  At  200-km  altitude,  e.g.,  our  procedure, 
per  Eq.  (B7),  gives  1 . 1 8  sec  and  0.935  sec  respectively  for  the  polarization 
potential  and  effective  (r~6)-potential,  whereas  treating  the  0  and  O2  as  N2  in  Eq 
(B7)  -  in  the  spirit  of  the  Ref.  45  equation  -  gives  0.954  sec  and  0.752  sec, 
respectively;  our  collision  times  are  longer  by  a  factor  of  1.24.  Again  at  200-km 
altitude,  in  comparison  with  Ref.  45  (our  Eq.  (B3)),  our  values  differ  by  respective 
factors  of  1 .  1 3  and  0.899;  in  comparison  with  Ref.  46  (our  Eq.  (B3a)),  by  factors  of 
1.09  and  0.866. 
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Table  29.  Ba+  collision  times  in  CIRA-65  Model-6  1800-hr  atmosphere 
per  several  formulas.  (See  Table  31  for  atmospheric 
properties.) 


Altitude 

_ Collision  Time  r.  sec _ 

km 

Pol.  Pot. 

RSHI 

mam 

wmsm 

140 

1.03  E-01 

8.87  E-02 

1  00  E-01 

1  00  E-01 

160 

2.96  E-01 

2.44  E-01 

2.78  E-01 

2.82  E-01 

180 

6.38  E-01 

5. 12  E-01 

5.79  E-01 

5.92  E-01 

185 

7.47  E-01 

5.96  E-01 

6.72  E-01 

6.92  E-01 

200 

I.I8E+00 

9.35  E-01 

1.04  E+00 

1.08  E+00 

220 

2.0 1  E+00 

1.57  E+00 

1.70  E+00 

1.80  E+00 

240 

3.23  E+00 

2.50  E+00 

2.64  E+00 

2.82  E+00 

260 

496  E+00 

3.83  E+00 

3.92  E+00 

4.24  E+00 

280 

7.39  E+00 

5.67  E+00 

5.66  E+00 

6.18  E+00 

300 

1.07  E+01 

8.22  E+00 

7.97  E+00 

8.80  E+00 

Eq.  (B7) 


Eq.  (B7) 


Eq.  (B3) 


Eq.  (B3a) 
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Ba  AND  Ba+  COLLISIONS  WITH  ATMOSPHERIC  SPECIES 
C.1  (Ba+0)-C0LLISI0NS:  MOMENTUM-TRANSFER  CROSS-SECTION. 

For  such  collisions  we  are  aware  of  only  one  (simple)  approach  for 
estimating  the  momentum-transfer  cross-section,  that  described  for  atom-atom 
collisions  in  Section  2.1.1.  For  relative  speeds  g  below  a  few  times  1 05  cm/sec, 
the  asymptotic  form  of  the  momentum-transfer  cross-section  based  on  the  r~& 
potential  is  given  by  Eq.  (A16c),  where  rm  is  the  radius  at  which  the  minimum 
energy  (- e0 )  occurs  in  the  Lennard-Jones  (12,6)-potential.  For  (Ba+0)-co11isions 
we  adopt  the  potential  parameters  from  Huber  and  Herzberg  [Ref.  34,  p.  72]: 


r**5.79eV  (Cl  a) 

rm*r,.  1.9397  A  (Clb) 

and  the  reduced  mass  number  from  Table  27.  Then,  from  Eq.  (Al6c)  we  have 

‘rtBa+O^-n**  =2.16x10-' ' grw  cm2  (C2a) 

=  1.00x1 0-'4  cm2.  (C2D) 


Equation  (C2)  is  plotted  as  the  long-dash  line  in  Fig.  16.  (Equation  (C2)  is  very 
close  to  what  we  found  for  (U+0)-collisions,  based  on  the  (r6)-potential,  in 
Section  2.3.2:  2.11  x  10_n  g-2/3.)  It  will  be  explained  later  that  our  final 
suggestion  for  the  momentum-transfer  cross-section  for  (Ba+0)-collisions  is  not 
Eq.  (C2)  but  Eq.  (C4). 

C.2  (Ba^O)-COLLISIONS:  MOMENTUM-TRANSFER  CROSS-SECTION. 

For  such  collisions  we  are  aware  of,  and  consider,  two  (simple)  approaches 
for  estimating  the  momentum-transfer  cross-section.  Unaware  of  a  strong  basis 
for  choosing  between  them,  we  will  offer  a  compromise,  which  will  lead  us  to 
modify  our  formula  for  (Ba+0)-collisions, 

C.2.1  Polarization  Potential. 

For  a  polarization-potential  Interaction,  the  momentum-transfer  cross- 
section  is  given  by  Eq.  (A 1 1b),  which,  with  parameters  from  Table  27,  gives 
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Figure  16.  Velocity-dependent  momentum-transfer  cross-sections  for  (Ba+0)- 
collisions,  based  on  (r6)-  and  effective  (r6)-potentials  and  for 
(Ba++0)-  collisions,  based  on  polarization  and  effective  (r6)- 
potentials. 

^BaVO)^  =6.08x10-'V'  cm2  (C3a) 

=  6.08x1  O’15  cm2.  (C3b) 

Equation  (C3)  is  plotted  as  the  short-dash  line  in  Fig.  16. 

C.2.2  (r_6)-Potential. 

in  Section  2  we  indicated  a  basis  for  reservations  about  the  applicability  of 
the  polarization  potential  for  UCC-collisions,  especially  for  the  higher  velocities 

(s  5  km/sec)  being  considered  by  IANL  personnel  early  in  these  studies.  We  have  a 
similar  reservation  with  respect  to  Ba+  collisions.  The  parameter  y  in  the 

( 12.6.4) -potential  [cf.  Section  2. 1.2.2]  is  a  measure  of  the  relative  strength  of  the 
r6  and  rA  energies.  For  y  =  1,  the  (12,6)-potentiai  obtains;  for  ^*0,  the 

(12.4) -potential  obtains. 

Not  aware  of  potential  parameters  being  available  for  BaO+,  we  f°ntative!y 
assume  they  can  be  approximated  by  those  for  BaO,  given  by  Eqs.  (C 1  a)  and  (C  l  b) 
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in  that  case,  y  would  be  equal  to  0.954,  close  to  the  value  (0  960)  we  found  for 
U0+  [Section  2.2.2.33  Such  a  value  for  y  suggests  the  importance  of  the  r-6  term 
compared  with  the  r-4  term,  if  one  ignored  the  r-4  term  completely,  and  also 
adopted  the  BaO  potential  parameters  for  BaO\  then  Eq.  (C2)  would  be  an 
approximate  expression  for  the  momentum-transfer  cross-section  for  (Ba++0)- 
collisions  based  on  a  "quasi-real istic"  potential.  (Note  that  Eq.  (C2)  is  extremely 
close  to  what  we  found  for  (U++0)-collisions  in  Section  2. 1.2.3,  2. 15  xlCH  1  x 

g_2/3) 

C.2.3  A  Compromise  (?):  An  Effective  (r_6)-Potential. 

Note  that: 

(a)  The  experimental  mobility  of  Ba+  in  N2  [Powell  and  Brata,  unpublished  (Ref 
58],  see  Section  C.3.2]  can  be  predicted  to  within  about  four  percent  by  use 
of  the  polarization  potential  (see  Section  C.3). 

(b)  A  similar  situation  may  hold  for  (Ba++0)-collisions. 

(c)  There  are  reservations  about  the  applicability  of  the  polarization  potential 
at  the  higher  velocities. 

(d)  There  is  concern  about  Eq.  (C2)  giving  a  larger  cross-section  for  (Ba+O)- 
colllslons  than  that  given  by  Eq.  (C3)  for  (Ba++0)-collls1ons  at  (room) 
thermal  energies. 

In  view  of  these  considerations,  we  use  Eq.  (A19)  in  Eq.  (A16c),  with 
parameters  evaluated  from  Table  27,  to  get  our  final  expression  for  the 
momentum-transfer  cross-section  for  [(Ba,Ba+),0]-collisions, 

<r[(Ba,Ba+)+ojEfJ(r^Hx)t  =1.40x10'"  g'2n  cm2  (C4a) 

=  6.50x1  O'15  cm2.  (C4b) 

Equation  (C4)  is  plotted  as  the  dotted  line  in  Fig.  16. 

C.3  (Ba++N2)~COLUSIONS:  MOBILITY. 

C.3. 1  Introduction. 

We  are  unaware  of  any  potential  parameters  for  (Ba++N2)-collisions. 
However,  there  is  a  mobility  measurement  which  we  will  compare  with  a 
theoretical  value  based  on  the  polarization-potential  interaction.  Owing  to 
reservations  about  the  applicability  of  the  polarization-potential  interaction  at 
the  higher  velocities  -  if  we  extrapolate  to  N2  our  considerations  of  (Ba++0)- 
collisions  -  we  will  use  the  theoretical  mobility  based  on  the  polarization 
potential  (essentially  equivalent  to  the  experimental  mobility)  to  set  the 
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parameter  (/Jr#)w1n  the  asymptotic  low-energy  cross-section  based  on  the  (r_&) 
potential. 

C  3.2  Experimental  Mobility. 

The  mobility  in  nitrogen  of  various  ions  has  been  reported  by,  e  g.,  Brata 
[Ref.  11]  and  by  Mitchell  and  Ridler  [Ref.  58],  A  graph  of  mobility  versus  ion  mass 
number  (erroneously  attributed  to  Ref.  1 1 )  is  shown  in  Fig.  1 0. 1 7  of  Hasted 
[Ref.  28];  basically  the  same  graph  (but  slightly  edited),  shown  in  Fig.  3.20  of 
Brown  [Ref.  12]  (and  also  in  Fig.  7m-87  of  Grey  [Ref.  26])  and  correctly  attributed 
to  Ref.  58,  is  reproduced  here  as  Fig.  17.  A  similar  graph  does  appear  in  Fig.  4  of 
Ref  1 1 ,  but  Ba+  does  not  appear  in  it.  The  only  comment  about  Ba+  made  in  Ref  58 
is  reproduced  below: 

”...  those  [mobilities  of  various  ions  in  nitrogen]  previously 
examined  [by  Tyndall,  Powell,  and  Brata] ...  including  Ba+  at  2.23 
obtained  by  Powell  and  Brata  previously  unpublished  ... ." 


Fig.  3.  20.  Mobility  in  nitrogen  of  various  ions  as  a 

function  of  mass  at  1  atmosphere  pressure. 

J.  H.  Mitchell  and  K.  E.  W.  Ridler,  Proc.  Roy.  Soc. 
(London)  A146,  911  (1934). 


Figure  17.  Reduced  mobility  in  N2  of  various  ions  as  a  function 
of  their  mass  number  (from  Brown  [Ref.  12,  p.  62]). 


The  temperature  is  not  stated  in  Refs.  28,  12,  26,  or  58;  Ref  1 1  variously 
refers  to  20  *C,  17  *C,  and  NTP.  We  will  assume  the  temperature  to  be  room 
temperature,  taken  as  T~  293  K. 


From  Eqs.  (A2b),  (A9f),  and  (A  13)  -  or  from  Eq.  (29a)  -  we  see  that,  for  the 
polarization-potential  interaction,  the  mobility  is  inversely  proportional  to  the 
number  density  of  the  gas,  or  equivalently  proportional  to  the  ratio  of  the 
temperature  to  the  pressure.  Hence,  to  adjust  the  experimental  value  of  Ba+ 
mobility  to  STP  conditions,  we  write 

'exp  As  TP 


=  223 


273  1 


293  I 

=  208  (cm/sec)/(vo!t/cm) 


(C5) 


C.3.3  Theoretical  Mobility. 

C.3.3. 1  Polarization  Potential.  From  Eq.  (A26b)  and  Table  27  we  have 

XS1P (Ba+ , N2 )| ^  =  2.16  (cm/sec) /(vo It/cm),  (C6) 

as  given  in  Table  28.  That  the  value  in  Eq.  (C6)  differs  from  that  in  Eq.  (C5)  by  less 
than  four  percent  is,  perhaps,  remarkable.  The  essential  agreement  could  probably 
justify  regarding  the  mobility  of  Ba+  In  N2  (for  the  temperature  in  the  experiment) 
as  being  properly  based  on  the  polarization-potential  interaction. 

C.3.3.2  Effective  (r'6)-Potential.  We  now  determine  an  effective  value  of  the 
parameter  (/Jr,)1*  such  that 

JtSTP(Ba-,N,)|^)<ol  =.  .  (C7) 

(From  a  practical  viewpoint  it  matters  little  whether  we  use  the  theoretical  value 
based  on  the  polarization-potential  interaction  or  the  experimental  value.  We 
choose  the  theoretical  value  to  simplify  some  subsequent  formulas.)  Thus,  from 
Eq.  (A  19),  we  have 

(^‘•.eWW-*,-4.94xlOr“a76)l« 

=  6.55x1  O'16  cm2  eVu3 


and  from  Eqs.  (A28b)  and  (C8), 


(C8) 


J£STP(Ba+,N2)| 


6.82x!  O'15 


Effir-Mx*  (23.26)1/2x  6.55x1  O'16 

=  2.1 6  (cm/ sec) /(volt/ cm). 


(C9) 


C.4  (Ba++02)-C0LLISI0NS.  MOBILITY. 

We  follow  the  procedure  used  for  (Ba++0)-collisions  but  limit  our 
considerations  to  mobility.  Thus,  from  Eq.  (A  19)  we  have 


(/*&  ~  y/3 
K/m  4  o  'til. 


=  4.94x1 0-,6(1.59r 
=  6. 23x1 0‘16  cm2  eVvl 


(CIO) 


and  from  Eqs.  (28b)  and  (CIO)  we  have 


XSIP(Ba*,Oj)|t„(r,)^)t « 


6.82x1  O' 


,-ts 


(25.95)l/2x  6.23x1  O'16 
=  215(cm/sec)/(volt/cm). 


(Cl  1) 
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APPENDIX  D 


AMBIPOLAR  DIFFUSION  COEFFICIENT  -  A  LIMITED  COMPARISON  OF 
EXPERIMENTAL  AND  THEORETICAL  RESULTS  FOR  SOME  BARIUM  RELEASES 

D.l  THEORETICAL. 

The  ambipolar  diffusion  coefficient  is  usually  defined  (cf.,  e.g.,  Dalgarno 
[Ref.  8,  p.  660];  McDaniel  [Ref.  55,  p.  5 1 3];  Fu  [Ref.  2 1  ])  as 

D,X.*D,X,  (D|) 

'j  t«/V# 

where  the  ion  and  electron  mobilities  X,  and  Xt  are  both  positive  numbers. 

With  use  of  the  (Einstein)  expressions  relating  mobility  and  diffusion 
coefficient  (cf.  Eq.  (A2a)), 


Eq.  (Dl)  becomes 
4- 


0 kr,/e)X, 

(D2a) 

(/ kTJe)X9t 

(D2b) 

DtDSr^T) 

DJ.+DJ,  ■ 

(D3) 

Equation  (D3)  is  the  form  in  which  the  ambipolar  diffusion  coefficient  appears, 
e.g.,  in  Holway  [Ref.  33,  Eq.  ( 1 8)]. 

If  one  combines  Eqs.  (Dl )  and  (D2)  and  then  assumes  Xe  »X,  due  to  the 
electron's  small  mass  and  its  concomitant  larger  contribution  to  the  conductivity, 
one  has 

(D4a) 

vD^WTjr,),  (D4b) 

where  Eq.  (D4b)  is  (or  essentially  is)  the  form  for  the  ambipolar  diffusion 
coefficient  presented  in  some  atmospheric  physics  texts,  e.g.,  Banks  and  Kockarts 
[Ref.  7,  Pt.  B,  p.  164],  Ratcliffe  [Ref.  63,  p.  108],  and  Bauer  [Ref.  9,  p.  97], 

D.2  MODEL  ATMOSPHERE  PROPERTIES. 

D.2. 1  Neutral  Particles. 


Later  we  will  need  selected  neutral-particle  densities  for  two  model 
atmospheres  presented  in  Tables  30  and  31. 
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The  ClRA-1965  Model-6  1800-hr  [Ref.  15]  atmosphere  (cf.  Table  3 1 )  appears 
to  be  the  model  atmosphere  used  by  at  least  some  members  (see,  e  g.,  Ref.  46,  p. 
23)  of  the  SECEDE  community  as  being  appropriate  for  the  solar-flux  conditions 
during  the  second  half  of  January  1971  when  the  SECEDE-II  events  occurred. 

Table  30.  Selected  properties  for  CIRA-65,  Model-5,  0600-hr. 


cm'3 

A  km 

T,  K 

[0] 

[NJ 

[02] 

190 

823 

4.25  E+09 

4.83  E+09 

5.45  E+08 

200 

852 

3.33  E+09 

3.32  E+09 

3.45  E+08 

210 

877 

2.63  E+09 

2.19  E+09 

2.23  E+08 

220 

898 

2.1 1  E+09 

1.51  E+09 

1.46  E+08 

cm-3 

A  km 

/?(tC' al) 

f[0] 

rtoj 

190 

9.62  £+09 

0.442 

0.502 

0056 

200 

6.99  E+09 

0.476 

0.475 

0.049 

210 

5.04  E+09 

0522 

0.434 

0.044 

220 

3.77  E+09 

0.560 

0.401 

0.039 

Table  31. 

Selected  properties  for  CIRA-65,  Model-6,  1800-hr. 

A  km 

T,  K 

cm'3 

[0] 

INJ 

[02] 

140 

662 

1.98  E+10 

6.05  E+10 

9.48  E+09 

150 

779 

1.36  E+10 

3.49  E+10 

5. 16  E+09 

160 

878 

9.36  E+09 

2.02  E+10 

2.81  E+09 

180 

1029 

5.49  E+09 

8.92  E+09 

1.13  E+09 

185 

1058 

4.93  E+09 

7.54  E+09 

9.39  E+08 

200 

1 135 

3.59  E+09 

4.55  E+09 

5.33  E+08 

220 

1210 

2.49  E+09 

2.52  E+09 

2.74  E+08 

240 

1264 

1.79  E+09 

1.47  E+09 

1.49  E+08 

250 

1286 

i.  54  E+09 

1. 14  E+09 

1.12  E+08 

260 

1304 

1.33  E+09 

8.85  E+08 

8.37  E+07 

280 

1333 

9  96  E+08 

5.46  E+08 

4.84  E+07 

300 

1355 

7.58  E+08 

3  43  E+08 

2.85  E+07 
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Table  31,  Selected  properties  for  CIRA-65,  Model-6,  1800-hr.  (Cont'd) 


/?,  km 

cm-3 
n(  total) 

f[0l 

f[Nj 

f[02] 

140 

8.98  E+ 10 

0.220 

0.674 

0.106 

150 

5.37  E+  10 

0.253 

0.650 

0.096 

160 

3.23  E+10 

0.289 

0.624 

0.087 

180 

1.56  E+ 10 

0.353 

0.574 

0.073 

185 

1.34  E+ 10 

0.368 

0.562 

0.070 

200 

8.67  E+09 

0.413 

0.525 

0.062 

220 

5.28  E+09 

0.471 

0.477 

0.052 

240 

3.41  E+09 

0.526 

0.430 

0.044 

250 

2.79  E+09 

0.552 

0.408 

0.040 

260 

2.29  E+09 

0.578 

0.386 

0.036 

280 

1.59  E+09 

0.626 

0.343 

0.031 

300 

1.13  E+09 

0.672 

0.303 

0.025 

D.2.2  Electron-  to  Ion-Temperature  Ratio. 

Unson  et  al.,  guided  by  results  in  Wand  [Ref.  67],  have  used  the  following 
prescription  for  the  ratio  of  electron  temperature  to  ion  temperature: 

T  [l.0+^^0.6,  1 30  £  /?(km)  <l  2 1 0 

-=  80  (D5) 

'  1.6--~^|— 0.4,  210s/tfkm)<300. 

Equation  (D5)  appears  to  be  an  approximate  fit  to  the  ratio  TtITi  (for  a  period 
spanning  an  evening  twilight)  appearing  in  Fig.  3c  of  Ref.  67,  reproduced  here  as 
Fig.  1 8. 

Figure  32  in  Bauer  [Ref.  9],  showing  Te  and  7]  for  different  solar-cycle 
conditions,  is  reproduced  here  as  Fig.  19;  the  time  of  day  is  not  specified,  but  it  is 
presumably  daytime  (noon?).  We  have  collated  in  Table  32  values  of  Te  and  T,  (and 
the  ratio  Tt/Tt)  at  200-km  altitude  read  from  Fig.  19. 
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Height  ,  km 


Table  32.  Temperatures  at  200-km  altitude  per  Bauer  [Ref  9,  Fig.  32] 


A  (Solar  Max) 

B(lntermediate) 

C (Solar  Min,  Mag.  Act.) 

T< 

!  125 

1450 

1950 

T, 

700 

775 

875 

TJT, 

1.61 

1.87 

2.23 

Electron  density  (ne) ,  cm‘3 

10s  10®  10s  10® 


Fig.  3c. 


Figure  1 8.  Electron  and  ion  temperatures  versus  altitude  for  1 730- 1 838 
AST,  10  August  1966  (from  Wand  [Ref.  67]). 


Fig.  32.  Electron  and  ion  temperatures  in  the  terrestrial  ionosphere  for  different 
conditions.  Also  shown  are  the  corresponding  electron  density  profiles.  Curves  A 
reflect  a  high  peak  electron  density  (e.g.,  solar  maximum)  and  high  atomic  oxygen 
concentration  at  the  lower  boundary;  curves  C  reflects  a  low  peak  electron  density 
(e.g.,  solar  minimum)  and  an  atomic  oxygen  concentration  reduced  by  50  %  while 
curves  B  reflects  an«yntermediate  condition.  Similarly,  curves  C  are  also  represent¬ 
ative  of  a  situation  occurring  during  magnetic  storms.  (After  Herman  and  Chandra) 


Figure  19.  Electron  and  ion  temperatures  versus  altitude 
for  different  solar-flux  conditions  (from 
Bauer  [Ref.  9]). 


The  smoothed  solar  flux  at  10.7-cm  wavelength  is  shown,  for  the  period 
from  1958  to  1968,  in  Fig.  6  of  Jacchia  [Ref.  36],  reproduced  here  as  Fig.  20.  The 
general  trend  of  the  same  flux  (/To?)  is  shown  in  Fig.  I  of  Hedin  et  ai.  [Ref  30], 
reproduced  here  as  Fig.  21.  From  Figs.  20  and  21,  one  would  judge  Fm  *1 10  at  the 
date  (10  Aug  66)  corresponding  to  the  data  in  Fig.  18.  (It  would  be  of  interest  to 
compare  this  value  with  the  daily  values  for  August  1966,  but  those  data  are  not 
readily  available.  Solar  data  for  an  approximate  period  including  1966  through 
September  1970  are  not  in  the  J.  Geophys.  Res.  but  were  published  in  "Solar- 
Geophysical  Data,"  U.  S.  Dept,  of  Commerce,  Boulder,  Colorado,  80302.  Data  for 
October  1970  are  in  J.  Geophys.  Res.  7£,  1097  (1  Feb  71),  with  subsequent  data 
following  regularly.) 

Solar  flux  data  for  January  1971  are  given  in  Table  33.  The  average  for  the 
last  half  of  January  1971, 

i  / -3 1  Jon  71 

T7-  =  173.3,  (D6) 

1  D  /-16J«i7l 
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is  somewhat  higher  than  the  average  for  the  month,  162.6,  and  considerably  higher 
than  the  value  (»152)  shown  as  the  "general  trend  of  fW7“  in  Fig.  1  of  Ref  30  (our 
Fig.  21). 


Figure  20.  Mean  10.7-cm  solar  flux  index  for  the  period  1958-1968  (from 
Jacchia  [Ref.  36]). 


YEAR  11900*) 

Fig.  I .  MSIS  model  data  coverage  in  altitude  and  year-day  of  year  coordinates.  Also  shown,  as  the  dashed  line,  is  the  gen¬ 
eral  trend  in  the  F„,  inde*  during  this  period. 

Figure  21.  Mean  10.7-cm  solar  flux  index  for  the  period  1966-1975  (from 
Hedin  et  al.  [Ref.  30]). 

Generally  speaking,  the  data  in  Fig.  18  are  for  a  relatively  low  value  of  fw7-, 
at  higher  values  of  F07  one  would  probably  expect  lower  values  of  the  ratio  Te/T,  , 
based  on  the  trend  from  Fig.  18  (and  Table  32). 
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Table  33.  Observed  solar  flux  at  2800  MHz  during  January  1971  [Ref  43] 


Day 

^107 

Day 

^10.7 

Day 

^07 

1 

135.0 

1 1 

154.1 

21 

184  6 

n 

4 

139  4 

12 

155  6 

22 

186  4 

3 

139.3 

13 

153.2 

23 

188.8 

4 

145.0 

14 

159  2 

24 

182  7 

5 

151.1 

15 

158.9 

25 

174.0 

6 

151.9 

16 

161.6 

26 

170.2 

/ 

155.3 

17 

160.4 

27 

171  9 

8 

154.2 

18 

165.6 

28 

1713 

9 

158  1 

19 

1713 

29 

165  7 

10 

157.3 

20 

174.8 

30 

163.0 

31 

175.2 

D.3  COMPARISONS. 

D  3  I  Kivel  [Ref .39]. 

For  Dogwood  (h  -  191  km  (Ref  62]),  the  neutral  cloud  grows  with  a 
diffusion  coefficient  D  *  0. 1  km2/sec.  The  ambipolar  diffusion  coefficient  for  the 
ion  cloud  along  the  Earth's  magnetic  field  is  £},  *  0.2  km2/sec. 

Comments.  The  value  0  2  km2/sec  for  £],  is  larger  than  the  0  1 36  km2/sec 
and  0. 1 1  km2/sec  also  reported  for  Dogwood  respectively  by  Lumenelio  et  ai  [Ref 
50]  and  by  Boquist  et  al.  [Ref.  10]  in  following  Sections  D  3.2  and  D  3  3 

To  make  predictions  to  compare  with  the  quoted  Dogwood  (and  Apple) 
ambipolar  diffusion  coefficients,  we  need  an  atmosphere.  For  convenience  and 
simplicity,  we  use  Figs.  20  and  21  to  read  respectively  *145  and  f]07  *142  [m 
units  of  I  O'22  W/(m2  cyc/sec)]  for  May  1968.  ( Apple  occurred  on  2  May  68; 
Dogwood,  1 3  May  68.)  These  values  of  ,S07  lead  us  to  use  Model-5  (for  which 
/To7  =  150)  in  CIRA-1965  [Ref  15].  For  0600  hrs,  we  have  atmospheric  data  in 
Table  30;  results  for  diffusion  coefficients  are  given  in  Table  34  The  ordinarv 
diffusion  coefficient  D,  for  Ba+  was  evaluated  from  Eq.  (A33),  the  ambipolar 
diffusion  coefficient  Dn  from  Eq.  (D4b),  and  Tef  Tt  from  Eq.  (D5).  The  value  0. 1 24 
km2/sec  for  the  effective  (i — €>)-potential  at  191-km  altitude  in  Table  34  does  not 
agree  well  with  the  value  0.2  km2/sec  quoted  by  Kivel,  but  does  agree  well  with 
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rne  values  0  136  and  0  1 1  quoted  respectively  Dy  Lumenello  et  al  and  Dv  Bouuisr  er 
a).  ^  Sections  D.3.2  and  D.3.3) 


Table  34  Ba+  diffusion  coefficients  for  CIRA-65,  Model-5,  0600-hr 


h 

Dt ,  km2/ sec 

L-\\  ,  Kfi  1  / 

'sec 

km 

Pol. 

Ef  f.  (r-6) 

1  *TJT, 

Pol. 

Eff  (r*) 

LB-7! 

LB-77 

190 

5  84  E-02 

486  E-02 

2  45 

0  143 

0  1 19 

0  114 

1  Of', 

1  L-  V 

191 

2.46 

0.149 

0.124 

0.1  19 

0  123 

194 

2.48 

0.167 

0  139 

0  !3~ 

iV  t  ar , 

196 

2.49 

0.181 

0. 1 50 

0  >47 

0  1 5  1 

200 

8  43  E-02 

6  98  E-02 

2.52 

0212 

0  176 

0  167 

u  1  ~~ 

201 

253 

0221 

0  183 

0  174 

0  1 34 

204 

2.55 

0.250 

0.207 

0  195 

0  207 

205 

2.56 

0.260 

0216 

0  203 

0  2 1 6 

207 

2.58 

0.283 

0.234 

0.219 

a  273 

210 

1.23  E-01 

1.01  E-01 

2.60 

0.320 

0.264 

0.246 

0  263 

215 

2  58 

0  375 

0.309 

0  286 

0  306 

220 

1.72  E-01 

1.41  E-01 

2.56 

0.440 

0.361 

0.332 

0  35™ 

11) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

/  r» 

•  0  ' 

Eq.  (A33) 

Eq.  (A33) 

Eq.  (D5) 

Col.  2 

Col.  3 

Eq.  <D7» 

c.t  /r,7 , 
um  ■  ^  ' 

X 

X 

(D8a) 

(B3  d 1 

Col. 4 

Col.  4 

(D5) 

1 05 1 

H 

U 

tl 

ii 

(D4b) 

(D4b) 

(DlOa) 

(D 1  ob) 

D.3.2  Lumenello,  Davis,  and  Freedman  [Ref.  50], 

Measured  diffusion  coefficients  of  the  ion  cloud  parallel  to  the  field  were 
0.238  km2/sec  for  Apple  ( /?  =  1 96  km)  and  0  1 36  km2/sec  for  Dogwood  n  88  km, 
or  191  km  per  Table  I  in  Ref.  62). 

Comments.  For  Apple,  the  interpolated  value  0.150  km2/sec  for  the 
effective  (r-$)-potential  at  196  km  in  Table  34  does  not  agree  well  with  the  value 
0.238  km2/sec  quoted  by  Lumenello  et  al.  [Ref.  50] 
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For  Dogwood ,  the  value  0  136  km2/sec  quoted  by  Lumenello  et  a!  (Ref  50] 
is  (a)  smaller  than  the  0.2  km2/sec  reported  by  Kivel  (Section  D.3. 1 ),  (b)  larger 
than  the  0.1 1  km2/sec  reported  by  Boquist  et  al.  (Section  D.3.3),  and  (c)  slightly 
larger  than  the  0. 1 24  km2/sec  computed  for  the  effective  (r_6)-potentia!  at  191- 
km  altitude  in  Table  34 

Although  the  report  by  Lumenello  et  al.  [Ref.  50]  is  an  EG&G  report,  it  is  not 
referenced  in  the  later  EG&G  report  by  Boquist  et  al.  [Ref.  10]. 

D.3.3  Boquist,  Overbye,  Klesling,  and  Eves  [Ref.  10]. 

These  authors  report  a  value  of  0. 1 1  km2/sec  for  D\\  for  both  Apple  and 
Dogwood,  each  occurring  during  morning  twilight. 

Comments.  For  Dogwood,  the  value  0. 1 1  km2/sec  is  (a)  less  than  the  0  2 
km2/sec  reported  by  Kivel  (Section  D.3  1 ),  (b)  less  than  the  0  1 36  km2/sec 
reported  by  Lumenello  et  al.  (Section  D.3.2),  and  (c)  less  than  the  computed  values 
(Table  34)  for  D, ,  of  0. 1 24  km2/sec  and  0. 1 49  km2/sec  respectively  for  the 
effective  (r~6)-potential  and  polarization  potential. 

For  Apple ,  the  value  0.1 1  km2/sec  is  (a)  less  than  the  value  0.238  km2/sec 
reported  by  Lumenello  et  al.  (Section  D.3.2)  and  (b)  the  computed  values  (Table 
34)  for  Dt  of  0.150  km2/sec  and  0.181  km2/sec  respectively  for  the  effective 
(r6)-potential  and  polarization  potential. 

D  3.4  Fu,  Marram,  Ponder,  and  Breedlove  [Ref.  22]. 

The  data  in  Table  35  have  been  collated  from  Tables  l ,  2,  and  8  of  Ref  22 

Comments.  Why  is  D^  ( Titmouse )  so  small?  Why  is  C\\  ( Titmouse )  < 

D*i Sapsucker ),  since  hi  Sapsucker)  <  hi  Titmouse)  ? 

Why  is  Du  {Mulberry)  so  small?  Why  is  {Mulberry)  <  D\  iKumcjuat),  since 
hiKumquat)  <  hi  Mulberry)  ? 

D.3. 5  Linson  and  Baron  [Ref.  45]. 

This  short  paper  includes  the  statements: 

"...  Beyond  20  sec,  the  ion  cloud  expands  by  ambipolar  diffusion 
parallel  to  the  magnetic  field.  An  ambipolar  diffusion  coefficient  of 
D  ■  0.1  km2/sec  at  185-km  altitude  is  consistent  with  an  ion-neutr3l 
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TaOie  35.  BIRDSEED  and  PRE-SESEDE  data  from  Fu  et  a!  [Ref  22] 


Pm 

Date 

Name 

Altitude 

km 

Time 

UT 

A 

km2 /sec 

*152 

26  May  70 

BIBDSEED 
Sapsucker b 

204+3 

0552.25 

0.3 

*152 

6  June  70 

Titmouse D 

2!0±5 

0555:25 

02 

*145 

5  Sep  69 

PRErSECEDt 

KumpuaP- 

194 

0.26 

*145 

6  Sep  69 

Lime c 

201  5 

0  29 

*145 

1 3  Sep  69 

Mulberry c 

195.8 

0  13 

3  As  read  by  the  writer  from  the  "general  trend  F[01"  in  Fig.  I  of 
Ref  30  (see  our  Fig.  21) 

b  The  results  for  these  events  are  the  primary  purpose  of  Ref  22. 
c  The  results  for  these  events  are  included  in  Ref.  22  for  comparison 
purposes. 


collision  time  of 

t,  =  9x10V;1 

where  Nn  is  the  neutral  density  in  cm-3,  with  an  electron-plus-ion 
temperature  of  2300  *K." 

"Excellent  agreement  between  the  radar  data  and  theoretical 
expressions  is  obtained  for  the  first  10  minutes,  indicating  that  the 
value  for  (0. 1  km2/sec)  is  appropriate  at  early  times  for  these  two 
[  Spruce  and  Olive]  releases." 

Comments.  The  quoted  value,  0.  l  km2/sec  at  185-km  altitude,  appears  to  be 
in  contradiction  with  a  later  value,  0. 1 5  km2/sec  at  1 85-km  altitude,  quoted  by 
Linson  [Ref.  44,  p.  5]. 

it  is  desirable  to  understand  why  185-km  altitude  is  selected  for  a  quoted 
{\  when  explicit  results  are  given  (on  the  next  page  of  Ref.  45)  for  Spruce  and 
Olive  (at  altitudes  of  1 87  3  and  193.4  km,  per  Ref.  62,  Table  l,p  3).  Could  it  be 
that  the  text  should  have  referred  to  Olive  Pit  (at  185  km)  instead  of  Olivel 


We  should  like  to  verify  the  0  1  km2/sec  by  using  formulas  likely  to  have 
been  used  by  Ref.  45.  Thus,  by  combining  Eqs.  (Al )  and  (D4b)  we  have 


^  ~  m,  v,  *  m,  " 

(D7) 

from  the  Ref.  45  text, 

r,=  9xl0V/V„ 

(D8a) 

T,  +  Te  =  2300  K, 

(D8b) 

and  for  185-km  altitude  in  Table  31, 

/V„  =  1. 34  xlO10  cm'3 

(D8c) 

Hence, 

l\  =  0. 093  km2 /sec, 

(D9) 

which  is  only  seven  percent  smaller  than  the  quoted  0  I  km2/sec. 


If,  instead  of  using  Eq.  (D8b)  for  the  temperatures,  we  use  for  Tf  the  neutral 
gas  temperature  of  1 054  K  at  185-km  altitude  in  Table  31  and  Eq.  (D5)  giving  a 
value  of  141  for  the  ratio  Te!  T, ,  we  then  obtain  Z}  =  0.1 04x1 010  km2/ sec 


In  Section  B.  I  we  disagree  with  the  use  of  Eq.  (B3)  («Eq.  (D8a))  for  the 
collision  time  -  though  the  disagreement  is  more  in  principle  than  in  practice, 
specially  with  consideration  of  the  improved  formula  in  Ref.  46.  The  proper 
formula  is  given  by  Eq.  (B7)  (or  by  Eq.  (B3b),  explicitly  for  the  polarization 
potential),  with  (interpolated)  results  from  Table  31  giving  r,  =0.747 and  r,  =0  5% 
respectively  for  collisions  based  on  the  polarization  and  effective  (r~6)- 
potentials.  Evaluation  of  Eq.  (D8a)  («Eq.(B3))  with  use  of  Eq.  (D8c)  gives 
t,  =0.67  sec,  a  value  intermediate  to  what  we  would  propose  by  factors  of  0  90  and 
I  1 3  for  the  polarization  and  effective  (r_6)-potentials. 

Use  of  the  ordinary  diffusion  coefficient  in  Table  36  for  185  km  and  use  of 
Eq.  (D5)  for  the  temperature  ratio  give  values  for  4  of  0.12  and  0.092  km2/sec 
respectively  for  the  polarization  and  effective  (r_6)-potentials. 

D.3.6  Linson  (Ref.  44], 

This  paper  [p.  53  states : 

"...  appropriate  ambfpolar  diffusion  coefficient  A =0.04, 0.15,  and 0.7 
for  the  three  altitudes  (150,  1 85,  and  250  km].“ 


values. 


It  is  not  clear  whether  these  are  theoretical  or  experimental 


Instead  of  0. 1 5  km2/sec  as  quoted  here  for  185  km,  Ref.  45  quotes  0.  i 
km2/sec. 


Our  computed  values  are  given  in  Columns  5  and  6  of  Table  36.  As  a  set,  the 
better  agreement  with  the  Ref.  44  values  is  for  those  obtained  with  the 
polarization  potential. 


In  Table  36  we  also  record  the  values  of  A  computed  from  Eqs.  (D7)  and 
(D8a)  (...  Eq.  (B3)],  i.e.,  the  Ref.  45  prescription,  and  with  TJT,  given  by  Eq.  (D5). 

A*5. 45x10s -7  !  +  -£  ,  km2/sec.  (DIGa) 

Nn  L  T< 

Similarly,  we  record  in  Table  36  the  values  of  l\  computed  from  Eqs  (D7) 
and  (B3a),  i.e.,  the  Ref.  46  prescription,  and  with  Te/Tt  given  by  Eq.  (D5) 


>5  20xl05  [n2]  +  (o,]  +  0.8[0]  U  T, 


km2/ sec 


(D I  Ob) 


We  see  that  the  results  from  Eq.  (DlOb)  agree  slightly  better  with  the  Ref  44 
values  than  do  those  from  Eq.  (0 10a). 


D.3.7  Summary  Table  and  Comments. 


A  summary  of  all  the  comparisons  is  provided  in  Table  37. 

In  the  Llnson-Baron  [Ref.  45]  and  Llnson-Baxter  [Ref.  46]  formulas  for  the 
coefficient  for  ambipolar  diffusion  along  the  magnetic  field, 

-  *r,[t  rA 

mi  [ +  T,\r‘”’  (11) 

the  expression  for  the  ion-neutral  collision  time  r,„,  given  by  Eq.  (B3)  [Ref  45]  and 
Eq.  (B3a)  [Ref.  46]  -  which  is  nominally  based  on  the  experimental  mobility  of  Ba+ 
in  N2  -  improperly  accounts  for  0  and  O2  if  a  polarization-potential  interaction  is 
assumed.  Correction  of  this  relatively  small  error  by  properly  accounting  for  tne 
polarizabilities  of  O2  (which  actually  does  not  matter)  and  0  (which  does  matter) 
does  not  significantly  change  the  overall  comparison  between  theoretical  values  of 
the  ambipolar  diffusion  coefficient  [z\  ^  ^  l  based  on  the  polarization  potential 
and  (sometimes  ambiguous  or  uncertain)  experimental  values  [fy  exp )  derived  from 
12  Ba-releases.  Values  of  the  ambipolar  diffusion  coefficient  [on  lfi  ] 

computed  from  an  effective  (rS)-potentlal  -  which  are  about  79%  to  85%  of 
A.  poi  p<*  for  ionospheric  ion  temperatures  ranging  from  1 135  K  to  750  K  -  are 
(accidentally)  very  nearly  the  same  as  those  obtained  with  the  Ref.  46  formula  for 
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Table  36.  Ba*  diffusion  coefficients  for  CIRA-65,  Model-6,  1800-nr 


/)  D, ,  km2  /sec  _ AJ<m2/sec 


km 

Pol. 

Eff.  (r6) 

UTJT, 

Pol. 

Eff  (r6) 

ReT,  44  Ref,  45  Ref  46 

140 

4. 1 3  E-03 

3,56  E-03 

2.08 

0.009 

0.007 

0  008 

0  008 

ISO 

8.28  E-03 

6.95  E-03 

2.15 

0.018 

0.015 

0.04 

0.017 

0.017 

160 

1.58  E-02 

1  30  E-02 

2.23 

0.035 

0.029 

0033 

0  033 

180 

3.97  E-02 

3. 19  E-02 

2.38 

0.095 

0.076 

0.085 

0.088 

185 

4.80  E-02 

3.83  E-02 

2.41 

0.116 

0.092 

0.15 

0  104 

0  107 

200 

8. 17  E-02 

6.44  E-02 

2.53 

0.207 

0.163 

0.180 

0  188 

220 

1.48  E-0! 

1.16  E-01 

2.56 

0.379 

0.296 

0.319 

0.337 

240 

2.48  E-01 

1  92  E-01 

2.47 

0.612 

0.474 

0  499 

0  532 

250 

3.13E-01 

2.42  E-01 

2.42 

0.757 

0.585 

0.7 

0.607 

0652 

260 

3.94  E-01 

3  03  E-01 

238 

0.937 

0.722 

0  738 

0  794 

280 

5.98  E-01 

4.59  E-01 

2.29 

1.37 

1.05 

1.05 

1.14 

300 

8.83  E-01 

6.76  E-01 

2.20 

1.94 

1.49 

1.44 

1  59 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

Eq.  (A33) 

Eq.  (A33) 

Eq.  (D5) 

Col.  2 

Col.  3 

Eq.  (D7 

)Eq.(D7) 

X 

X 

(D8a) 

(B3a.) 

Col. 4 

Col.  4 

(D5) 

(D5) 

If 

II 

If 

If 

(D4b) 

(D4b) 

(D 1 0a) 

(D10D) 

one  of  the  model  atmospheres  used.  An  overall  comparison  of  Dn  Efr  with 

A,  *xp.  is  about  the  same  as  for  C\\,  ^  p*,  so  that  a  strong  preference  cannot  be 
established  on  this  basis,  unless  one  should  have  reasons  (unknown  to  the  writer) 
to  select  or  exclude  certain  experimental  data. 

It  would  be  desirable  to  determine  whether  or  not  there  exists  a  definitive 
set  of  experimental  values  for  all  the  Ba-releases,  including  the  six  evening 
twilight  Ba-releases  for  Project  STRESS  ( Anne,  early  December  1976;  Betty, 
Carolyn,  Dianne,  Ester,  and  Fern  in  late  February  and  March  1 977  [Ref.  47 
(substantially  in  Ref.  48)].  Reference  53,  while  providing  an  overview  of  Project 
STRESS,  gives  no  information  regarding  early-time  ambipolar  diffusion 
coefficients,  since  the  interest  was  in  late-time  striations,  A  similar  comment 
applies  to  Project  PLACES  [Ref.  54], 
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Table  37  Summary  of  measured  and  predicted  Ba+  diffusion  coefficients  for  barium  releases 
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APPENDIX  E 


RELATIVE  SPEED  FOR  MONOENERGETIC  BEAM  IN  A  MAXWELLIAN  GAS 


E.  1  APPROXIMATE  FORMULAS. 


In  developing  the  energy  loss  rate  for  a  particle  traversing  a  Maxwellian  gas 
(Section  3,3. 1 ),  we  were  not  initially  aware  of  a  literature  formula  for  the  mean 
relative  speed  between  a  monoenergetic  beam  and  the  particles  of  a  traversed 
Maxwellian  gas.  We  expediently  started  with  the  formula  for  the  relative  speed 
for  the  particles  of  two  separate  Maxwellian  speed  distributions  [Ref  6,  Ref  7,  Pt 
A,  p.  1 89J  and  adapted  it  to  the  current  need  by  assuming  that  the  beam  speed  ( vb ) 

was  the  rms  speed  corresponding  to  the  fictitious  temperature  assigned  to  the 
beam  ( fbct),  i.e., 

(E-i) 

Note  that  the  relative  speed  for  two  Maxwellian  gases  (cf.  Eq.  (48b)  in 
Section  3.2.2)  is  just  the  square  root  of  the  sum  of  the  squares  of  the  mean  speeds 
for  the  two  individual  gases.  Thus,  as  the  temperature  of  one  gas  approaches  zero, 
the  relative  speed  tends  toward  the  mean  speed  of  the  second  gas.  Owing  to  this 
fact,  we  now  express  the  relative  speed  of  the  beam  traversing  the  Maxwellian  gas 
(#*)  in  units  of  the  mean  speed  of  the  Maxwellian  gas  Wg)  and  generalize  the 

relation  of  Vb  to  Tbc\  i.e., 

_  [(8/3*)  for  Vb  =>  rms  speed  corresponding  to  7;,ict 
1  for  Vb  =>  mean  speed  corresponding  to  T"a 
R-vb/vg. 

E.2  EXACT  FORMULA. 


(E-2a) 

(E-2P) 

(E-2c) 


Fite  et  al.  (Ref.  20,  Appendix  A.C]  have  derived  the  formula  for  the 
distribution  of  relative  speed  (i/r)  for  a  monoenergetic  particle  of  speed  w  passing 
through  a  Maxwellian  gas, 


fWr,w)  di^  = 


2  A 


3/2 


jf'2 


exp[-^(i/,2+w'2)] 


sinh(2ixr  pw) 
-jTw 


where 


K->  0 


(E-3a) 


E-l 


and  p,  the  inverse  square  of  the  most  probable  speed,  is 


P  =  Mgl2kTg 

The  mean  relative  speed  would,  of  course,  be  defined  as 


_  r°° 

vr  =  vrrwr>w)  dt-v 

'0 


(E-3b) 


(E-4) 


A  formula  for  the  mean  relative  speed  is  presented  in  Ref.  20,  without 
intermediate  steps.  However,  one  of  the  terms  in  the  formula  contains  the  error 
function  with  the  obviously  incorrect  (dimensional)  argument  [pU2-w],  whereas 
elsewhere  in  the  formula  the  (dimensionless)  product  *  a  p'nw  appears  Thus,  one 
would  expect  that  the  argument  of  the  error  function  should  probably  be  r,  a 
surmise  consistent  with  the  equations  in  Appendix  B  of  Ref.  20  and  with  formula 
3.562.5  on  p.  365  of  Ref  25: 


-oo 

I  x2zw(-px2)s\vk\{yx)tix  = 

Jo 


^(2  p+y2) 
Q/p'2 


exp 


AP 


Erf 


2  p 


1/2 


Hence,  one  can  write 


(*fi) 


1/2 


r  J7I/2 

2j(\+2x2)lrUx)+exp(-x2) 


'  A  P2 
(E-5) 

(E-6) 


if  one  introduces  the  mean  speed  of  the  gas  %), 

ll/2 


then 


8  kTg 

ft  n„ 


and 

with 

Thus, 


n 

U  py"2  ={7g/ 2 
R=w!Vg 


. ,1/2 

-x—(U2x2)Ert(x)+exQ(-x2) 


(E-7) 

(E-8a) 

(E-8b) 

(E-8C) 

(E-9) 


it  is  readily  shown,  by  using  the  series  expansion  of  Erf(x)  given  as  Formula 
7.1.6  in  Ref.  1,  that 

-ErfU)— 

*  jp2  (E- 1 0) 

Hence, 

9 bg/Vg  — '  •  (E-11) 

£.3  COMPARISON  OF  APPROXIMATE  AND  EXACT  RESULTS. 


Hence,  the  exact  and  both  of  the  approximate  formulas  give  the  same  result 
at  very  low  beam-speeds  (/?-» 0).  At  very  high  beam-speeds,  the  exact  formula 
for  gtqlVg  gives  R\  the  general  approximate  expression,  Eq.  (E-2a),  gives  a'bnR , 
so  that  the  first  approximate  formula  [Eq.  (E-2b1 )]  gives  the  slightly  smaller  value 

of  (8/3->r)  R  (  =  0.921  R)  and  the  second  approximate  formula  [Eq.  (E-2b2)j  gives 
R  Results  from  evaluating  the  exact  and  two  approximate  formulas  for  R  = 

0(0. 1  )3(  1)10  are  given  in  Table  38.  For  low  beam-speeds  (R  £  1 )  the  first 
approximate-formula  results  are  in  excellent  agreement  with  the  exact  results, 
but  the  second  approximate-formula  results  provide  a  better  overall  fit  to  the 
exact  results. 


E-3 


Table  38.  Comparison  of  two  approximate  methods  with  the  exact  method  of 
computing  relative  speed  for  monoenergetic  beam  in  (ambient) 
Maxwellian  gas. 


Relative  speed ,g^tvg  Relative  speed lgbgIV9 

Approximate  Ea.  Exact  Eq.  Approximate  Eo.  Exact  Ea 


R 

(E-2P1) 

(E-2P2) 

(E-9) 

(E-2b1) 

(E-2P2) 

(E-9) 

0.0 

1.0000 

1 .0000 

1.0000 

1.9 

2.0160 

2.1471 

2.1066 

0.1 

1 .0042 

1 .0050 

1.0042 

2.0 

2,0965 

2.2361 

2.1963 

0.2 

1.0168 

1.0198 

1.0169 

2.1 

2. 1 779 

2.3259 

2.2870 

0.3 

1.0375 

1 .0440 

1.0378 

2.2 

2.2602 

24166 

2.3785 

0.4 

1,0657 

1.0770 

1.0666 

2.3 

2.3431 

2.5080 

2.4707 

05 

1.1010 

1.1 180 

1.1029 

2.4 

2.4268 

2.b.'C: 

2.5636 

0.6 

1 . 1 426 

1.1662 

1  1 462 

2.5 

2.51  10 

2.6926 

2  6571 

0.7 

1.1899 

1.2207 

1.1961 

2.6 

2.5958 

2.7857 

2.7510 

0.8 

1.2423 

1.2806 

1.2518 

2.7 

2.6810 

2.8792 

2.8454 

0.9 

1.2991 

1.3454 

1.3129 

2.8 

2.7667 

2.9732 

2.9402 

1.0 

1.3597 

1.4142 

1.3787 

2.9 

2.8528 

3.0676 

3.0354 

1.1 

1.4238 

1.4866 

1.4487 

3.0 

2.9393 

3.1623 

3  !  309 

1.2 

1 .4907 

1.5620 

1.5224 

40 

3.8185 

4. 1 23 1 

4  0982 

1.3 

1  5603 

1  6401 

1.5993 

5.0 

4.7139 

5.0990 

5  0785 

1.4 

1.6321 

1.7205 

1.6789 

6.0 

5.6176 

6.0828 

6.0654 

1.5 

1.7058 

1.8028 

1.7609 

7.0 

6.5263 

7.071  1 

7.0561 

1.6 

1.7813 

1.8868 

1.8449 

8.0 

7.4381 

8.0623 

8.049! 

1.7 

1.8583 

1.9723 

1.9307 

9.0 

8.3519 

9.0554 

9.0436 

1.8 

1.9365 

2.0591 

2.0180 

10.0 

9.2673 

10.0499 

10.0393 

Relative  speed  is  in  units  of  the  ambient  mean  speed.  R  is  the  ratio  of  beam 
speed  to  ambient  mean  speed. 


APPENDIX  F 


MAXWELLIAN  ENERGY  DISTRIBUTION 


The  (dimensionless)  Maxwellian  differential  energy  distribution  function 
is 

[ 7  1  Ml/2 

f(f)  =  2^exp(f),  (Fla,  lb) 

The  integral  distribution  f(f)  is 

n  ( ) = £>(  f )  o  f  =  -  n  ( )  *  Erf [V#  |  (F2) 


where  ErfCAr)  is  the  error  function, 

Erf(*)-^  J>P<  -t2)  dz\ 


(F3) 


Figure  22  shows  the  differential  and  integral  distribution  functions 


Figure  22.  Differential  [/(*•)]  and  integral  [  F(e)\  Maxwellian  energy 
distributions. 
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